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Abstract 



In this thesis four separate problems in general relativity are considered, divided 
into two separate themes: coordinate conditions and perfect fluid spheres. Regard- 
ing coordinate conditions we present a pedagogical discussion of how the appropriate 
use of coordinate conditions can lead to simplifications in the form of the spacetime 
curvature — such tricks are often helpful when seeking specific exact solutions of the 
Einstein equations. Regarding perfect fluid spheres we present several methods of 
transforming any given perfect fluid sphere into a possibly new perfect fluid sphere. 
This is done in three qualitatively distinct manners: The first set of solution generat- 
ing theorems apply in Schwarzschild curvature coordinates, and are phrased in terms 
of the metric components: they show how to transform one static spherical perfect 
fluid spacetime geometry into another. A second set of solution generating theorems 
extends these ideas to other coordinate systems (such as isotropic, Gaussian polar, 
Buchdahl, Synge, and exponential coordinates), again working directly in terms of the 
metric components. Finally, the solution generating theorems are rephrased in terms 
of the TOV equation and density and pressure profiles. Most of the relevant calcula- 
tions are carried out analytically, though some numerical explorations are also carried 
out. 
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Preface 



This thesis looks at four problems in general relativity: 

• Coordinate conditions in general relativity: A method of looking for nice solu- 
tions of the Einstein equations. 

• Develop several new transformation theorems that map perfect fluid spheres into 
perfect fluid spheres. 

• Develop a systematic way of classifying the set of all perfect fluid spheres. 

• Rephrase all these theorems which apply in Schwarzschild coordinates, directly 
in terms of the pressure profile and density profile. 

General relativity (GR) or general relativity theory (GRT) is a theory of gravitation 
discovered by Albert Einstein. The fundamental physical postulate of GR is that 
the presence of matter causes curvature in the spacetime in which it exists. This 
curvature is taken to be the gravitational field produced by the matter. Einstein's 
field equation gives the mathematical description of how the matter and curvature are 
related. Moreover, once this curvature is given, GR describes how other objects (such 
as planets and light beams) move in this gravitational field via the geodesic equation. 

In addition, general relativity states that clocks run slower in strong gravitational 
fields (or highly accelerated frames), predicting a gravitational redshift. It also pre- 
dicts the existence of gravitational lensing, gravitational waves, gravitomagnetism, the 
Lense-Thirring effect, and relativistic precession of orbiting bodies. 
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The first problem is looking for nice solutions of the Einstein equations. However, 
Before stating Einstein's equation, we need to briefly describe the concept of special 
relativity, otherwise general relativity will be hard to understand. 

Structure of the thesis 

This thesis has been written with the goal of being accessible to people with a basic 
background in general relativity, especially in coordinate conditions in general relativ- 
ity and perfect fluid spheres. 

This thesis is made up of six chapters and three appendices. Two of appendices are 
papers published on work relating to this thesis. One was produced as a collaboration 
with Dr. Matt Visser, Silke Weinfurtner, Celine Cattoen and Tristan Faber. The 
second was produced in collaboration with Dr. Matt Visser and Silke Weinfurtner. At 
the time of writing all papers have been published. 
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Chapter 1 

General introduction 



1.1 Introduction 

General relativity (GR) is the theory of space, time and gravitation formulated by 
Einstein in 1915. General relativity is a beautiful theory, however it is often regarded 
as a very deep and difficult theory. It has been considerably developed since the 
late 1950s. Furthermore, in the mid 1960s, the modern theory of gravitational col- 
lapse, singularities, and black holes has been developed. In this thesis, we present 
the coordinate condition in general relativity, we will discuss the recent advances and 
developments in methods of looking for a nice solutions of Einstein equation. In ad- 
dition, we also present the perfect fluid spheres in general relativity. Understanding 
the concept of perfect fluid spheres are important. They are first approximations to 
attempt at building a realistic model for a general relativistic star. We introduce the 
"algorithmic" techniques that permit one to generate large classess of perfect fluid 
spheres in a purely mechanical way. Furthermore, we will extend these algorithmic 
ideas, by proving several solution-generating theorems of varying levels of complexity. 
We shall then explore the formal properties of these solution-generating theorems and 
then will use these theorems to classify some of the previously known exact solutions. 
Indeed, we develop several new transformation theorems that map perfect fluid spheres 
into perfect fluid spheres. Additionally, we develop a systematic way of classifying the 
set of all perfect fluid spheres. The remaining part of this thesis will be devoted to 
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rephrasing all these theorems, which originally apply in Schwarzschild coordinates, 
directly in terms of the pressure profile and density profile. 

Now we would like to give a mathematically precise formulation of the ideas in 
Einstein's equation. Before stating Einstein's equation, we need to describe in brief 
the concept of special relativity, otherwise general relativity will be hard to understand. 



To start we give an exposition of elementary topics in the geometry of space, and 
the spatial tensorial character of physical laws in prerelativity physics. We formulate 
general relativity and provide a motivation for Einstein's equation, which relates the 
geometry of spacetime to the distribution of matter in the universe. 

1.2.1 General and Special Covariance 

Assume that space has the manifold structure of 1R 3 and the association of points of 
space with elements (x 1 , x 2 ,x 3 ) of M 3 can be accomplished by construction of a rigid 
rectilinear grid of metersticks. We call the coordinates of space derived in this manner 
as Cartesian coordinates [21]. The distance, S, between two points, x and x, defined 
in terms of Cartesian coordinates by 



1.2 Einstein's equation 





This formula is the distance between two points. 

Referred to equation (jl.lj) . the distance between two nearby points is 




(1.2) 



Therefore, the metric of space is given by 





In the Cartesian coordinate basis, we derive 





a,b 
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with h a b = diag (1, 1, 1). This definition of h ab is independent of choice of Cartesian 
coordinate system. 

When the components of the metric in the Cartesian coordinate basis are constants, 
we get 

d a h bc = 0. (1.5) 

The space is the manifold IR 3 which possesses a flat Riemann metric. We are able to 
use the geodesies of the flat metric to construct a Cartesian coordinate system. We use 
the fact that initially parallel geodesies remain parallel because the cuvature vanishes. 

1.3 Special Relativity 
1.3.1 Minkowski space-time 

We introduce a 4- dimensional continuum called space-time in which an event has 
coordinates (t,x,y,z) [TH] . 

Minkowski space-time is defined as a A- dimensional manifold provided with a flat 
metric of signature +2. By the definiton, since the metric is flat, there exists a special 
coordinate system covering the whole manifold in which the metric is diagonal, with 
diagonal elements equal to ±1. For convenience, we prefer to use the convention that 
lower case latin indices run from to 3. The special coordinate system is called a 
Minkowski coordinate system and is written 

(x a ) = 0°, x\x 2 , x 3 ) = (ct, x, y, z). (1.6) 

In special relativity [21], it is assumed that spacetime has the manifold structure of 

The spacetime interval, S, between two events x and x defined by 

S = -(x° - x ) 2 + (x 1 - x 1 ) 2 + (x 2 - x 2 ) 2 + (x 3 - x 3 ) 2 (1.7) 

in units where c = 1. 

From the equation (jl.7j) . we are able to defind the metric of spacetime r] ab by 

3 

ds 2 = Vab dx a dx b (1.8) 

a,b=Q 
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with 7]^ — diag (-1, 1, 1, 1), where x a is any global inertial coordinate system. 

Therefore, the ordinary derivative operator, d a , of the global inertial coordinates 
satisfies 

OaVbc = (1.9) 

The curvature of r\ a \j vanishes. In addition, we can parameterize timelike curves by 
proper time, r, defined by 



T 



-r) ab T a T b y/ 2 dt, (1.10) 



where t = arbitrary parameterizaton of the curve, and T a = the tangent to the curve 
in this parameterization. 

The tangent vector u a to a timelike curve parameterized by r is defined by the 
4-velocity of the curve. The square of any 4-vector is an invariant and so 

u a u a = -c\ (1.11) 

where c = 1 so we get, 

u a u a = -l. (1.12) 
In the absence of external forces, its 4-velocity will satisfy the equation of motion, 

u a d a u b = 0, (1.13) 

where d a is the derivative operator associated with rj ao . In addition, when forces are 
present, the equation (jl.lHj) . u a d a u b is nonzero. Furthermore, all material particles 
have a parameter known as "rest mass", m, which appears as a parameter in the 
equations of motion when forces are present. We can define the energy momentum 
4-vector, p a , of a particle of mass m by 

p a = mu a . (1.14) 

Finally, we can define the energy of particle which is measured by an observer - present 
at the site of the particle - as 

E = -p a v\ (1.15) 
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where v a is the 4-velocity of the observer. 

In special relativity, the energy is the "time component" of the 4- vector, p a . At the rest 
frame, a particle with respect to the observer, equation ()1.15|) reduces to the familiar 
formula E = mc 2 . When the spacetime metric, r} a f, is flat, and the parallel transport 
is path independent, we are able to define the energy of a particle as measured by an 
observer who is not present at the site of the particle and has 4-velocity parallel to 
that of the distant observer. 



1.3.2 The Stress-Energy Tensor 



In special relavity, we define the energy-momentum 4-vector of a particle of mass m 
as in equation (jl.l4|) . 

From the Minkowski line element takes the form 



ds 2 = -dt 2 + dx 2 + dy 2 + dz 2 
We can write this in tensorial form as 

ds 2 = r] ah dx a dx b , 
We take rj a b to denote the Minkowski metric 



Vab 



-1 














+1 














+1 














+1 



diag(-l,+l,+l,+r 



In relativistic units the equation for the proper time satisfies 

dr 2 = -ds 2 . 



(1.16) 



1.17) 



1.18] 



1.19) 



Now we present proper time r relates to coordinate time t for any observer whose 
velocity at time t is v, where 



dx dy dz 



(1.20) 
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So we have, 



dr 2 



-ds 2 = -(-c 2 dt 2 + dx 2 + dy 2 + dz 

2 / , \ 2 



-dt' 



dx V 

dty 



dy 
dt 



-c 2 dt 2 {-l + 



dx 
dt 



dy 
dt 



2 (dz" 2 



c 2 dt 2 |l 



dr 



v 



cr 

2\ V2 



1.21) 



The time-component of the energy-momentum vector does represent the energy of the 
particle 



P 



o 



E 
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2\ -1/2 



1 - 



and E = mc 2r y. 



c \ c z 

The space-components are the components of the three-dimensional momentum 



1.22) 



p = rwyv. 



1.23) 



A perfect fluid is defined to be a continuous distribution of matter with stress 
energy tensor T & of the form 



Tab = PUaUb + p(Vab + U a U b ), 



1.24) 



where u a is the 4-velocity of the fluid, p is the mass-energy density in the rest-frame 
of fluid, and p is the the pressure in the rest-frame of the fluid. 

When there is no external forces, the equation of motion of a perfect fluid is simply 



d a T ah = 0. 



1.25) 
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Consider d a T ab , we can write this in terms of p, p, and u a as 

d a T ab = d a [pu a u b + p(r} ab + u a u b )] 

= (d a p)u a u b + (d a p) (rj ab + u a u b ) 

+ (p + p)(d a u a )u b + (p + p)u a (d a u b ) 
= [(u a d a p) + (p + p)d a u a ]u b 

+ [0 + p)u a d a u b + d a p(r] ab + u a u b )} (1.26) 

For equation ()1.26|) we can project the resulting equation parallel and perpendicular 
to u b , we find: 

[(u a d a p) + (p + p)d a u a ] = 0, (1.27) 



[{p + p)u a d a u b + d a p( Vab + u a u b )} = 0. (1.28) 

;ivistic 1 

(jl.27|) becomes, 



In the non-relativistic limit, when p < p , m m = (1, V), and v ^jj <C |Vp|, equation 



u a d a p + pd a u a = 
<9 t p + ~v ■ V p + p V • 1? = 

<9 t p+V-(p^) = (1.29) 



and equation (|1.28|) becomes 



1.3.3 Relativist ic hydrodynamics 

Ordinary hydrodynamics is dealing with two basic equations |22j : 
The equation of continuity 

d t p + V-{plT) = (1.31) 

This relates to the density and velocity of the fluid. It is equivalent to the conservation 
of mass. 



The Euler equation is the fluid mechanics equivalent to Newton's second law, it 
relates the acceleration of a particle following the flow 

We can write it in term of force density / and mass density p: 

~a=^- + (~v-V)l? = ^- (1.33) 
at p 

These laws have relativistic generalizations which are: 

V a (pv a ) = (1.34) 

A a = V b V b V a (1.35) 

where A a is now a 4-vector field of 4-accelerations. 
The relativistic continuity equation yields 

^ = V-( P7 ^) (1-36) 

We interpret p as proportional to the number density of particles as measured by an 
observer moving with the fluid. Indeed, the 7 factor is corresponding to the fact that 
Lorentz contraction "squashes" in the direction of motion, therefore, as seen by an 
observer moving with respect to the fluid the number density of particles is pj. 
Now for the 4-acceleration 

A* = (r/dt + jF?- V])[ T ^] (1.37) 

and 

A° = (jd t + j[-v- V])[ 7 ] (1.38) 

The standard Newtonian results is reproduced at the low velocity, where 7 — ► 1. 

The special theory of relativity only deals with flat spacetime and the motion of 
objects is usually treated in terms of Lorentz transformations and translations from 
inertial frame to another. The general theory of relativity extends it to deal with 
non-inertial frames, and via Einstein's equations with curved spacetimes as well. We 
need to understand the basic concepts of special relativity before turning to gen- 
eral relativity — as otherwise the mathematical constructions used in general relativity 
would appear rather unmotivated. 
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1.4 Conclusion 



1.4.1 Postulational formulation of special relativity 

There are two sets of postulates which are useful to generalize the general theory ^5] • 

Postulate I. Space and time are represented by a 4- dimensional manifold provided with a 
symmetric affine connection, T a bc , and a metric tensor, g ab , which is satisfied as 
follows: 

(i) g a b is non-singular with signature — h ++; 

(ii) V c g ab = 0; 

(iii) R a bcd = 0. 

The Postulate states that r a f, c is the metric connection and that the metric is flat. 
Postulate II. There exist privileged classes of curves in the manifold singled out as follows: 

(i) ideal clocks travel along timelike curves and measure the parameter r defined by 
dr 2 = — g ab dx a dx b ; 

(ii) free particles travel along timelike geodesies; 

(iii) light rays travel along null geodesies. 

The first part of the second postulate makes physical the distinction between space 
and time in the manifold. In Minkowski coordinates, it distinguishes the coordinate x° 
from the other three as the "time" coordinate. Furthermore, it states that the proper 
time r which a clock measures is in accordance with the clock hypothesis. The rest of 
Postulate II singles out the privileged curves that free particles and light rays travel 
along. 
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1.4.2 The correspondence principle 

Any new theory is consistent with any acceptable earlier theories within their range of 
validity. General relativity must agree on the one hand with special relativity in the 
absence of gravitation and on the other hand with Newtonian gravitational theory in 
the limit of weak gravitational fields and low velocities by comparing with the speed 
of light. This gives rise to a correspondence principle, as in figure (jl.lj) . where arrows 
indicate directions of increased specialization. 



General relativity 




Newtonian theory 
of gravitation 












Special relativity 




Newtonian mechanics 
in the absence of 
gravitation 





Figure 1.1: This structure shows the correspondence principle for general relativity. 

1.4.3 Discussion 

Special relativity is a physical theory published in 1905 by Albert Einstein [2|. The 
theory is a special case of Einstein's theory of relativity where the effects of acceler- 
ation and of gravity can be ignored. Special relativity can be cast into an abstract 
geometrical form based on the use of a flat 4-manifold called Minkowski space. In- 
deed, most of this can easily be handled using high school algebra, however this is 
not necessarily the most useful way of doing things — especially when you then want 
to generalize to the curved spacetime of the general relativity. Most of the special 
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relativity, especially the foundations, can be adequately treated using matrix algebra 

m 

The theory of special relativity asserts that spacetime is the manifold M 4 with a 
flat metric of Lorentz signature defined on it. In special relativity, the principle of 
general covariance states that the spacetime metric, r] a b, is the only quantity relating 
to spacetime structure can appear in any physical laws. In addition, the laws of physics 
in special relativity satisfy the principle of special covariance with respect to the proper 
Poincare transformation. 

In special relativity, energy is described in terms of the "time component" of 4- 
vector p a . The continuous matter distributions in special relativity are described by a 
symmetric tensor called the stress energy momentum tensor. 

In this chapter we introduce the theory of special relativity, which deals with the 
motion of objects as seen from inertial frames. General relativity extends this to cover 
non-inertial frames as well. As the theory of general relativity relies on the theory of 
special relativity, we present the concept of special relativity first, so general relativity 
is easier to understand. 
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Chapter 2 

Coordinate conditions in general 
relativity 

2.1 Introduction 



General relativity is a beautiful physical theory, invented for describing the gravita- 
tional field and the equations it obeys. General relativity is phrased in the language 
of differential geometry. It is by construction coordinate independent. 

General relativity does not have a preferred coordinate system, and the coordinates 
must be constructed along with the spacetime. The mathematics and physics cannot 
depend on the particular choice of coordinate system you choose to place on the 
manifold. However, the choice of the coordinate system can be useful in other ways. 
It may simplify the mathematical calculation and the physics interpretation of what 
is going on. 

For example, if the physics is time independent and spherically symmetric, it is 
extremely useful to choose a coordinate system that respects this property. 

Reference [TH] describes several coordinate conditions that we can consider to solve 
the vacuum field equations 
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for g a b. At first sight, the problem seems well posed. Indeed, there are 10 equations for 
the 10 unknowns g a t,. However, the equations are not independent, they are connected 
by 4 differential constraints through the contracted Bianchi identities 

V b G ah = 0. (2.2) 

We therefore seem to have a problem of under-determinacy, when there are fewer 
equations than unknowns. We cannot expect complete determinancy for any set g a t, 
since they can be transformed with fourfold freedom by a coordinate transformation: 

X a _ X' a = X' a (x) (a = 0,1, 2, 3). (2.3) 

In fact, we can use this coordinate freedom to impose four conditions on the g a b- 
These are known as coordinate conditions or gauge conditions. We can give an example 
by introducing Gaussian normal coordinates in which 

0oo->+l, Qoi^O (2 = 1,2,3). (2.4) 

The remaining six unknowns g a p can be determined by six independent equations 
in (1). 

The assumption of spherical symmetry of a spacetime means that each event be- 
longs to a two dimensional submanifold with spherical symmetry. This two dimensional 
submanifold can be given the metric 

ds 2 = R 2 (d8 2 + sin 2 d0 2 ) (2.5) 

where 9 and <fi are the standard angular coordinates on the sphere. There must exist 
(at least locally) two additional coordinates t and r labelling the different spheres. 

They can be chosen to be orthogonal to the sphere and also mutually orthogonal, 
t to be time-like and r space-like. The metric thus reads 

g ik = diag (-c 2 e 2< ^, e 2A , R 2 , R 2 sin 2 0) (2.6) 

where A = A(r, t), <p = <f>(r,t) and R = R(r,t). If ^ 0, it is possible to rescale R 
to be identical with r. However, this special choice of the coordinate system need not 
always be useful in particular cases. 
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In general relativity, we express physics in terms of tensors. General relativity 
is a geometrical theory. Specifically, we will deal with a Riemannian manifold which 
defines a metric. We used the metric to define a distance and length (norm) of vectors. 
We have infinitesimal distance or interval which is defined as 

(ds) 2 = ds 2 = g ab dx a dx b (2.7) 

where, ds in (|2.7|) is called a "line element", g ab is called a metric form. In addi- 
tion, g^b is called "covariant" metric and it has an inverse of the form g ab which is a 
"contravariant" metric. 

2.2 Independent components: Riemann's argument 

The quick way to see how much we can get away with is to simply count the number 
of independent components of the metric tensor. In n dimensions, since the metric is a 
symmetric T(0,2) tensor, it has n(n+l)/2 independent components. Since there are n 
coordinates, we can eliminate up to n of the metric components by suitably choosing 
our coordinates. This leaves 

n(n+l)/2 - n = n(n-l)/2 

true "independent" functions in the metric. This argument is due to Riemann. 

The number of free components of the Riemann tensor is strictly limited by all 
symmetries when we deal with a [torsion-free] metric connexion. 

2.2.1 For n=l 

If n — 1, then this implies there are zero truly independent components in the metric 
(see the proof below): In any 1-dimensional Riemann manifold we can always choose 
coordinates so that the metric is flat. (A space-time is said to be flat if R l jki = at 
every point in it.) 
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Proof. 



ds 2 = f(x) dx 2 = g xx (x) dx 2 . 



Let 



y = 



I 



Vf( x ) dx i 



then 




which implies 



dy 



Vfdx. 



So finally 



dy 2 = fdx 2 = ds 2 . 



□ 



i.e. we always can choose a coordinate y such that ds 2 = dy 2 in any 1-dimension 
manifold 



If n = 2, then this Riemann's argument implies there is exactly one truly independent 
component in the metric: In any 2-dimensional Riemann manifold it is a standard 
result that locally we can always choose coordinates to make the metric conformally 
flat. 

Indeed, when n = 2 the Riemann tensor has only one independent component, 
essentially the Ricci Scalar, R. We have (in any coordinate system) 



2.2.2 For n = 2 



2 R (9ac 9bd - g a d gbc) 

■^R 9ab 



Theorem: Any 2-dimensional Riemann manifold is locally conformally flat 



gij = exp [20(a)] 5, 
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2.2.3 For n = 3 



If n = 3, then Riemann's argument implies there are exactly three truly independent 
components in the metric: In any 3-dimensional Riemann manifold it is a standard 
result (typically attributed to Darboux) that locally we can always choose coordinates 
to make the metric diagonal: 



9ij 



9ii 
g 22 
#33 



That is, Riemann 3-manifolds have metric that are always locally diagonalizable. If 
n — 3, then the Riemann tensor has only six independent components, essentially R a t,. 
Indeed 

Rabcd = —2{ga[dRc]b + 9b[cRd]a} — R 9a[c9d]b (2-8) 

That is 



Rabcd — \9ac Rbd + Qbd Rac — 9 ad Rbc ~ Qbc R a d} 
1 

Lfi'ac 9bd ~ 9ad gbc} 



(2.9) 



2.2.4 For n = 4 

If n = 4 Riemann's argument implies there are exactly six truly independent compo- 
nents in the metric. This means that for a generic manifold in n = 4 you cannot assert 
that 





' 9n 
















922 








9ij = 








#33 
















#44 . 



because this metric has too few free functions (four rather than six). So in general 
we cannot in 4 dimensions choose patches to make the metric patchwise diagonal. In 
local coordinate patches, there are often non-zero off-diagonal elements. For a truly 
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general metric there must be at least 2 off diagonal elements. For instance, it might 
be true that we can always set: 



9ij 



#11 9i2 

#12 922 











#33 #34 

#34 #44 



This at least has the correct counting properties, and it would be nice to know if this 
sort of thing is always possible. Another thing that can always be done locally in 4 
dimensions is to set: 

"10 

922 923 924 
#23 #33 #34 
#24 #34 #44 

This is called a Gaussian normal coordinate patch. 

If n = 4 the Riemann tensor has only twenty independent components. Ten of 
them are the Ricci tensor R ao and the other ten are hidden in the "Weyl tensor" . 



9ij 



Cabcd — Rabcd + {9a[dRc\b + 9b[cRd]a} + -R9a[c9d]b 



(2.10) 



That is 

Rabcd ' 



Cabcd + ~{9acRbd + 9bdRac ~ 9adRbc ~ 9bcRad] ~ ~R{_9ac9bd ~ 9ad9bc} (2.11) 

2 o 



2.3 Lorentzian manifolds 

The Lorentz metric is used to formulate the free-space Maxwell equations in a 4- 
dimensional pseudo Riemann space and also used in describing mechanical phenomena 
in a mass-free Riemann space. The line element and metric are then assumed to be 

ds 2 = _ (dx ) 2 + (dx l ) 2 + (dx 2 ) 2 + (dx 3 ) 2 

= -c 2 dt 2 + (dx 2 + dy 2 + dz 2 ) (2.12) 
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This particular metric tensor will be called the Lorentz metric tensor, and the corre- 
sponding metric, the Lorentz metric. 



dab 











1 
1 
1 



The signature of a metric tensor is the number of positive and negative eigenvalues of 
the metric. That is, the corresponding real symmetric matrix is diagonalized, and the 
diagonal entries of each sign counted. If the matrix is n x n, the possible number of 
positive signs may take any value p from to n. Let q = n — p, then the signature may 
be denoted either by a pair of integers such as (p, q), or as the bracketed sum (p + q), 
or as a explicit list such as — h ++, or even as the number s = \p — q\. (Unfortunately 
all of these different notations are in common use). A Lorentzian metric is one with 
signature (p, 1) = (n — 1, 1) or sometimes (l,q) = (l,n— 1). For a Lorentzian signature: 

(0+1) dimension 

For a Riemannian geometry in (0+1) dimensions, this is trivial to get Lorentzian 
geometry 

gu = -exp[20(t)] 



(1+1) dimensions 

For a Riemannian geometry in (2+0) dimensions we can always locally choose coordi- 
nates such that 

g tj = exp[20(t,x)] 5 tj 
For Lorentzian geometry in (1+1) dimensions this becomes 



gij = exp[20(t,x)] 77^ 



where 



Vij 



-1 
+1 
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so 



9ij 



gij = exp[20(t, x)\ 

— exp[20(t, x)] 




-1 
+1 



exp[20(£, x)] 



We can find the Ricci scalar n in term of <fr, its derivative, rj, the Ricci tensor and 
Riemann tensor. The Ricci scalar is 



R 



dx' 



(/)(t,x] 



df 



<f>(t,x] 



exp(20(t, x)) 

And we can calculate the Ricci tensor, its non-zero components are: 



(2.13) 



Rn 
R22 



dx 2 
d 2 



(p(t,x] 



dx' 



4>(t, x 



dT 2 ^ x) 



d 2 

dT 2 ^ x - 



(2.14) 
(2.15) 



Similarly, we can calculate the Riemann tensor, its non-zero components are: 

R 1212 = - exp(20(t, x)) (- <f>(t, x)) + 0(t, x)) ) (2.16) 



In Lorentzian geometry, it is possible to make another coordinate choice by setting 

-[1-v 2 ] -v 



9ij 



—v 



where v = v(t,x), so 

" -[1 - v(t,x) 2 ] -v(t,x) 
-v(t,x) +1 

In terms of this coordinate system we have: The Ricci scalar is 



9ij 



R = -2 ( ^ v(t, x)) - 2 (l v(t, x)) - 2v(t, x) (* v(t, x) 



(2.17) 



20 



And we can calculate the Ricci tensor, its non-zero components are: 



R u = -{-l + v{t,x)) 



(2.18) 



Rl2 = v(t,x) 



2 



dxdt 







dx 



v(t, x)) + [ ir v(t, x) + v(t, x) — v(t, x) } (2.19) 



2 



dx 2 



R 



22 



d 2 \ (d \ 2 ( d 2 

— V {t, x )\-1- V {t,x) \ - v{t, x) I — 2 v{t, X] 



[2.20) 



Similary, we can calculate the Riemann tensor, its non-zero components are: 



R 



1212 



d 2 



dxdt 



v(t, x) 



d y ( d 2 

m v(t,z)\ -v(t,x)(—v(t,x) 



(2.21) 



(2+1) dimensions 

The simplest special coordinate patch that it is always possible to set up yields 



9ij 



-h 2 
h x 2 
h 2 2 



where hi = hi(t,x,y), so 



9ij 



-h (t,x,y) 2 

h^x^y) 2 

h(t,x,y) 2 



Alternatively we can always locally set up patches of the form 



9ij 



— [N 2 — v 2 ] vi v 2 

V! 1 

v 2 1 



21 



where V{ = Vi(t, x,y),N — N(t, x, y), v 2 = v\ + v\. 
So the metric is 

-[N(t,x,y) 2 - (vf(t,x,y) +v%(t,x,y))] Vt(t,x,y) v 2 (t,x,y) 
vi(t,x,y) 1 

v 2 (t,x,y) 1 

Coordinate patches of this type are called Painleve-Gullstrand coordinates. Another 
alternative, in (2+1) dimensions, is that we can always locally set up patches of the 
form 

-[1-v 2 ] V\ v 2 
9ij = Ui exp[20] 

v 2 exp[ 

where = </>(t, x, y),v 2 = exp[— 20] (vf + v 2 ). 
The inverse metric is then: 

— 1 vi exp [—20] v 2 exp [—20] 

g 13 — t>iexp[— 20] (exp[20] — vf) exp[— 40] —ViV 2 exp[— 40] 

f 2 exp[— 20] —ViV 2 exp [—40] (exp[20] — v|) exp[— 40] 

In (2+1) dimensions we can calculate the Ricci scalar, the Ricci tensor and the Rie- 
mann tensor. However, all results are too complex to analyze in detail. 



(3+1) dimensions 

This is the physically interesting case, and this is the case that will be most difficult 
to analyze in detail, and the remainder of the Chapter is devoted to this case. 



2.4 The block diagonal ansatz 

• It is certainly true that we can always locally set: 



9u 


012 








012 


022 














033 


034 








034 


044 
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• This at least has the correct counting properties. 

Now we can take the above as an ansatz and calculate the Ricci tensor, and then the 
Einstein tensor, respectively. 

Specifically, We write the four coordinates as = (y a , z l ), with 



x 



x 



y 



X 



X 



z . 



Then 



9vA x ) = g^(y,z) 



9u(y,z) 


9u(y,z) 










922(y,z) 














93s(y,z) 


934 (y, z) 








934(V,Z) 


gu(y,z) 



That is 



9nu(y,z) 



h a b(y, z) 








fa(y, z) 



This means the 4 dimensional metric has split up into four 2x2 blocks, two of them 
trivial (zero). This will make the Ricci tensor relatively easy to calculate. 
The inverse tensor is very simple 



/<!/ 



h ab (y,z) 








f j (y,z) 



9 

where h ab and / y are two dimensional inverses of the two 2x2 matrices h a b and /y. 

For definiteness, we assume the h subspace has signature (1+1) [one time direc- 
tion plus one space direction] while the / subspace has signature (2+0) [two space 
directions]. But this subtle point does not really affect the subsequent analysis. 

The Christoffel symbols can be computed from 

1 



9 



9v{a,p) 



:9ap,u 



'2.22) 



and now we can spilt them up into little 2x2 blocks 



ap 



(y,z) 



' ^ n (y,z) 
^12 (y, z) 


^22(y,z) 


^is(y,z) 
^2s(y,z) 


^u(y,z) ' 
T^(y,z) 


1^13 (y, z) 
^ u (y,z) 


^2 3 (y,z) 
^2*(y,z) 


^3 3 (y,z) 


T^(y,z) 
r>*u(y,z) 
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Some of these blocks should have a simple interpretation in terms of the Christoffel 
symbols of the 2 dimensional metrics h ab and fy. 
For example 

T a bc (y, z) = g au |^(6, C ) - -fl'&c.z/J = 9 ad ^9d(b,c) - ^9bc,d^ 

= h ad \h d{b , c) - l -h bc ^ = [T(h)] a bc (y, z). (2.23) 

That is, this portion of the 4 dimensional Christoffel symbol can be calculated simply 
by looking at the 2 dimensional Christoffel symbol corresponding to the two dimen- 
sional metric h ab . Similarly, we can calculate this: 

F l jk(y, z) = g w {g v m - ^9jk,uj = g d {gim - \gjk,i 

= W l {him ~ \hjk,^ = in/WAV, z). (2.24) 
The complications of course come with the other pieces, such as 

r\ c (y, z) = g lv {g u (b,c) - 7}9bc,v^ = P jfl^M - 7}9bc,j 

= r{o-\h bc ^ = ~ph bc , r (2.25) 
Similarly, and we can calculate this: 

F a jk(y, z) = g av j#i/(j,fc) - -j9jh,v | = h ab {gbm ~ 2 gjh ' b 

= h ab jo - 1/ A6 J = -\h ab f 3Kb . (2.26) 

This part of the 4-dimensional Christoffel symbol is rather simple but does not have 
a simple interpretation in terms of the 2 dimensional sub-geometries. 
There are four other sub-blocks to calculate: 

F a bk(y, z) = g av (g„(b,k) ~ ^9bk,u^j = g ac (g c {b,k) ~ ^9bk,^j (2.27) 

= g ac (^-g c b,k + -^gck,b — -gbk,cj = h ac (-h^k^ = -h ac h cb ± 
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Whence 



F a jc (y,z) = -h°»h bc>j 



(2.28) 



Similarly 



T l JC {y, z) = g iv [ g v{hc) - ^g jc A = g ik \g k {j,c) - \gjc,hj (2.29) 

ik I | fik ( f | cik r 

9 \ lj9kj,c + ~y9kc,j ~ 7j9jc,k J — J I ~^Jkj,c ) = 2 3jk,c 



Whence 



T\ k (y,z)= l -p f jk>b 



(2.30) 



Various symmetries should make all four of these results obvious once any one of 
them has been calculated. 

2.4.1 Collected results for the connexion 

There are a total of 8 sub-blocks to deal with: 





= [T{h)]\ c {y,z 


F l jk(y,z) 


= nmAviz 


r l bc(y, z) 




r a jk(y,z) 


= ~2~ h<lb fjkfi- 


T a bk (y,z) 


= -h ac h cbk . 


r a JC (y,z) 


= hhr 7- 

2 J 


r l jc(2/, z) 


fik v 
— 2^ Jjk,c- 


r\k(y,z) 


= l^f 3 fjk,b- 



2.4.2 Back to the Ricci tensor: qualitative results 

Once we have all eight sub-blocks, we have the entire 4 dimensional Christoffel symbol. 
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The Ricci tensor of can then be calculated from ^7], page 224 



R 



pv 



rfj _ per 

pi;, a J- pa,u 



per -pp _ -pp per 
pa ^ pv up L pv 



(2.31) 



It is convenient to first calculate T a 



pa 



In fact, it is a standard result that 



T % ff = (In y/=g) 



2 9 



(2.32) 



(see, for instance, equation (8.51a) on page 222 of fTTj). Then we can simplify: 

1 



R 



pv 



-g 



(^ r V), ff -(inv^),„-r% r* 



pvi 



(2.33) 



(see, for instance, equation (8.51b) on page 222 of [T7J). This is the explicit formula 
for the Ricci tensor. 

To interpret the Ricci tensor, write it in 2 x 2 sub-blocks 



Rpu(y,z) 



Rn(y,z) 

Rn{y,z) 


Rn(y,z) 
R22(y,z) 


Ri3(y,z) 
R23(y,z) 


Ru(y,z) 
R2i(y,z) 


Ri3(y,z) 
Ru(y,z) 


R23(y,z) 
R2i{y,z) 


Rs3(y,z) 
Ru(y,z) 


R^{y,z) 
R 4 4(y, z) 



That is 



RpAy,z) 



R a b(y, z) 


Rai(y,z) 


_ Rjb(y, z) 


Rij{y,z) 



Some of these blocks should have a simple interpretation in terms of the Ricci 
tensor of the 2 dimensional metrics h a b and For instance we expect 

[R(g)) ab = [R(h)U+?T! 

[R(g)h = [«(/)]«+??? 

where the extra terms (the ???) involve y derivatives of / and z derivatives of h. That 
is, we expect the full 4 dimensional Ricci tensor of these particular 2x2 sub-blocks 
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to be related to the Ricci tensors of the 2 dimensional metrics / and h, but with some 
additional terms. 

For the off-diagonal block 

[R(g)]ia =??? 

It will probably involve y derivatives of / and z derivatives of h. 
Once you have the Ricci tensor, the Ricci scalar is easy 

R(g) = g^R^ = h ab R ab + = R(h)+m + £(/)+??? 

That is 

R(g)=R(h) + R{f)+m 

That is, we expect the full 4 dimensional Ricci scalar to be related to the sum of the 
Ricci tensors of the 2 dimensional metrics / and h, but with some additional terms. 
Then for the four dimensional Einstein tensor 

Gyv = Rfiu — ^ ^ 9nv (2.34) 
we can sub-divide the 2x2 sub blocks as 

[G(g)] ab = [G(h)} ab - ~ R(f) h ab +m 

[G(g)h = [G(f)h - \ R(g) fa+m 
But the Einstein tensor of any two dimensional manifold is zero; this is one of the 
special features of two dimensions, so 

[G(h)] ab = = [G{f)} l3 

and we have 

[G(g)} ab = - 1 - R(f) h ab +m 

[G(g)h = ~\ Rig) fa+m 

Finally, for the off-diagonal 

[G(g)] aj = [R(g)] ai =??? 

Of course it is all these ??? pieces that give the important mathematics and physics. 
Again, we are not really sure what to expect, they will probably involve y derivatives 
of / and z derivatives of h. 
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2.4.3 Calculating the Ricci tensor - diagonal pieces 



Start from the general expression for the Ricci tensor of the metric g^: 



Note that, since the metric is block diagonal, 

g = de%^] = det[h ab ] det [/;.,■] = h f. 

Now consider 



(2.35) 



(2.36) 



1 



-9 



[R(g)U - 

which we expand as 

mg)U = ■4(/ ; i7""J+H 7 (/5Ji'-),-M, 



(2.37) 



V=h7 v v J ""Ac ^hj 

pc pd pc p/c pj pA: pj pd 



ab 



Then 



ab 



a6 



,k 



(.n^S) -(lnV7 



rc pd pc pfc pj pfc pj pd 

da ^ cb J- fca J- c6 J- fca J- j'6 J- da J- j'fe 



which we regroup as 



i 



ab 



l-h 
lie 

2 / 



-hY c 



ab 



\nV-h 



.ab 



■pc -pd 
1 da 1 cb 



T C ab + 



hfT\ b ) - W/ 



rc pfc pj pfc pj pd 

fca J- eft J- fca J- j'ft J- da J- jfe 



,aft 



(2.38) 



But the first line is just [R(h)] ab , the 2-dimensional Ricci tensor calculated using the 
2-metric h ab . The other terms can be slightly rearranged as 



[R(g)U = [R{h)]ab-\{{\nf), ab - f ^T\ b } + -^= (y^h]v k ab ) 



rC "pfc r\j -p/c t-ij pid 

/ca -L cb J- fea -L j& -L da J- j"6 



(2.39) 
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Now use colons : to denote 2-dimensional covariant derivatives using the 2-metric h. 
Then 

[R( g )U = [R(h)U- 1 -(\nf). ab + 1 -^r k ab + ^= (v^r 



n \ v J ao i 

fca J- cb J- ka J- j'b ^ da J- j'fe 



-r c fc a r fc cfc - r J fca r fe Jb - r J da r ,- 6 (2.40) 



Now we start to insert the various pieces. 
Note that 

We can write this more compactly as 

;j?(^ r *-), = -5 A '*- (2 ' 41) 

where by A j we mean the two-dimensional Laplacian in the / subspace, treating the 
h indices as though they were scalars. 
Also note that 



which we can re-write as 



r^fca r jb .f^ fmk,a f fnj,b (2.42) 



-pj -pk _ rij £ rkl f 

L ka L jb J jk,a J J li,b 

= itr [r 1 f a r 1 f 6 ] (2.43) 



Next, we start looking at 



r C fca r fc co — ~~^J lCd hda,k) (/^ hbc,j] 



Similarly, 
Finally 



\h cd (h ad , P h hc4 ) (2.44) 



T j da T d ]b = -h cd (h ad , P h bc>j ) (2.45) 



17 1 ^ 
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Now pull all these pieces together. We get 

[R(g)] ab = [R(h)] ab - ~ (ln/) :o6 - i A, h ab - ~tr [f" 1 f a f 1 f J 

fX&j + V (h ad ,i f h bcJ ) (2.47) 
There are still simplifications we could make. For instance 

(I"/):,,, (Jf) y-yy (2.48) 
Using this we get rid of the logarithm. Then finally 

[R(g)] ab = [R(h)] ab - \^ + 1^ ^ - \ A f h ab - ^tr [f 1 f a f" 1 f b ] 

~Th fijhab ' j + \ hCd {had ' 1 fJ hbc ' j) (2,49) 

For [R(g)]ij we can appeal to symmetry to simply write: 

[R(g)h = [R(f)h -l^ + ^T^-l^^-t 1 ^ hi h " h ^ 

~\j h ab h, a + \f kl (f ik ,a h ab f jljb ) (2.50) 

The vertical bar | now denotes 2-dimensional covariant derivatives in the / direction, 
treating the h indices as scalars. 
To obtain the Ricci scalar we note 

1 A f If f :a 1 1 
h ab [R(g)] ab = h ab [R(h)] ab ---JLL + -lj J —- -h ab A f h ab - -h ab tr [f" 1 f a f 1 f b ] 

lll + ^' tr [ h " h!h " h ^ (2 ' 51) 
with a similar formula holding for h l i [R^g)]^: 

~rf y + \ hah tr I rl f - a f-1 ^ (2 - 52) 
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Then in the Ricci scalar there are a lot of partial cancellations 

Rig) = h ab [R(g)] ab + P [R(g)] l3 

y v 1AJ lA f ft If - a f :a lhu h\ { 

-V Af h ab - ^f ij A h fij 

+^h ab tr [r 1 f a f" 1 f 6 ] + l/« tr [h- 1 hi h- 1 hj] (2.53) 
2.4.4 Off-diagonal part of the Ricci tensor 

Now we consider what we can say about the off-diagonal piece [R(g)] a i- Again we will 
want to start from the general formula for the Ricci tensor of the metric g^ v \ 

*V = ^= W^g rV), CT - (in v^,, - r% v\ v . 

Now consider 

[Rig)]* = -j= (V^9 T\ t ) a - (In ^) a . - YP aa T° pi (2.54) 
which we expand as 

rc pd -pc pfc pj pfc pj pd 

da ci J- fca J- ci J- fca J- j'i J- da J- 

That is 

~\ M-h)) M - \ (ln(/)),„, 

r ^ r c j r \ a r " C j r ^ R'da r jj (2.55) 

So that 

[#(#)]<« = -^(M/)), a , + ^(ln(/)) )6 r 6 ai + - 7 L= (V^hT b ai ) b 
~\ (H~h)), ai + ~ (H-h)),j P ai + -±= (v7 1%,) . 

r da r c j r r C j R'fca r jj r^da r (2.56) 
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Now consider 



p6 



1 



ai:b 



On 



-hT" 



ab 



(ry 



(2.57) 



Here T b a i is treated as though it were a tensor with respect to the h subspace, the 
% index [which lies in the / subspace] is treated as though this quantity were a scalar 
with respect to coordinate transformations in the / subspace. This notation, together 
with the analogous 

- rj (r**«) 



J- 7. 



ia\j 



77* 



(2.5J 



allows us to simplify the off-diagonal Ricci tensor as: 

[R(g)U ~ 



^ai: b + r j ia \ j -l(f, a /f)^-l(h i /h) a 

+\ (f,b/f) r b ai + l - (hj/h) r ai - l -h a ^h cd f^f 3hd 



To simplify things, we can define two tensors by 

K a 
U 



■~h n 
2 ° 



2^3,0- 



ly.a 



(2.59) 

(2.60) 
(2.61) 



with the convention that the abc... indices are raised and lowered using h, and the 
ijk... indices are raised and lowered using /. 

(The reason for this notation is that these quantities are some sort of generalization 
of the notion of extrinsic curvature.) 
With these conventions we have: 

T a bc (y,z 

T a bk (y,z 
T a jc (y,z 



F l jc (y, z 
F l bk(y,z 



[T(h)} a bc(y,z). 

K bc \ 

T a 

—K a bk- 
—K a 



f.) ■ 

- 77 ■ 
-L % kb- 
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Which then implies 



—K ai -, b + K bi:a ~ 




(2.62) 



2.4.5 Discussion 

This form of the metric is physically interesting because there is a good chance that it 
really is general for 4 dimensional spacetime, and general enough to encode the Kerr 
metric for a rotating black hole. 

For this chapter we simply counted the number of independent components of the 
metric tensor. In 1-dimension, there are zero truly independent components in the 
metric. Furthermore, in any 1-dimensional Riemann manifold we are able to choose 
the coordinates to make the metric flat. For 2-dimensions, Riemann's argument im- 
plies that there is only one truly independent component in the metric. Indeed, any 
2-dimensional Riemann manifold has a standard result that locally, we can choose 
coordinates to make the metric conformally flat. The Riemann tensor also has one 
independent component, essentially the Ricci scalar, R. In addition, for n — 3, then 
Riemann's argument implies that there are exactly three truly independent compo- 
nents in the metric. The Riemann 3-manifolds have metrics that are always locally 
diagonalizable. For n = 4, the Riemann tensor has only twenty independent compo- 
nents. Ten of them are the Ricci tensor R a b, while the other ten are hidden in the 
"Weyl tensor" . 

We use the Lorentz metric to formulate the free-space Maxwell equations in a 
4-dimensional Riemann space, and to describe mechanical phenomena in a vacuum 
Riemann space. It is trivial to get a Lorentzian geometry from a Riemannian geometry 
in (0+1) dimensions. For a Riemannian geometry in (1+1) dimensions we are able to 
find the Ricci scalar, the Ricci tensor and the Riemann tensor. While, for a Riemannian 
geometry in (2+1) dimensions it is more complicated to calculate the Ricci scalar, Ricci 
tensor and Riemann tensor. 

The case for (3+1) dimensions is physically interesting, but it is the most difficult 
to analyze in detail. However, we are able to calculate the Ricci tensor and Einstein 
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tensor in the 4 dimensional metric. We split the metric into four 2x2 blocks, two of 
them trivial (zero), which allows us to calculate the Ricci tensor easily. Furthermore, 
the Christoffel symbols can be computed from equation (|2.22|) . For the 4-dimensional 
Christoffel symbol, there is a total of 8 sub-blocks to deal with (see section (j2.4.1|l ). 
Now we are able to calculate the Ricci tensor of g^ v from equation (J2.31)) . We expect 
the full 4 dimensional Ricci tensor to be related to the sum of the Ricci tensors of 
the 2 dimensional metrics / and h, with some additional terms. In addition, we can 
calculate the 4 dimensional Einstein tensor from equation (|2.34j) . This is useful as a 
way of looking for nice new solutions of the Einstein equations. 

2.5 The ADM (Arnowitt-Deser-Misner) ansatz 

Both ADM and reverse- ADM (Kaluza-Klein) decompositions seek to split "spacetime" 
into "space" plus "time" in some sort of natural way, and then use this split to describe 
spacetime curvature in terms of the curvature of space and some additional structure. 

In addition, this decomposes the four-metric of spacetime as a time- varying three- 
metric on a three-dimensional domain together with the lapse and shift, which are 
scalar-valued and three- vector-valued functions of time and space, respectively. 

2.5.1 ADM Decomposition 

Definition 1. In the case of the ADM decomposition one takes 







' -[c 2 -v 2 ] 




-Vi 





Notation: Latin indices from the beginning of the alphabet (a,b,c,...) run from 0-3 
and refer to space-time; Latin indices from the middle of the alphabet (i,j,k...) run 
from 1-3 and refer to space. We use X- a to denote a space-time covariant derivative 
and X :i to denote a three-space covariant derivative. 

Interpretation: The light cones, defined by ds 2 = g a b dx a dx b = take the form 

hij {dx l - v l dt) (dx j - v j dt) = c 2 dt 2 
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Here hu is the 3-metric of space. We take 



to be the inverse 3-metric, and 



h 13 Vi Vj = hij v l v J 



Then v i is the 3-velocity of a coordinate system with respect to which c is the "coor- 
dinate speed of light" . (Note that in comparison to usual conventions v i is minus the 
"shift vector", while c is equal to the "lapse function".) 
We are able to calculate the inverse metric: 



ab 



' -1/c 2 


v 3 /c 2 


—v l /c 2 


h t] — v l v : '/c 2 



(J 



Note that 3 contravariant vectors that point in the "space" directions are 

X a = (0; 1,0,0), Y a = (0; 0,1,0), and Z a = (0; 0, 0, 1). 
So if we want n to be a vector that is "orthogonal" to space, then we need to pick its 
covariant components to satisfy 

n a cx (1; 0,0,0) 

If we want n to be a unit vector that is "orthogonal" to space, then we need to pick 
its covariant components to satisfy 



n a = ± (c; 0,0,0) 
Its contravariant components are then 

\c c J c 

Choosing n a to be a future directed unit vector orthogonal to the space slices then 
fixes 

'l t v*\ _ (lju*) 
, c c 

and 

n a = (-c; 0,0,0). 



n 
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2.5.2 Reverse- ADM (Kaluza-Klein) Decomposition 

Definition 2. The "reverse-ADM decomposition" can be found in Landau- Lifshitz and 
consists of setting 



' -c 2 










hij - Vi Vj/c 2 



Notation: Again, Latin indices from the begining of the alphabet (such as a,b,c,...) 
run from 0-3 and refer to space-time; Latin indices from the middle of the alphabet 
(i, j, k, ...) run from 1-3 and refer to space. We use X. a to denote a space-time covariant 
derivative and X ;i to denote a three-space covariant derivative. 

Interpretation: The light cones, defined by ds 2 = g a b dx a dx b = now take the form 

dx l dx j = c 2 (dt + Vi dx 1 ) 2 

We again say that h^ is the 3-metric of space, though it's a different 3-metric than in 
the ADM case. We take 

h ij = [h.j]- 1 

to be the inverse 3-metric, and 

v 2 = h 1 -* ViVj = h^ 

Then c is the "two-way coordinate-averaged speed of light" , which is by construction 
isotropic. See §84 on pages 251-255 of [T3| for details. The physical interpretation of 
the 3- velocity v % is not particularly clear, and there seems to be no clean interpretation 
of the current formalism in terms of "analogue models" ■ Nevertheless it is one of 
the standard metric decompositions, so let us adopt it and see what happens. 
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2.6 Static spacetimes: ADM and reverse- ADM (Kaluza- 
Klein) Decomposition 

Our strategy is straightforward: 

• Take any static spacetime, and use the natural time coordinate to slice it into 
space plus time. 

• Use the Gauss-Codazzi and Gauss- Weingarten equations to decompose the 
(3+l)-dimensional spacetime curvature tensor in terms of the 3-dimensional spatial 
curvature tensor, the extrinsic curvature of the time slice [zero], and the gravitational 
potential. 

In any static spacetime one can decompose the spacetime metric into block diagonal 
form [1, 2, 3]: 

ds 2 = dx^dx u 

= -c 2 dt 2 + dx l dx j . (2.63) 

Here c and hij are assumed to be time-independent, they depend only on position in 
space. This is equivalent to taking 



9ab = 



-c 2 








h^ 



Theorem: In any static spacetime coordinates can be found in which the ADM and 
reverse- ADM (Kaluza-Klein) decompositions are identical. The proof is by inspection. 
Simply set v = in the general formalism. 

Being static tightly constrains the space-time geometry in terms of the three-geometry 
of space on a constant time slice, and the manner in which this three-geometry is 
embedded into the spacetime. For example, from jTTj, page 515, and applying page 
518, we have the standard results (which we will re-derive in a slightly more general 
context below): 



(3+1) A--*., 


^ ^ Rijkl- 


(2.64) 


^' R-tijk 


= 0. 


(2.65) 


Rtitj 




(2.66) 
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This will generalize when we allow time dependence for c and hij. 
Now taking suitable contractions, 



(3+1) p.. 

{3+1) Ru 



(3) p.. _ 

C 

0. 



( 3 ) Ac. 



So that 



( 3+1 ) J R= ®i?-2 



W Ac 



(2.67) 

(2.68) 
(2.69) 

(2.70) 



Finally, for the spacetime Einstein tensor (for bits of this see [T7|, page 552, noting 
that in our present case = 0, while in the corresponding formulae of ^7] they chose 
c= 1) 



(3+1) G^ 
(3+1) G t , 
{3+1) G tt 



•3 (• _ 

c 

0. 



+7 l3)R - 



9ij 



W Ac 



(2.71) 
(2.72) 
(2.73) 



This decomposition is generic to any static spacetime. (You can check this decom- 
position against various standard textbooks to make sure the coefficients are correct. 
For instance see [Tl IT2l I2T] ) . 



2.7 Zero— shift spacetimes: ADM and reverse- ADM 
(Kaluza-Klein) Decomposition 

If the spacetime is not static but the "shift vector" is zero we still have 



9ab 



' -c 2 











(2.74) 



(Such a decomposition is always possible in a globally hyperbolic spacetime.) 

But now we will have to permit time dependence for both c and h^. Note that as 
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long as the shift vector is zero the ADM and reversed-ADM decompositions remain 
identical. 

In this case it is useful to define a quantity called the extrinsic curvature [Hj, the 
equation (21.67) 

= 1 dh^ 

3 2c dt v ' 

and then 



{3+1) R zjkl = WRiw + KnKji-KuKfr. (2.76) 
{3+1) R tljk = c{K lk:] -K ir . k }. (2.77) 
^R titj = cc :ij + "mess". (2.78) 



See [T7j, equations (21.75) and (21.76). Note (n ■ n) = —1, and the "n-coordinate" of 
equation (21.76) is really "iVt" [that is, "ct"]. The last equation for ^ +1 ^R tit j is not 
given by [IZj, except implicitly on page 518, and this really is a bit of a mess which 
we shall explicitly calculate below. 

But some pieces of the Riemann, Ricci, and Einstein can be immediately read off 
from ^7j. For instance from [TJj, equation (21.8) 



R tl = G ti = N{K lj - j -triK)..,} 



while from [T7], equation (21.77) 



G tt = -N 2 { {3) R + (trK) 2 - tr(K 2 )} 
2 
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2.7.1 Systematic derivation 

Christoffel symbols of the first kind: 

1 1 



r«t = g t (t,t) - ^ 9tt,t = ^ gtt,t = {c 2 ), t = -cc >t (2.79) 



Fat = gi(t,t) - 2 9tt,i = -\gtt,t= +^ (c 2 ),t = +cc t (2.80) 

11 1 

r«i = T m = g t{ t,i) - - gu,t = ^ 9tt,i = ( c2 )/i = ~ cc ,* ( 2 - 81 ) 

r*ij = gt(i,j) - - gij,t = gij,t = -- h ijjt = cKij (2.82) 

111 

r iy = ^ijt = gi(t,j) - - gtj,i = - gij,t = - h ijtt = -cKij (2.83) 



— o S'i'M ~~ ht(j,k) — 7, hjk,i — ^Tyfc (2.84) 



Raising the first index: 

Christoffel symbols of the second kind: 



r*« = C± (2-85) 

c 

r ft = +ch ij C)j = +cc j (2.86) 

r** = r* 4i = ^ (2.87) 
c 

r*y = -^f (2.88) 

r* ti = r jr = -ciCj (2.89) 



{3) r jk (2.90) 



2.7.2 Riemann tensor 

Now insert this into the formula for the Riemann tensor. 

R bcd = - 2 {F a b{c,d\ - r a e[c r e | 6 | d] } (2.9i) 

where the \b\ indicates we do not anti-symmetrize over the b. Lowering the a index 

Rabcd = —2 {r ab [ Cjd ] + g a e\dF e \b\c] ~ ^ \b\d] } (2.92) 
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and rearranging 

Rabcd = —2 {r ab [ Cjd ] + T ea [ c g ef V\fb\d\ } (2.93) 

This is equivalent, for instance, to the unnumbered equation on page 21 of ^UJ , which 
can further be simplified to 

Rabcd = - {9ad,bc ~ 9bd,ac ~ 9ac,bd + 9bc,ad} + ^ ead 9^ T/fcc _ ^ eac 9^ ^ fbd (2.94) 

This is equation (11.6) of ^H]. So far this is completely general. If g a b block diagonal- 
izes in the manner g a b = — c 2 @ hij then we have 

Rabcd 77 \9ad,bc 9bd,ac 9ac,bd 9bc,ad\ — 2 ^tad ^tbc 

H — n Tiac F t bd + h mn T mat i T n b c — h mn T mac T n bd (2.95) 

Rijki : If we now look at the space components of the 4-dimensional Riemann tensor 
we have 

Rijkl = ~ {giLjk — 9jl,ik — 9ik,jl + 9jk,il} J ^ til ^ t i k 

— c 

+ ^2~ Ttifc r tj -( + h mn T mi i T n jk — h mn T mi k T n ji (2.96) 
But most of this is simply the 3-dimensional Riemann tensor 

Rijkl = ^ Rijkl — ^til^tjk + Ttifc F t ji (2.97) 

and the rest we recognize in terms of the Kij we defined earlier, so that 

Rijkl = ^ Rijkl ~ K tt K jk + K ik Kji (2.98) 

which is the explicit way of doing things. 

R-tijk : in a similar manner we can evaluate Rujk as 

Rtijk — 7^ {dtk,ij ~ 9ik,tj ~ 9tj,ik + 9ij,tk} ~ ~^ ^ttk ^tij 

+ ^2" r«j r*jfc + h mn T mt k T ni j — h rnn T mt j Tnik (2.99) 
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which we can re-write as 



Rtijk — ( c Kik),j ~ \pKij\k ~ c ,jKik + C^Kij — cK U ] z T n ij + cK U j Tnifc 

so that 

Rtijk cKikj cKij k C K ^nij C K j ^nik 

that is, in terms of 3-space covariant derivatives 

Rt ijk = c{K ik:j -K ij:k } (2.100) 
This is the second part of what we were trying to prove. 

Rtitj : It only remains to calculate Rutj- (After all Rtuj, Run, Rtttt and so on are all 
zero by symmetries of the Riemann tensor.) 
So let us calculate 

1 1 

Rtitj — — {dtj,it — 9ij,tt — 9tt,ij + 9it,tj} 2 ^ ut 

— c 

+ ~ r*ti F t ij + h rnn T mt j T n u — h mn T mt t T n ij 



So that 



Rutj — - {-hij,tt + (c 2 ),*?} 



1 1 

(- cc j) (-cc,i) + — (-cc t ) (cKij) 

+ h mn (- cKmJ ) (- C K ni ) - h mn (cC m ) T nij 

That is: 

Rtitj (^-^V;),t ~t~ CC^j ~\~ C^i Cj Cj C{ C^Kij ~\~ C (Ki m h K n j} (cC ;m ) r ij 

So finally 

Rtuj = cK iht + <■<■,. j + {K 2 )ij (2.101) 
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Summary: Whenever the metric g a b block diagonalizes in the manner g a t, = —c 2 hy, 
then even if we allow arbitrary time dependence we have the relativity simple results 

RijH = ^Rijki — Ku K jk + K ik Kji (2.102) 
Rtijk = c {K ik:j - K ij:k } (2.103) 
R titj = cK ijit + cc :i:j + (K 2 ) ij (2.104) 



where 



1 dhij 



LiJ 2c dt 

Of course once we have deduced these results for the Riemann tensor, results for the 
Ricci tensor and Einstein tensor can be obtained by simple contraction. 



2.8 General stationary spacetime: ADM Decom- 
position 

For a general stationary spacetime we need to use the full ADM decomposition 





-[c 2 


-v 2 } 


~ V 3 








Vi 








-1/c 2 


-vJ/c 2 






-v l /c 2 


h ij - 





We do however have the simplification that all quantities are time-independent. (All 
time derivatives vanish). Note that the ADM decomposition is in this case quite 
distinct from the reverse-ADM decomposition. 
We still have, thanks to Dirac 

Rabcd = - {9ad,bc ~ 9bd,ac ~ 9ac,bd + 9bc,ad] + ^ ead 9^ T/fcc _ ^ eac 9^ ^ fbd (2.105) 
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But now when we insert the inverse metric 



,bc 9bd,ac 9ac,bd 9bc,ad} 

9 ^tad ^tbc H o r^oc Ttfod 
C C 

— (f m /c ){r tad r m6c + r ma(J ri &c } + (t> m /c ){r mac r iM 

_ <,"v/ c 2) r mad r nfec - (/r™ - ^"/c 2 ) r mac r n6d (2.106) 

Rijki : Start by looking at the space components Rijki- We have 
1 

Rijki — -\9il,jk — 9jl,ik — gik,lj + 9jk,il\ 

{V I C ){^til ^mjk ~l~ ^ mil I" \jk tik ^mjl J 

-(v m v n /c 2 ) Y ma T njk - (v m v n /c 2 ) Y mik r njl 

+ h mn T m u T n jk — /l™ T mik T n ji (2.107) 

But most of this is simply the 3- dimensional Riemann tensor, and the rest can easily 
be regrouped as 

Rijki = ^ Rijki 2^ til Vm ^mil) (r*jfc + V n T n j k ) H ^{T tik + V m T mik ) (T t j[ + v n T n ji) 

(2.108) 

Now suppose we define 

K l3 = l -{T tlJ + v m T miJ } (2.109) 

This is a natural generalization of our previous definition, and would allow us to again 
write 

Rijki = {3) Rijki - K u K jk + K ik K jX (2. 110) 

but now for this "improved" definition of K^. Let's check what happens. Having de- 
layed things as long as possible, we'll have to start by calculating Christoffel symbols. 
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2.8.1 Christoffel symbols of first kind 

(Remember, stationary, no t dependence): 

11 . 
tttt = 9t(t,t) ~ j 9tt,t = ^ 9tt,t = (2-lH) 

r itt = s-i(t, t ) - - #t t ,t = -- g- tt ,i = +^ (c 2 - w 2 ),; (2.112) 

r«i = r« t = ^ t (t,o - 2 = o ^ = _ 2~ ^ 2 ~ (2.ii3) 

r«j = - - 0y,t = +^ (gu,j + 0tj,i) (2.114) 

1 / 

Titj = r^t = ^(tj) - - g t j t i = - (g it ,j - gtj,i) 

= -r - v j,i) = - v [i,j] (2.115) 
^ijk = <7i(j,fc) ~~ 2 ^' fc '* = ^Cfi*) ~~ 2 = (2.116) 

Then we have 

l -{T tij + v m T mij } = l - {-v {iJ) + Vm^ij} = — v m (2.117) 
so that in this situation 

K tj = — v m (2.118) 

This is compatible with the general definition of extrinstic curvature, see ^7j equation 
(21.67), note that our v l is minus the shift vector. 

Rtijk: In a similar manner we can evaluate Rujk as 

Rtijk = 2 {dtk,ij ~ gik,tj — gtj,ik + gij,tk} 
J 1 1 

— [v m I C 2 ) {Tftk Tmij + Tmtfc ^tij} + (f" 1 / (?) {^mtj ^ \ik + ^ttj ^mik} 

+ ( h ™ _ v ™ v n /c 2j Tmtk ^ _ {h mn _ ^nj^ p^ . p^ (2.119) 
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which we can re-write as 
1 

Rtijk — ~^\~ v k,ij + v j,ik\ 



(r ttk + v m v mtk ) (r tij + v m r mij ) + - (r ttj + v m r mtj ) (r« fc + v m r mik 

h ^mtk^nij h ^mtj^nik (2.120) 



Now 



r« fc + v m T mtk = -^(c 2 - v 2 ) tk - ^v m (v m>k - v k , m ) = -^(c 2 - v 2 )-. k - ^v m (v m:k - v k:m ) 

(2.121) 

so that 

Tttk + v m T mtk = -^(c 2 ) jfe + l -v m {v m± + v k:m ) = -c {c k + v m K rnk ] (2.122) 
So now 

Rtijk = „{ — u *w + — (c^ k + v m K rnk )Kij + (cj + v m k m j)K ik — V[ m>k ]T m ij +V[ m j]T m ik 

(2.123) 

This can be rewritten in several ways: If we remain in a coordinate basis 

Rtijk = - v [m,k] rm ji ~ v [jM rm ki - (c, fe + v m K mk )K i:j + (cj + v m K mj )K ik (2. 124) 
whence 

= - (c fe + v m K mk )K ij + ( Cj - + v m K mj )K ik (2.125) 

This is relatively compact, but not so easy to compare to our previous result. 
But we can also write 

2^ [j:k]:i ^j:k:i ^k:j:i Vj:i:k ^ ^ R jki Vk'.i'.j ^ R kji 

Vj:i:k ^i-.j-.k ^k:i:j ^i:k:j ^Vrt ( R jki ^ R kji ^ ^ R ijk) 
= 2 (u(iy) : fc — l>(i:fc) : j — f m ^ R™ ijk) 

= 2([-cK ij \, k -[-cK ik \, j -v m ^R m ijk ) (2.126) 
which lets us write 

R tljk = c {[K lk ] :k - [K^]*} - u m pi^jfc - K mfc i^ + K mj i^] (2.127) 
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that is, in view of our previous result, 

Rtijk = c[K ik .j — Kij. k ] — v m R mi j k (2.128) 

This is the tidiest we can make things in a coordinate basis. This is the analogue of 
PU, equation (21.76). 

For future reference define 

% = V [i,j] = V [i-j] (2.129) 

and note the identity 

n jk:i = [cK ik } :J - [cKijlj, - v m ®R m ijk (2.130) 

Note that it sometimes easier to work in what is referred to as a "non-coordinate basis" 
by defining: 

Rnijk = n a Raijk (2.131) 

in which case 

Rnijk — — (1; ^ m ) {Rtijk] Rmijk) = ~ [Rtijk + ^Rmijk] = Kik:j ~ K-ij-.h- (2.132) 

This means you are not really looking in the "t" direction, but are instead looking in 
the "n" direction normal to "space" . 

Rtiy: Now for the Rutj pieces. This is where things might get a little messy. From 
Dirac's general formula, applied to the present situation 

Rtitj — — {gtj,it ~~ 9ij,tt ~~ 9tt,ij + 9it,tj} 

— {v m / c 2 ) {T tt j T mit + T mt j T tit } + (v m / c 2 ) {T mtt T ti j + Y ttt T mi j} 

+ ( h ™ _ v m v n /c 2^ Fmtj _ {h mn _ A »^| ^ ^ . ^33) 



Because of time independence this greatly simplifies 
1 I l 1 

- \-gtt,ij} - - 

.m I „2 



Rutj — - {-gtt,ij} - {r«j F U t} 



— (v m I c ) {T tt j T mit + T mt j Tut ~ ^mtt r%} 

+ ( h ™ _ v m v n /c 2^ mt . _ (r _ v m v n /c 2 )TmU ^ . (3 ^ 
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so that 



Rutj — 2 {~ 9tt,ij} ~ O^ttj + v m T mt j) (Tut + v m T m it) 

+ {v m / c 2 )T mtt T ti j + h mn T mt j T nit — h mn T mtt T ni j (2.135) 

Inserting the known values of the Christoffel symbols (calculated above) 
1 1 

+/i m > m ,^ M - l - h mn (c 2 - v 2 ), m T nij (2.136) 

That is 

Rtitj = \(c 2 - v 2 ).. l3 - ( CJ + v m K mj ){c^ + v m K mi ) + \ (v m /c) (c 2 - v 2 ), m Kj 

+h mn v [m>j] v [n>i] (2.137) 

But then 

Rutj c c : jj — (f ) ; jj t' [cj K m { -\- C j -ft. m j] f K m j v K n i 

+v m c m Kn - \ (v m /c) (v 2 ) >m ICy + h mn v [mJ] v [nA (2.138) 

Which still needs to be simplified further. 
For instance 

^y 2 U = \(v k v k)-.ij = {v k u fc :i)y = [v k ..j v^) + {v k v k±j ) (2.139) 

But 

v k v k . A = h kl v k . 4 v k . A = h kl (- cK kj + v^) (- cK ' u + v [l:i] ) 

= c 2 (K 2 )ij — c(KQ — QK)ij + h mn V[ m j] V[ nA (2.140) 

where with obvious notation Qij = V[ij] and 

{KQ - nK)ij = K ik g kl % - n ik g kl K Xj (2. 141) 
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Aside: Write (KQ-Q K)^ as 

Kg^Vt - n g- x K 

then 

(Kg- 1 Sl-ng- 1 K) T = Q T [g'Y K T - K T [g~ 1 } T Q T 

= -n g- 1 K + K g' 1 (2.142) 

so this object is symmetric, as long as the minus sign is included. 
Then 

Rutj = cc. Aj -(v k v k±j )-c 2 (K 2 ) ij + c(Ktt-ttK) ij 

+v m c )m Kij - (v m /c) (v k v k:m ) K l3 (2.143) 

so that 

Rutj = c c-ij - (v k v k ,i,j) - c 2 (K 2 ) i:j + c(KVt - £lK)ij 

+v m c m Kij + (v k v l K kl ) (2.144) 

Now by the Ricci identities the combination v k v k -_i-j is symmetric in ij, so that we can 
write: 

v k v k:i:j = v k v k .(i tj ) (2.145) 

But now consider 

Vk:(i:j) = 2 {''/.:-:/ + 1 = 7> {{~cK ki + Q ki ).j + (~cK kj + Q kj ) :i } 

= --{(cK ki ) :j + (cK kj ). A } + -{^ kilj + n kj . A } (2.146) 
Then recall the identity 

n jk . A = [cK lk } : , - [cK^ - v m {3) R m ljk (2.147) 
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to deduce 

v k:{tj) = [cKijU + l -v m ( {3) R m ljk + (3) iT» (2.148) 

whence 

v k v Hi:j) = v^cKy]* - v k v l ^R kllj (2.149) 
Putting this back into Rutj we have 

Rutj = c c.ij - v^cKijU + v k v l (3) R k uj - c 2 {K%j + c(KQ - QK)^ 

-v m K mj v n K ni + v m c m K l3 + (v k v l K kl ) K iS (2. 150) 

which simplifies to 



Rutj = c c :i j - v k cK ij:k - c 2 (K 2 )ij + c(KQ - QK)^ 

-v m K mj v n K ni + (v k v l K kl ) Kij + v k v l ®R kilj (2.151) 



We can rewrite the last line in terms of the 4-curvature as 

Rutj = c c :i j — v k cKij. k — c 2 (K 2 )ij + c(KVt — VtK)^ 

+v k v l R kilj (2.152) 

This is the simplest expression we can get in a coordinate basis. 

Rninj : Now define in a "non-coordinate basis" . 

Rninj = n a n b Raibj = i\\V % Y (l',V l ) b R a ibj 

= ^2 [Rutj + v k R uk j + v k R kitj + v k v l R ki ij] 

= [Rutj - cv k [K ik ,j - K ij± ] - cv k [K jk ,i - K ij:k \ - v k v l R kilj ] 

[R titj - cv k [K iktj + K jk:i - 2K ij:k ] - v k v l R kilj ] (2. 153) 



c 2 



Then 



c 2 Rninj = c c-ij - v k c K ij:k - c 2 (K 2 )ij + c(KQ - QK)ij 

-c v k [K ik:j + K jk . A - 2K ij:k ] (2.154) 
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There are now some simplifications 

c 2 Rninj = cc-ij + cv k K ij:k - c 2 (K 2 )ij + c(KQ - QK)ij 

-cv k [K ik:j + K jk:i ] (2.155) 

To simplify this we consider the Lie derivative 

CvKij = v k d k K i3 + v\i K kj + v k j K ik (2.156) 
which can also be written in terms of 3-covariant derivatives as 

CvKij = v k K ijtk + v k . A K kj + v k . tj K ik (2.157) 

so that 

C v Kij = v k K ij:k + h k \-cK kl + n ki )K tj + (-cK kj + Q kj )K u ] (2.158) 

and 

CvKij = v k K ij:k - 2c {K%j + (KQ - QK)^ (2.159) 

Therefore 

C 2 Rninj = C C.ij + cC v Kij + c 2 {K 2 )ij - c v k [K ik:j + K jk:i ] (2.160) 

Summary : 

For an arbitrary stationary spacetime. taking 

K a = -- v (i--j) 

we have 

Rijki = {3) R i j k i-K il K jk + K ik Kj l (2.161) 

Rtijk = c [K ik ,j — Kij. k ] — v m Rmijk (2.162) 
Rmj = cc :ij -v k cK ij:k -c 2 {K 2 ) ij + c{KQ-QK) lj + v k v l R kilj (2.163) 
We can also write this as 

Rnijk = n a Raijk = Kij, k — K ik ,j. (2.164) 

and 

Rninj = n a n b R aibj = c c tij + cC v Kij + c 2 {K 2 )ij - c v k [K ik .j + K jk:i \. (2.165) 
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2.8.2 Completely general ADM decomposition 

For the general ADM case coordinate techniques are too messy to contemplate. The 
easy results are on page 514 of [Hj and boil down to 



K, 



1,1 



2<% 9ij + 



(2.166) 



note that Kij is now a combination of the zero-shift and stationary case. 
Then 

Rijkl = Rijkl — Kil Kjk + Kik Kji 

and 

Rnijk = ^ Raijk ^ij:k ^ik:j 

The difficult one is then encoded in page 518, ^7] which seems to assert that 



(2.167) 
(2.168) 



Rnink = n a n b R aibj = [C n K]ij + [K 2 ]^ + a (i;i) + a { aj (2.169) 
where Oj is the space projection of the 4-acceleration a = (n • V) n. 

2.8.3 The inverse- ADM decomposition 



with 



follows 




' -1/iV 2 


v t /N 2 


_ v r /N 2 


hij - ViVj/N 2 


= Vi(t,x,y,z), and % = 



(2.170) 



The inverse metric is now 



-[N 2 - WviVj 



h 13 



(2.171) 



where v l = h 1 ^ Vj and h %: > = [hij] 1 
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This is the reverse of ADM formulation 



and 



9u 


912 


913 


9u 




-1/N 2 




v 2 /N 2 




V3/N 2 




v A /N 2 


921 


922 


923 


924 




v 2 /N 2 


h 2 2 


- v 2 v 2 /N 2 


h 2 3 


- v 2 v 3 /N 2 


h 2i 


- v 2 v 4 /N 2 


931 


932 


933 


934 




V3/N 2 


h 3 2 


- v 3 v 2 /N 2 


h 3 3 


- v 3 v 3 /N 2 


h 3 A 


- v 3 v 4 /N 2 


9u 


912 


943 


944 




v 4 /N 2 


h>42 


- v 4 v 2 /N 2 


fl43 


- v iV3 /N 2 


I144 


- v A v A /N 2 



-[N'-Et^ih^v,} 


v 2 


V 3 


v 4 ' 


v 2 


h 22 


h 23 


h 24 




h 32 


h 33 


h 34 


v 4 


h 42 


h 43 


h 44 



For the time being, let Vi = then 



9p.v 



9u 


912 


913 


914 




-1/N 2 











921 


922 


923 


924 







h 22 


h 23 


h 2 4 


931 


932 


933 


934 







h 3 2 


h 3 3 


hu 


941 


942 


943 


944 







h i2 


h 4 3 


h<44 



and 



9 U 


9 12 


9 13 


9 U 




-N 2 











9 21 


g 22 


g 23 


9 24 







h 22 


h 23 


h 24 


9 31 


g 32 


g 33 


9 U 







h 32 


h 33 


h 34 


9 41 


9 A2 


9 43 


9 Ai 







h 42 


h 43 


h 44 



(2.172) 



(2.173) 



(2.174) 



(2.175) 



The 4 dimensional metric can be spilt up into one 3x3 block, and two 1x3 blocks 
with two of them trivial (zero). That is 



' -l/N 2 








hij 



The inverse tensor is very simple 



(2.176) 



' -AT2 








h ij 



(2.177) 
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where — N 2 and h l i are two dimensional inverses of the lxl matrix and 3x3 matrices 
— 1/N 2 and hij, respectively. 

The Christoffel symbols can be computed from 

r% = gT{g^, p) - \g* PlV } (2.178) 

and you can spilt them up into little sub-blocks 





r" 12 


r" 13 


1^14 " 


r^i 2 




r" 23 




r^i 3 


r^ 23 


r^ 33 


r^ 3 4 




1^24 


r^ 34 


r^ 44 



(2.179) 
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2.8.4 Christoffel symbols of the 2 metrics — 1/N 2 and h 



ij 

Some of these blocks should have a simple interpretation in terms of Christoffel sym- 
bols of the lower dimensional metrics —1/N 2 and hij. 

When t = time and i, j = space, for example 

r*« = g tv {gu(t,t) - ^9u,u} = g tt {gt(t,t) - -xgttA 

ttf 1 1 1 i 

— g i^gtt.t + -gu,t — 2 

= g tt {\gu,} = -\n 2 {(-^- 2 )A = (2.18O) 
= g tv {gu{i,j) - -^g^A = g tt {gt(i,j) - ^9vA 

= g t \0- 1 -h lJ A = N 2 { 1 -h. lJ ,}= 1 -N 2 {h^} (2.181) 

1 1 

r u = g tu {gu{t,i) - -^guA = g u {gt(t,i) - -^gt%A 

ttf 1 1 1 1 

— g \-^gtt,i + -^gti.t — -^gti.t) 

= 9 tt {\ guA = \ iV 2 {(-^),} = ^ (2-182) 

. iv 1 .. 1 

rv = g iv {gu{t,t) - -gtt,u} = g l3 {gj{t,t) - -gtt,j} 

= gHO - I huA = ~\ /''•'{( = ~h ij ^ (2.183) 

1 1 

r 4 jfc = g w {g v {j,k) - - gjk,A = g l {gi(j,k) - ^ 9jk ' 1 ^ 

= h a {h im - l - h m } = [Tih^jk (2.184) 

r V = g lu {gv(j,t) ~ 2 9jt ' u ^ = glk { gk{ jt) ~ 2 9jt ' k ^ 

lk A 1 1 1 

— g i^gkjt + 2 9kt,j - 2 9jt,ki 

= \g lk g kh t = \h lk h kht (2.185) 
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2.8.5 Collected results for the connection 



There are a total of 6 sub-blocks to deal with: 




diN 



(2.187) 



(2.186) 



(2.188) 



(2.189) 
(2.190) 
(2.191) 



2.9 Conclusion 

In this section, we analysed a generic spacetime using both ADM and reverse- ADM 
(Kaluza-Klein) decompositions which split "spacetime" into "space" plus "time" in a 
natural way. We use the split to describe spacetime curvature in terms of the curvature 
of space and some additional structure. Furthermore, this decomposition encodes the 
four-metric of spacetime in terms of a time- varying three-metric on a three-dimensional 
domain together with the lapse and shift. 

In addition, we used the Gauss-Codazzi and Gauss-Weingarten equations to decom- 
pose the (3+l)-dimensional spacetime curvature tensor in terms of the 3-dimensional 
spatial curvature tensor, the extrinsic curvature of the time slice [zero] , and the grav- 
itational potential. The ADM and reverse- ADM decompositions are identical in any 
static spacetime. In addition, we derived the standard results in equation ()2.64|) to 
(J2.66j) for the spacetime geometry in terms of the three geometry of space on a con- 
stant time slice, in which this three geometry is embedded into the spacetime. Finally, 
we derived equation (J2.71j) to (J2.73J) for the spacetime Einstein tensor. 

We consider time dependent but block diagonal metrics in the ADM and reverse- 
ADM decompositions. When the "shift vector" is zero, the ADM and reverse-ADM 
decompositions are identical. Futhermore, we defined a quantity called the extrinsic 
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curvature in equation (|2.75j) . We are able to calculate the Christoffel symbols of 
the first kind and second kind from section (|2.7.1|) . Furthermore, we calculated the 
Riemann tensor by using the formula from (J2.91|) to ()2.93|) . Whenever the metric g a b 
block diagonalizes in the manner g a b = — c 2 @ hij, then even if we allow arbitrary time 
dependence we have the relatively simple results as equation ()2.1()2|) to (j2.1()4j) . 

For a general stationary spacetime we need to use the full ADM decomposition. 
However, we have the simplification that all quantities are time independent. Indeed, 
the case for ADM decomposition is quite distinct from the reverse-ADM decomposi- 
tion. The Christoffel symbols first kind are given in section (|2.8.1j) . We derived the 
simplest expression for the Riemann tensor in terms of a coordinate basis. 

For the general ADM case coordinate techinques are too complicated to contem- 
plate (page 514 of ^Zj). Using the inverse- ADM decomposition we can calculate the 
Christoffel symbols from equation (j2.178|) . Furthermore, we can split them up into sub- 
blocks in equation (j2.179j) . Some of these blocks can interpreted in terms of Christoffel 
symbols of the 2 metrics —1/N 2 and see ([2.8.4)) . 
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Chapter 3 



Perfect fluid spheres in general 
relativity 

3.1 Introduction 

Perfect fluid spheres are interesting because they are first approximations to any at- 
tempt at building a realistic model for a general relativistic star |9"j ITT ] IT5 | I2()j. 

Although they represent a real and important astrophysical situation, explicit and 
full general solutions of the perfect fluid constraints in general relativity have only 
very recently been developed. 

The first static spherically symmetric perfect fluid solution with constant density 
was found by Schwarzschild in 1918. 

He found two exact solutions to Einstein's field equation — the "exterior solution" 
(relevant outside the star) and the "interior solution" (an approximation to what 
goes an inside the star). The so-called Schwarzschild solution is amazingly simple, 
especially in view of the field equations which are very complicated. Indeed, the 
exterior Schwarzschild solution is not a degenerate case for over-simplified situations 
but physically most meaningful. It is this solution by means of which one can explain 
most general relativistic effects in the planetary system. The reason is that it describes 
the gravitational field outside of a spherically symmetric body like the planets and the 
sun. 
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Over the past 90 years a confusing tangle of specific perfect fluid spheres has been 
discovered, with most of these examples seemingly independent from each other. 

Many "ad hoc" tricks have been explored to solve this differential equation, often 
by picking special coordinate systems, or making simple ansatze for one or other of 
the metric components. The big change over the last several years has been the 
introduction of "algorithmic" techniques that permit one to generate large classess of 
perfect fluid spheres in a purely mechanical way. Perfect fluid spheres may be "simple" , 
but they still provide a surprisingly rich mathematical and physical structure. 

In this chapter we will extend these algorithmic ideas, by proving several solution- 
generating theorems of varying levels of complexity. We shall then explore the formal 
properties of these solution-generating theorems and then will use these theorems to 
classify some of the previously known exact solution, and additionally will generate 
several previously unknown perfect fluid solutions. 

3.2 Six theorems for generating perfect fluid spheres 

Perfect fluid spheres, either Newtonian or relativistic, are the first approximations 
in developing realistic stellar models. Because of the importance of these models, 
fully general solutions of the perfect fluid constraint in general relativity have been 
developed over the last several years. 

The central idea is to start solely with spherical symmetry, which implies that in 
orthonormal components the stress energy tensor takes the form 



p 














Pr 














Pt 














Pt . 



(3.1) 



and then use the perfect fluid constraint p r = p t . This is making the radial pressure 
equal to the transverse pressure. By using the Einstein equations, plus spherical 
symmetry, the equality p r = pt for the pressures becomes the statement 

G ee = G fr = G$$- (3-2) 
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In terms of the metric components, this leads to an ordinary differential equation 
[ODE], which then constrains the spacetime geometry, for any perfect fluid sphere. 

Over the last 90 years, many "ad hoc" approaches to solving this differential equa- 
tion have been explored. We often solve it by picking special coordinate systems, or 
making simple ansatze for one or other of the metric components |S1 IH El 13 IHl ESj • (For 
recent overviews see [OIEIED]- For a table of "well-known" perfect fluid solutions see 
Table on page O) The big recent change has been the discovery of "algorithmic" 
techniques that permit one to generate large classes of perfect fluid spheres in a purely 
mechanical way fSJEHEI]- in addition, an alternative algorithmic approach is used 
to generate any spherically symmetric static solution of Einstein's equation without 
recourse to the details of the material that makes it up. 

Solution generating theorems 
Algorithmic techniques 

We would like to present "Algorithmic" techniques ^H] that permit one to generate 
large classes of perfect fluid spheres in a purely mechanical way. The condition of 
hydrostatic eqilibrium for Newtonian gravity is 



where g(r) is the local acceleration due to gravity, m(r) is the mass contained within 
a radius r and p{r) is the density. 

We can simply calculate the mass within a radius r in terms of the density: 



Jo 

When we used p(r) as a known function, we have to solve a second-order differential 
equation: 
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We can solve for P r (r) in terms of an integral of the density: 



P r (r)= / J — — p(r)dr + C 2 . 



We let g(r) be an arbitrary function to solve simple algebra which allows us to find 
m(r). The differentiation allows us to find p(r) and then P r (r), respectively: 



In both the Newtonian and relatistic cases we have to ensure that the density is non- 
negative. The pressure is finite and non-negative. 

Spherical symmetry 

Spherically symmetric means "having the same symmetry as a sphere". Indeed, 
"sphere" means S 2 , not spheres of higher dimension. Rather than the metric on a 
differentiable manifold, we are concerned with those metrics that have such symme- 
tries. We can characterize symmetries of the metric. By the existence of Killing 
vectors, which we can see from the definition below. By demanding spherical sym- 
metry we can simplify the equations of motion considerably. It turns out that for 
vaccuum spherical symmetry of the metric automatically means the metric is static 
and asymptotically flat. This is Birkhoff's theorem, a proof of which can be found in 
Furthermore, Spherical symmetry can be defined rigorously in terms of Killing 
vector fields as follows. 

Definition 3. A space-time is said to be spherically symmetric if and only if it admits 
three linearly independent spacelike Killing vector fields X a whose orbits are closed 
(i.e. topological circles) and which satisfy 




=> p{r) 



2rg(r) + r 



2dg 



4irr 2 




[X\X 2 ] =x 3 



[X 2 ,X 3 ] = X 



1 



[X 3 ^ 1 ] = X 2 
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Then there exists a coordinate system in which the Killing vectors take on a stan- 
dard form as expressed in the following result. 

Proposition 1. In a spherically symmetric space-time, there exists a coordinate sys- 
tem (x a ) (called quasi Cartesian) in which the Killing fields X a are of the form 



X a — > uj a pX , LO aj 3 — —UJp a . 

In addition, we know that what the Killing vectors of S 2 are, and that there are 
three of them. Therefore, a spherically symmetric manifold is one that has three 
Killing vector fields which are the commutator of the Killing vectors which is the same 
in either case of those on S 2 . 

Schwarzschild solution 

The gravitational field of a homogenous spherically symmetric body is derived in 
firstly Newton's and then in Einstein's gravitational theory. The gravitational field of 
a simple Einsteinian model star consists of the exterior and the interior Schwarzschild 
solutions. They are joined together at the surface of the star. 

The use of arbitrary coordinates is permitted in general relativity. Indeed, the 
physical significance of statements about tensor or vector components are other quan- 
tities which are not always obvious. However, there are some situations where the 
interpretation is almost as straightforward as in special relativity. The center point of 
a local inertial coordinate system is the most obvious example. Its coordinate system, 
where the principle of equivalence comes from, allows one to treat all local quantities 
exactly as in special relativity. Furthermore, Schwarzschild coordinates for a spherical 
system turn out to be a second example [T7] . 

Schwarzschild geometry 

The geometry of a spherical symmetric vacuum, i.e. vacuum spacetime outside the 
spherical black hole is the "Schwarzschild geometry" can be described in terms of the 
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Schwarzschild metric, 

^-(l-^tf + j-^ + Atf, (3.3) 

which was derived originally as the external field of a static star. The Schwarzschild 
metric seems to have a singularity at the surface where r = 2M due to its coordinates, 
in which space and time change their meanings. 

Schwarzschild geometry illustrates clearly the highly non-Euclidean character of 
spacetime geometry when gravity becomes strong. Furthermore, it illustrates many 
of techniques one can use to analyze strong gravitational fields. When appropriately 
truncated, it is the spacetime geometry of a black hole and of a collapsing star as well 
as of a wormhole [T7j . 

However, that the Schwarzschild exterior solution is the unique spherically sym- 
metric vacuum solution, is known as Birkhoff 's theorem. This is a theorem of general 
relativity which states that all spherical gravitaional fields, whether from a star or 
from a black hole, are indistinguishable at a large distances. A consequence of this is 
that purely radial changes in a spherical star do not affect its external gravitational 
field. 

That the Schwarzschild geometry is relevant to gravitational collapse follows from 
"Birkhoff's theorem". 

Birkhoff 's theorem: Let the geometry of a given region of spacetime be spherically 
symmetric, and be a solution to the Einstein field equations in vacuum. Then that 
geometry is necessarily a piece of the Schwarzschild geometry |17j . 

In particular, Birkhoff 's theorem implies that if a spherically symmetric source like 
a star changes its size, however does so always remaining spherically symmetric, then 
it cannot propagate any disturbances into the surrounding space. From figure 13.11 
shows that a pulsating spherically symmetric star cannot emit gravitational waves. If 
a spherically symmetric source is restricted to the region r ^ a for some a, then the so- 
lution for r > a must be the Schwarzschild solution, or can be called the Schwarzschild 
exterior solution. However, the converse is not true: a source which gives rise to an 
exterior Schwarzschild solution is not necessarily spherically symmetric. Some counter- 
examples are known. Therefore, in general, a source need not inherit the symmetry of 
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its external field. If we take the limit of equation 13.31 as r — > oo, then the flat space 
Minkowski metric of special relativity in spherical polar coordinates is obtained, which 
is 

ds 2 = -dt 2 + dr 2 + r 2 (d6 2 + sin 2 6 d<f) 2 ). (3.4) 

This shows that a spherically symmetric vacuum solution is necessarily asymptotically 
flat Eg. 




Figure 3.1: A Pulsating spherical star cannot emit gravitational waves. 

Firstly, start with some static spherically symmetric geometry in Schwarzschild 
(curvature) coordinates 

ds 2 = -C(r) 2 dt 2 + + r 2 dn 2 (3.5) 

B(r) 

and then assume it represents a perfect fluid sphere. That is, = G f f = Git. While 
Gqq = Gtt is always fulfilled due to spherical symmetry, we get a nontrivial constraint 
from Gqq = Gff- 
We calculate 

2Br(' -( + (B 
Gff = t? (3.6) 
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and 

1 B'C + 2BC + 2Br(" + r('B' 
G§§ ~ ~2 < 

We set Gf-r = Gs§, which supplies us with an ODE: 

[r(r()']B' + [2r 2 C" - 2(rC)']£ + 2( = . (3.8) 

This reduces the freedom to choose the two functions in equation (J3.5j) to one. 

This equation is a first order-linear non-homogeneous equation in B(r). Thus, once 
you have chosen a ((r) — this equation can always be solved for B(r). Solving for 
B(r) in terms of £(r) is the basis of [T3*| lift], (and is the basis for Theorem 1 below). 

On the other hand, we can also re-group this same equation as 

2r 2 B(" + (r 2 B' - 2rB)(' + (rB' - 2B + 2)C = , (3.9) 

which is a linear homogeneous 2nd order ODE for C( r ); which will become the basis 
for Theorem 2 below. Thus, once you have chosen a B(r) — this equation can always 
be solved for C( r )- The question we are going to answer in this section is, how can one 
systematically "deform" this geometry while still retaining the perfect fluid property? 
That is, suppose we start with the specific geometry defined by 

ds 2 = -Co(r) 2 dt 2 + + rW (3.10) 

and assume it represents a perfect fluid sphere. In the following section, we will show 
how to "deform" this solution by applying five different transformation theorems on 
{Co, B }, such that the outcome still presents a perfect fluid sphere. The outcome 
of this process will depend on one or more free parameters, and so automatically 
generates an entire family of perfect fluid spheres of which the original starting point 
is only one member. 

In addition, we analyze what happens if we apply these theorems more than once, 
iterating them in various ways. We also try to find the connection between all six 
different transformation theorems. 

66 



(3.7) 



3.2.1 The first theorem 

The first theorem we present is a variant of a result first explicitly published in [To] , 
though another variant of this result can also be found in j2Q] . We first re-phrase the 
theorem in ^Bj. This is slightly different formalism, and demonstrate an independent 
way of proving it. 

Definition (Idempotence): We call a transformation T idempotent if T 2 = T. 
Speaking loosely, we will call a theorem idempotent if whenever we apply the trans- 
formation more than once, no further solutions will be obtained. 

We use our proof to show that the transformation in theorem 1 is (in a certain sense 
to be made precise below) "idempotent". By applying theorem 1 more than once, no 
new solutions will be found. Furthermore, the idempotence property of theorem 1 and 
theorem 2, which we will soon establish, enables us to divide the class of perfect fluid 
spheres into seed metrics and non-seed metrics. 

Theorem 1. Suppose {Co( r ), B (r)} represents a perfect fluid sphere. Define 
Co(r) V , „„„ L fC^r) Co(r)-r<S(r) 



\(o{r)+r( Q (r)J { J (o{r)(o{r)+rQ{r) J 

Then for all \, the geometry defined by holding Co( r ) fixed and setting 

d, 2 = -Co(r) 2 dt 2 + btTXYXm + rW ( 3 ' 12 ) 
B (r) + A A (r) 

is also a perfect fluid sphere. That is, the mapping 

U\) : {Co, A)} - {(o,B + AA (Co)} (3.13) 

takes perfect fluid spheres into perfect fluid spheres. Furthermore a second application 
of the transformation does not yield new information, Ti = TjQ^ Aj) is "idempotent" , 
in the sense that 

HX n ) o ■ • • o T 1 (X 2 ) o T 1 (X 1 ) = Ai) : {Co, B } » {Co, B + (^ =i A,) A (Co)} 

(3.14) 
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We also note that Tj always has an inverse 

IUXT 1 = 7~i(-A). (3.15) 

Proof for Theorem 1. Assume that {Co( r )? B^ir)} is a solution for equation ()3.8|) . 

We want to find under what conditions does {Co( r ) ; -^i( r )} a l so satisfy equation 
(J3.8)) ? Without loss of generality, we write 

By{r) = £? (r) + A A (r) . (3.16) 

Equation ()3.8|) can now be used to determine A (r). Substitute B\{r) in equation 
(EH) 

[r(rCo)1(S + A A )' + [2r 2 ^ - 2(rCo) / ](i?o + A A ) + 2( = 
we can also re-group this same equation as 

[r«o)K + [2r 2 Co - 2«o)']£o + 2( + [r«o)'] A' + [2r 2 ^ - 2«o)'] A = 

That ordinary inhomogeneous first-order differential equation in B now simplifies to 

KrCo)'] A' + [2r 2 C - 2(rCo)'] A = , (3.17) 

which is an ordinary homogeneous first-order differential equation in A . 
Now we can calculate Aq by following, 



[r(rCo)1 A' = - [2r 2 ^ - 2«o)'] A . (3.18) 



Rearrange it into 



AJ, - [2r 2 C - 2(rCo)1 
A [r(rCo)'] 
-2rC | 2 
(rCo)' r 

A straightforward calculation, including an integration by parts, leads to 



(3.19) 



A (r) = — - — ? exp ( / |> . (3.20) 

[Ko)f U Ko)' 
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We can also write this same equation as 

,v,2 

exp { 2 / 



[«o)f 



7 Co (Co + rCo) 7 Co (Co + rCo) J 



Adding and subtracting ±2(rC ) / (Co( r Co)') to the argument under the integral leads 
to 

A ° w = ^ e -{ 2 /f(^*'- 2 /l(arb* 



[(rCo)f I 7 Co (Co + rft) J Co (Co + <5) 

We can simplify to, 

-^ exp { 2 /^(c. : <) dr + 2 /ci dr 



[K„)f I y Co (Co + <J) J Co 

Furthermore, 

^ 2 „„„ fo /"Cfi (Co-rCo) 



[«o)1 

Note that, exp(ln(Co)) = Co ' so we can calculate 



An _ ( CoW y f Ko(r) Co(r)-r Co (r) 1 

°" Uo(r)+rCo(r)J P \ 7 Co(0 Co(r) + r Co (r) ' ^ 

as previously illustrated. 

If we apply this transformation a second time we obtain no additional information. 
To see this, consider the sequence 

{Co, B } ^ {Co, B 1 } i ► {Co, B 2 }... (3.24) 
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But at the second step (and all subsequent steps), since Co has not changed, so Aj(r) = 
A (r). More generally, at all stages of the iteration Aj(r) = A (r). We can write this 

as 

n / n \ 

i=l \i=l J 

or in the more suggestive form 



Jt x = Ti (3.26) 



i=i 

where the symbol = indicates "equality up to relabelling of the parameters" . That 
is, transformation Ti is "idempotent" up to relabelling of the parameters (see figure 

E2J. □ 



A version of Theorem 1 can also be found in [201 - Specifically, after several manip- 
ulations, changes of notation, and a change of coordinate system, the transformation 
exhibited in equation (16.11) of can be converted into Theorem 1 above. 

Applying theorem 1 to a fixed {( ,B } generates a one dimensional space of perfect 
fluid spheres, which leads to the corollary below: 

Corollary 1. Let {(q, B a } and {£o,-B&} both represent perfect fluid spheres, then for 
all p 

{( ,pB a + (l-p)B b } (3.27) 

is also a perfect fluid sphere, furthermore all perfect fluid spheres for a fixed £o can be 
written in this form. 

Proof. The result is automatic once you note that for fixed (o the ODE for B is linear, 
(though not homogeneous, which is why the two coefficients p and 1 — p are chosen to 
add up to 1) . □ 

We defer extensive discussion of the application of this theorem and its corollary 
until section 13.41 where we use this and our other generating theorems as a basis 
for classifying perfect fluid spheres. At this stage we mention, only as a very simple 
example, that Ti applied to Minkowski space results in the geometry of the Einstein 
static universe. 



70 




01 23456789 10 



r 

Figure 3.2: The solid lines show B(r) for 10 reapplications of Theorem 1 onto the Minkowski 
metric. The dashed line corresponds to a single application with a specific choice for Aonce- 
It can be seen that 10 applications of Theorem 1 can be expressed by one application. 



3.2.2 The second theorem 

The second theorem we present is a new transformation theorem. This is a different 
formalism, and demonstrates our independent way of proving it. Furthermore, we 
would like to show the "idempotence" property of theorem 2 in this section. 



Theorem 2. Let {(o, Bo} describe a perfect fluid sphere. Define 

rr i \ I r dr 
Z {r) = a + e I 

J Co(r) 2 v^oW 
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seed metric Minkowski 



old solution: 
Einstein static 



Figure 3.3: This structure shows that when we apply theorem 1 to Minkowski, we get an 
Einstein static universe as a new metric. 

Then for all a and e, the geometry defined by holding B (r) fixed and setting 

ds 2 = -Co(r) 2 Z (rf dt 2 + — ^ + r W (3.29) 
is also a perfect fluid sphere. That is, the mapping 

T 2 (a, e) : {Co, B } ^ {Co Z (( 0l B ), B } (3.30) 

takes perfect fluid spheres into perfect fluid spheres. Furthermore a second application 
of the transformation does not yield new information, T 2 is "idempotent" in the sense 
that 

T 2 (cr n , e n ) o • • • o r 2 (cr 3 , e 3 ) o T 2 (cr 2 , e 2 ) o T 2 (<Ti, e x ) = T 2 (a n . . . a 3 a 2 cri, e n ... 32 i), (3.31) 
where 

£n.. .321 = - ■ -0- n ) + (0-^ 1 € 2 0- 3 - ■ + (o-{ 1 (T2 1 € 3 - ■ -0- n ) 

+ • • • + (<rr V 2 - V3- 1 •••e n ). (3.32) 
Furthermore, theorem 2 is invertible (as long as a ^ 0): 

[T^e)}- 1 ^ T 2 (l/a,-e). (3.33) 
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Reduction of order: 

This is a method to find a general solution to a linear differential equation, provided 
you already have one particular solution. In this method, we start with an n^-order 
linear differential equation. This method is especially useful in solving a second-order 
linear differential equation. It reduces the problem to one for solving a first-order 
linear differential equation. 

Proof for Theorem 2. The proof of theorem 2 is based on the technique of "reduction 
in order". Assuming that {Co( r ), B (r)} solves equation ()3.9j) . write 

Ci(r) = Co(r) Z (r) . (3.34) 

and demand that {(i(r),B (r)} also solves equation ()3.9|) . We find 

2r 2 £? (Co Z o y + (r 2 B' - 2rB )(Co Z )' + {rB' - 2B + 2) (Co Z Q ) = (3.35) 

we can expand the above equation to 

2r 2 B (Co Z + 2(' Z' + Co Z'£) + (r 2 B' - 2rB )(C Z + Co Z' ) 
+{rB' -2B + 2)(( Z )=0 (3.36) 

we can also re-group this same equation as 

{2r 2 B Co + (r 2 B' - 2r£? )Co + (rB' - 2B + 2)Co} Z 

+ (r 2 C 5 + 4r 2 C B - 2r( B )Z' + {2r 2 C, Q B Q )Z^ = (3.37) 

This linear homogeneous 2nd order ODE for Z now simplifies to 

(r 2 C i?o + 4r 2 C 5 - 1r^B Q )Z' G + (2r\ Q B Q )Z'^ = , (3.38) 

which is an ordinary homogeneous second-order differential equation, depending only 
on Z' Q and Zq. (So it can be viewed as a first-order homogeneous order differential 
equation in Z', which is solvable.) Separating the unknown variable to one side, 

IHfHH- ^ 
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Re-write Zq/Z' q = dln(ZQ)/di, and integrate twice over both sides of equation (|3.39j) . 
to obtain 

/rdr . 
Co(r) 2 V 5 o(r) 

depending on the old solution {( (r) , B (r)} , and two arbitrary integration constants 
a and e. 

To see that the transformation T 2 defined in Theorem 2 is "idempotent" we first 
show 

T 2 oT 2 ^T 2 , (3.41) 
and then iterate. The precise two-step composition rule is 

T 2 (<7 2 , e 2 ) o T 2 (cti, ei) = T 2 (a 2 vi, ^a 2 + ■ (3.42) 

In this section we consider a two-step composition rule to find idempotence. In 
particular, the most outstanding features of these steps are: 

• To see "idempotence", note that for fixed B (r) equation (|3.9|) has a solution 
space that is exactly two dimensional. 

• Since the first application of T 2 takes any specific solution and maps it into the 
full two-dimensional solution space, any subsequent application of T 2 can do 
no more than move one around inside this two dimensional solution space - 
physically this corresponds to a relabelling of parameters describing the perfect 
fluid metric you are dealing with, not the generation of new solutions. 

To be more explicit about this note that at step one 

Co -> Ci = Co { <ri + ei / r T , } , (3.43) 

Co(r) 2 ' 

while at the second step 
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That is: 



C2 = Co < 0"i + e x 



x < cr 2 + e 2 



r dr 



Co(r) 2 v/^oW 



(3.45) 



r dr 



Co(r) 2 v^BoTO ki + ei Jr dr/(Co(r) 2 v 7 ^))? 



But this can be rewritten as 



Co ^0 ^2 



1 1 

Z (7i 



Therefore 



That is 





, 62 


1 " 




Hi* 


0" 2 H 




^0} 






°1_ 



£1 



cr 2 + 



0-2(7! + I a 2 ei 



£2 1 

£2 



r dr 



Co(r) 2 v / A^) 



from which the composition law 



T 2 (cr 2 , e 2 ) o T 2 ((Ti, ei) = T 2 a 2 a\, ^2 + 



£2 
01 



(3.46) 
(3.47) 

(3.48) 
(3.49) 



(3.50) 



follows immediately. (Note that the composition law for T 2 is actually a statement 
about applying reduction of order to 2nd-order ODEs, it is not specifically a statement 
about perfect fluid spheres, though that is how we will use it in the thesis). The 
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general composition law then follows by induction. To find the inverse transformation 
we choose <j 2 = and e\ = — e 2 , for which 

T 2 (1/ct 1 , -e x ) o J 2 (a u e,) = T 2 (1, 0) = I. (3.51) 

□ 

Comment: As other special cases of the composition law we will also mention the 
results that 

n / n \ 

J]T 2 (l,e i ) = T 2 ll,^e i ), (3.52) 

i=l \ i=l / 

and 

T 2 ((7, e) n = T 2 (a 11 , e [a n ~ l + a 11 " 3 ■■■ + a~ (n - 3) + a^ n ^] ) . (3.53) 
Now as long as a > 1 then for sufficiently large n we see 

T 2 (a, e)» « T 2 (a", a' 1 -^) = a"" 1 T 2 (a, e) 4 T 2 (a, e), (3.54) 

where at the last step we have used the fact that the overall multiplicative factor cr"" 1 
can simply be reabsorbed into a redefinition of the time coordinate. Because of this 
result, we see that for fixed a > 1 (and fixed but arbitrary e) repeated numerical 
applications of T 2 (cx, e) will have a well-defined limit. In figure EOl we have tested the 
composition law numerically. 

A strictly limited version of theorem 2, with little comment on its importance, can 
be found in (2D]- Specifically, after several manipulations, changes of notation, and a 
change of coordinate system, the transformation exhibited in equation (16.12) of [20] 
can be seen to be equivalent to the sub-case a = 0, e = 1 of theorem 2 above. 

For some purposes it is more useful to rephrase theorem 2 as below: 

Corollary 2. Let {( a , B } and {0,, B } both represent perfect fluid spheres, then for 
all p and q 

{p(a + q(b,B } (3.55) 

is also a perfect fluid sphere. Furthermore, for fixed B all perfect fluid spheres can be 
written in this form. 
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0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 



Figure 3.4: The solid lines show £(r) for 10 reapplications of Theorem 2 onto the Minkowski 
metric. The dashed line corresponds to Si single application with a specific choice for 0" O nce 

and 

e once . It can be seen that 10 applications of Theorem 2 can be expressed by one application. 



Proof. The result is automatic once you note that for fixed Bq the 2nd order ODE for 
£ is linear and homogeneous. □ 

We defer extensive discussion about the application of these theorems and corollary 
until section EOl which cover the following: 

• A simple example involving T 2 applied to either the Einstein static universe or 
the anti-de Sitter universe results in Schwarzschild's stellar solution (position- 
independent density). 

• Similarly, corollary 2 applied to a combination of the Einstein static universe and 
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anti-de Sitter space is another way of obtaining Schwarzschild's stellar solution. 

Note that when comparing the results of our computations with the extant literature 
we shall adopt the naming conventions of the review article by Delgaty and Lake j^]. 
For example, some of the less standard names are: 

• K-0 III: The third spacetime considered by Korkina and Orlyanskii. 

• M-W III: The third spacetime considered by Matese and Whitman. 

• B-L: The spacetime considered by Buchdahl and Land. 

• Kuch: One of several spacetimes considered by Kuchowicz. 

• Heint: One of several spacetime considered by Heintmann. 

References to the primary literature can be found in Delgaty and Lake [§]. We also 
add several new abbreviations: 

• Martin: One of several spacetimes considered by Martin and Visser. 

• SI: A special name for a particularly simple spacetime. 

• P1-P8: Special names for several apparently new perfect fluid spacetimes. 

See in particular tables (3.2) and (3.3) for the spacetime metrics corresponding to 
these geometries. 

3.2.3 The third and fourth theorems 

Having now found the first and second generating theorems, is now useful to define 
two new theorems by composing them. What we are going to show in the following 
section is how to deform the geometry while still retaining the perfect fluid property. 
To further illustrate the formalism, we will show how to classify perfect fluid spheres 
by using two new generating theorems. We present theorem 3 by taking a perfect fluid 
sphere solution {(q, Bq} and applying theorem 1 onto it which gives us a new perfect 
fluid sphere {Co, B\\. The new B\ is given in equation (|3.11j) . We now continue with 



78 



seed metric:Minkowski 



old solution: 
K-0 III 



Figure 3.5: This structure shows that when we apply theorem 2 to Minkowski, we get the 
K-0 III geometry. 

applying theorem 2, again we get a new solution {£, B{\, where ( now depends on 
the new B\. All together we can consider this as a single process by introducing the 
following theorem: 



Theorem 3. If {Co, -Bo} denotes a perfect fluid sphere, then for all a, e, and X, the 
three-parameters geometry defined by 

, T dr 1 dt 2 + . . dr ' . - + rW 

Co(r) 2 V^o(r) + A A (r) J B (r) + AA (r) 

(3.56) 

a/so a perfect fluid sphere, where A is 

CoW V , f /"aw Co(r)-rC (r) 



Ao(r)= 7 ' r^exp 2 ^ " ^ dr . (3.57) 

7 5 = 7 2 o7, :{C 0)J Bo}^{Co,5 + AA (Co)} 

' * {Co^o(Co,5 + AA (Co)), J Bo + AA (Co)}- (3.58) 

Now we would like to present another new theorem. Instead of starting with the- 
orem 1 we could first apply theorem 2 on {(q,Bq}. This gives us a new perfect fluid 
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sphere B }, where Ci = Co Z (( Q , B ) is given by equation (|3.40|) . We now continue 
with theorem 1 which leads to {£1, B} where B depends on the new £i- Again, we can 
consider this as a single process by introducing the following theorem: 



Theorem 4. If {( , B } denotes a perfect fluid sphere, then for all a, e, and X, the 
three-parameters geometry are defined by 

ds 2 = -Co(r) 2 I a + e / T % 1 dt 2 + - - , - , + r W (3.59) 

z's a/so a perfect fluid sphere, where Ao(Ci, is defined as 

)+rCi(r-)/ 
depending on (i = Z Q , where as before 



A (G,r)= ^ v ' r ex Pi 2 / 7TT7TT— ^H dr ' ( 3 ' 6 °) 

Ci(r) +r CiW/ I J Ci(H CiW + r Ci( r ) 



Z (r) = a + el r ^ (3.61) 

Co(r)2 



TTiat is 



T 4 = 7,o7 2 :{Co, J Bo}^{Co^o(Co, J B ), J Bo} 

^ {Co^o(Co,5o),5 + AA (Co^o(Co,5o))}. (3.62) 

In addition, we analyze what happens, if we apply theorems 3 or 4 more than once. 
Some general comments about both new transformation theorems: 

• Theorem 3 and theorem 4 are in general distinct, which can be traced back to 
the fact that theorem 1 and theorem 2 do not in general commute. 

• Theorem 3 and theorem 4 are in general not idempotent. That means when 
we apply theorem 3 more than once further distinct solutions will be obtained. 
Similarly for theorem 4. 
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In addition, these comments can be described in terms of the equations below. 

T 3 #T 4 ; T 3 oT 3 #T 3 ; T 4 o T 4 # T 4 . (3.63) 



The best way to verify this is to try a few specific examples. There may be some 
specific and isolated special metrics for which theorem 3 and theorem 4 happen to be 
degenerate, or idempotent, and finding such metrics is important for our classification 
programme. 



seed metric:Minkowski 



- v 



old solution: 
K-0 III 



Th1 



new solution: 
Martin 3 



old solution: 
Einstein static 



V 
V 



old solution: 
Schwar;. Int. 



Figure 3.6: This structure shows that Theorem 3 and theorem 4 are in general distinct. 
When we apply theorem 3 to Minkowski, we get the Schwarzschild Interior geometry, while 
applying theorem 4 to the Minkowski, we get Martin 3 . 



3.2.4 Formal properties of the generating theorems 

The solution generating theorems we have developed interact in a number of interest- 
ing ways and exhibit numerous formal properties that will be useful when classifying 
generic perfect fluid spheres. To start with, theorem 3 (T 3 ) and theorem 4 (T 4 ) can 
be expressed in terms of theorem 1 (Ti) and theorem 2 (T2): 

T 3 (g):=(J 2 oJ 1 )g = J 2 (J 1 (g)); 

(3.o4) 

Mg) := (Ti o T 2 ) # = T x (T 2 (<7)). 
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where g is a metric representing a perfect fluid sphere. Having this in mind, some of 
the new solutions generated by starting from some specific solution can be identified. 
For example: 

T 4 o Ti = Ti o T 2 o Ti = Ti o T 3 , (3.65) 

or 

T 3 o T 3 = T 2 o Ti o T 2 o Ti = T 2 o T 4 o T l . (3.66) 
The idempotence of Ti and T 2 in this formalism is: 



(Ti o Ti) g =Ti (Ti(y))4Tx(^) 
(T 2 oT 2 ) 9 =T 2 (T 2 (g))±J 2 (g) 



(3.67) 



Taken together, it is possible to simplify all formulae wherever Ti and T 2 appear next 
to each other more than once. The following examples should demonstrate how this 
works: 



(T 2 o T 3 ) g = (T 2 oT 2 oTi)5 = (T 2 oTi)j = T 3 (g) 
(Ti oT 4 )g = (Ti o Ti o T 2 ) g = (T 1 o T 2 ) g = T 4 (g) 



(3.68) 



and in the same way 



T 4 o T 3 = Ti o T 2 o T 2 o Ti = Ti o T 2 o Ti = Ti o T 3 = T 4 o Ti ; 



T 3 o T 4 = T 2 o Ti o Ti o T 2 = T 2 o Ti o T 2 = T 2 o T 4 = T 3 o T 2 . 



(3.69) 



These relationships can be used to structure the solution set generated starting from 
any particular perfect fluid sphere, and moreover to classify which metrics can be 
produced by our theorems, and which ones cannot. For example, the idempotence 
property of theorem 1 and theorem 2 enables us to divide the class of perfect fluid 
spheres into seed metrics and non-seed metrics. Seed metrics can never be generated 
by using one of the two theorems Ti or T 2 , while non-seed metrics are connected to 
other simpler metrics via one of these theorems. We formalize this in the following 
subsection. 
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3.2.5 Seed and non-seed metrics 

Definition (Seed metric): Take a metric g (or a parameterized class of metrics) 
and apply theorem 1 or theorem 2 on it. Three different cases are possible: 

• Each of the applications supplies us with a new solution. [T^g) ^ g ^ T 2 (g).} 
We define a metric with this behaviour as a seed metric. (We shall soon see 
several examples of this behaviour.) For example, figure 2.5 show that Minkowski 
is a seed metric. When we apply theorem 1 to it, we get Einstein static. While 
applying theorem 2 to Minkowski, we get K-0 III instead. 

• Only one of the applications supplies us with a new solution, while the other one 
gives us the same metric we started with. [T^g) = g or T2(g) = g.) These met- 
rics are non-seed metrics. (We shall soon see several examples of this behaviour.) 
For example, figure 2.5 show that Einstein static and K-0 III are both non-seed 
metric. The reason is only one of the applications gives us a new solution, while 
the other one still give us the same metric which are Einstein static and K-0 
III, respectively. 

• Both applications give us the same metric we started with. [Ti(g) = g = T2(g).] 
Metrics of this type are fixed points of the transformation process and we then 
also have T 3 (g) = g = T 4 (g). While we have encountered numerical examples 
that seem to exhibit this behaviour, we have no analytic proof for the existence 
of non-obvious fixed-point metrics. There is one obvious but not particularly 
useful example of a fixed point class of metrics. If we take the ODE in equation 
(|3.8|) . and write down its most general solution as a function of the arbitrary 
parameters C( r ); then any of our solution generating theorems applied to this 
most general solutions will at most move us around in the parameter space 
characterizing the most general solution — the most general solution of equation 
()3.8j) . or equivalently equation (j3.9J) . is thus an infinite-parameter fixed point of 
the generating theorems. But apart from this obvious example, it is unclear 
whether other fixed point classes of metric exist. 

Classifying seed and non-seed metrics are very important. In particular, the most 
outstanding features of this algorithm are: 
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• We developed a tool to generate new solutions for a perfect fluid sphere, which 
does not require to solve the Einstein equations. 

• We also established the relationships among the generating theorems. Before all 
metrics seem to have nothing more in common than presenting a perfect fluid 
sphere. 

• We invented the concept of a seed metric, which cannot be generated by our the- 
orems. Starting with such a metric, a class of non-seed metric can be produced. 

• This can be used to structure the solutions generated by a perfect fluid sphere 
and classify which metrics can and cannot be produced by our theorems. 



Examples of some seed and non-seed metrics 



Seed metric 


Non-Seed metric 


Minkowski 

Exterior Schwarzschild 
de Sitter 


B-L 

Interior Schwarzschild 
Einstein static 


Tolman IV (A = 0) 

M-W III 

Heint Ha (C=0) 
SI 


Tolman IV , Tolman V {n — +1) , Tolman VI 

K-0 III 

Heint Ha 

Wyman I la, Wyman III 



Table 3.1: This table shows several seed and non-seed metrics, which satisfied the definition 
of seed and non-seed metrics (see Table 3.2 for further details). 

We can nevertheless develop several formal lemmata regarding fixed-point metrics. 
For instance 

Lemma 1. Suppose we have a metric such that V a, e, A 

T 3 (a : e,\)g±g, (3.70) 
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where we recall that = denotes equality up to redefinition of parameters. Then in 
particular 

T 1 (\)g±g±T 2 (a,e)g, (3.71) 

and so 

T 4 (a 1 e,X)g^g. (3.72) 
Conversely, suppose we have a metric such that 

T 4 (a,e,X) g = g, (3.73) 

then 

HXjgAgAT^g, (3.74) 

and so 

T 3 (a,e,\)g±g. (3.75) 

Proof. Trivial, note that T 3 (cr, e, A = 0) = T2(cr, e) and T 3 (<j = 0, e = 0, A) = Ti(A). 
Then recall T 4 = Ti o T 2 . Similarly for the converse. □ 

Lemma 2. Suppose we have a metric g such that 

T 4 g=T 3 g (3.76) 

and then define g' by 

T 4 g±g'±T 3 g. (3.77) 

Then we have 

T 4 g'±g'±T 3 g' (3.78) 
so that g' is a 'fixed point" of both T 3 and T 4 . 
Proof. Note 

T 4 g = T 4 o T 3 g = T 1 o T 2 o Ti g = T x o T 3 g = J 1 o T 4 g = T 4 g = g (3.79) 
and similarly for T 3 : 

T 3 9 — T 3 o T 4 g = T 2 o Ti o T 2 g = T 2 o T 4 g = T 2 o T 3 g = T 3 g = g (3.80) 

□ 
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Several other formal theorems along these lines can be constructed, but these seem 
the most important results. 

Finally, among the formal properties enjoyed by the generating theorems, we men- 
tion the fact that theorems 3 and 4 are "conjugate" to each other in the following 
sense 

T 4 eT 1 oT 2 = T 1 oT 2 oT 1 o [Tx]- 1 = T x o T 3 o [Tx]- 1 , (3.81) 
and similarly (when the appropriate inverse [T 2 ] _1 exists) 

T 3 = T 2 o Ti = T 2 o Ti o T 2 o [T 2 ] _1 = T 2 o T 4 o [T 2 }~ 1 . (3.82) 

We can write this as 

T 3 (<T,e,\)~T A (a,e,\), (3.83) 

where ~ denotes the concept of "similarity" under conjugation by invertible generating 
theorems. 

In these properties we would like to ensure that 

• We wish to emphasise that similarity ~ is a statement that holds for particular 
and fixed values of the parameters (a, e, A), as opposed to = which denotes 
equivalence under redefinition of parameters. 

• If one is working numerically, it is much easier to ask questions involving sim- 
ilarity ~. For analytic work, it is typically easier to ask questions involving 
equivalence =. 

3.2.6 The fifth and sixth theorems 
Two linking theorems 

The last two solution generating theorems we shall present are slightly different from 
those developed so far: Consider a perfect fluid sphere solution {£ , -Bo} an d extend 
it to a new perfect fluid sphere {Co^o, B + A}. Previously, we had either set Z = 1 
and obtained theorem 1, or we had set A = and obtained theorem 2. In other words, 
we only changed one metric component at a time. (From this point of view theorem 3 
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and theorem 4 are, strictly speaking, not new theorems, in that they are replaceable 
by iterations of theorem 1 and theorem 2 and vice versa.) We now investigate what 
happens if we place no a priori restrictions on Z and A, and allow both metric com- 
ponents to vary simultaneously. The differential equation (|3.8J) for this problem now 
becomes 

[r«oV] A' + [2r 2 (Co V - 2(rCo^o)'] A 

+ [r 2 ( B' + Ar 2 C B - 2r( Q B Q ] Z' + 2r 2 ^B Z'^ = 0. (3.84) 

Note that if A = this becomes equation (J3.17|) . while if Z = 1 this becomes equation 
()3.38|) . In general, this ODE of first-order in A, and — as long Z is not a constant - 
inhomogeneous. In terms of A this ODE can be solved explicitly and the result stated 
as a new theorem: 

Theorem 5. Suppose {Co, A)} describes a perfect fluid sphere, and let Z (r) be arbi- 
trary. Define 



A(A,r) = A (r)l\-J 



n2 



[(Ar 2 C B + r 2 Co^ - 2rCo£? ) Z' + 2r\^Z>>\ {Co + 
r3 (rCo^o)' Co 2 Zl 



i (Co-^o)' Co^o — r (Co^o)' , 
a ^- 2 % (,Vr(W dr ' Jr '' (385) 



where 



Co^o V 2 f f (Co^o)' Co^o — r (Co^o)' i 
A o( r ) = U v , „ ^ >z v r ex P ^ 2 / /• v a ^ , : 77 v v dr t • ( 3 - 86 ) 



v Co^o + r (CoZo)'J " I J (oZo ( Z + r (CqZ )' 

Then for all X, the geometry defined by an arbitrary chosen -Zo(r) and setting 

ds > = -CMrY + flo(r) d ;l (Air) + (3.87) 
corresponds to a perfect fluid sphere. That is, the mapping 

T 5 (X) : {Co,B } i-> {Co Z ,B + A(A,Co)} (3.88) 
ia£;es perfect fluid spheres into perfect fluid spheres. 
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Note that if Z (r) = 1 this simply reduces to theorem 1. Re-arranging equation 
(|3.84|) in terms of Z leads to a second-order inhomogenous differential equation, which 
cannot in general be solved for a prescribed A, unless one imposes further constraints. 
So further exploration in that direction is moot. There is however a related theorem 
(which may be easier to understand) in terms of parametric derivatives: 

Theorem 6. Let {((fi), B(n)} denote a one-parameter class of perfect fluid spheres, 
so that the differential equation hS. <5j) is satisfied for all \i. Then 



|r(rCm _j +|2rr -2K)^-j + 2,.^ 



+ (r 2 B' - 2rB) (^) + (rB' - 2B + 2) = 0. (3.89) 



In particular if d(/dfi = this reduces to the ODE \3. 1 7| ), while if dB/dfi = this 
reduces to the ODE / EO)) . 

This is simply an alternative viewpoint on the previous theorem, emphasising the 
differential equation to be solved. 

We again defer extensive discussion to the next section, but that this stage point 
out that if we invoke theorem 5 and apply it to Minkowski space, then making the 
choice Zq = 1 + r 2 /a 2 leads to the general Tolman IV metric — that is: 

T 5 (Minkowski; Z = 1 + r 2 /a 2 ) = (Tolman IV). (3.90) 

Even before we systematically start our classification efforts, it is clear that the solution 
generating theorems we have established will inter-relate many of the standard perfect 
fluid spheres. 



3.3 Formal properties of the linking theorems 

Before turning to issues of systematic classification of perfect fluid metrics, we wish 
to establish a few formal properties of the linking theorems. To simplify the notation, 



N8 



seed tnetricMnkowski 
^ = 1 + 1*2/4*2 



Th5 



Tolrnan IV 



Figure 3.7: This structure shows that when we invoke theorem 5 and apply it to Minkowski 
space, then making the choice Zq = 1 + r 2 /a 2 leads to the general Tolman IV metric. 

let us define the differential expression 

D(C,B) = [r(rQ']B' + [2r 2 C," -2(rQ']B + 2t 

= 2r 2 B(" + (r 2 B' - 2rB)£' + (rB' - 2B + 2)(. (3.91) 

Then the condition for a perfect fluid sphere is simply 

D((,B) = Q (3.92) 

Now define 

D(Co, B + A ) = D(Co, B ) + L>i(Co; A ) (3.93) 
Then it is easy to check that 

Di(Co; A ) = [r(rCo)1 A' + [^Co ~ W] A (3.94) 
The ODE for theorem 1, where we assume {Co, -Bo} is a perfect fluid sphere, is then 

D 1 (Co;A )=0. (3.95) 

Now define 

D(( Z , B ) = Z D((o, B ) + L> 2 (Co, B ; Z Q ), (3.96) 
then it is easy to check that 

D 2 (Co, B ; Z ) = (r 2 ( B' + 4r 2 (' B - 2r( Bo)Z' + (2r 2 (oB )Z'J. (3-97) 
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The ODE for theorem 2, where we assume {(o, -Bo} is a perfect fluid sphere, is then 

D 2 (Co, B ; Z ) = 0. (3.98) 

Now let us consider any simultaneous shift in ( and B, as considered in theorem 5. 
We have 

D(( Z ,B + A ) = D(C Z ,B ) + D 1 (CoZ ;A ) (3.99) 
= Z D((o, B ) + D 2 ((o, B ; Z ) + £>i(Co^o; A ). 

But now let us write 

D l (( Z ; A ) = ZoD^Cq; A ) + £> 12 (Co; ^o, A ), (3.100) 
where a brief computation yields 

D l2 ((o; Z , A ) = r 2 Co^ A + [2r 2 Co^o + 4r 2 Co^o " 2rCo^ ] Ao- (3.101) 
Then all in all 

D(( Q Z Q ,B + A Q ) = Z Q D(Co,B ) + Z D 1 (Co;A Q ) 

+D l2 (Co; Z , A ) + £> 2 (Co, B ; Z ). (3.102) 

So if {Co, -Bo} an d {Co^o, -Bo + A } are both perfect fluid spheres we must have 

Z £>i(Co; A o) + DuiCo; Z , A ) + £> 2 (Co, B ; Z ) = 0. (3.103) 

This is 

Z x (ODE for theoreml) + (cross term) + (ODE for theorem2) = 0. (3.104) 

The cross term vanishes if either Z = constant or A = in which case we recover 
the usual theorem 1 and theorem 2. If we do things in the opposite order 

D(CoZ , Bo + A ) = Z D((o, B + A ) + £> 2 (Co, B + A ; Z ); (3.105) 
= Z D(Co, B ) + Z D 1 (Co, A ) + £> 2 (Co, B ° + A o; Z ). (3.106) 
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We now have to compute 

D 2 ((o, B + A ; Z ) = D 2 (( , B ; Z ) + An (Co, B ; Z , A ), (3.107) 

and after a brief computation 

D 21 (Co,So;^o,A ) = 2r 2 CoA ^ + [r 2 CoA , + 4r 2 ^A -2rCoAo]^ 

= D 12 (( Q ,B ;Z ,A ). (3.108) 

Thus the cross term is the same, no matter how you calculate it, and we still have the 
identity 

D(( Z ,B Q + A Q ) = Z D(( ,B Q ) + Z Q D l (( ;A ) 

+D 12 (Co; Z , A ) + D 2 (Co, B ; Z ). (3.109) 

Provided Co, B Q and {Co-^o, B + A } are both perfect fluid spheres we again deduce 

ZoD^Cq; A ) + D 12 (( ; Z , A ) + D 2 (( Q , B ; Z ) = 0. (3.110) 

Now this gives us another way of looking at theorem 3 and theorem 4. For theorem 3 
we first apply theorem 1 so we have the two equations 

J D 1 (Co;A )=0; (3.111) 

and 

£> 2 (Co, B + A ; Z ) = D 12 (C ; Z , A ) + £> 2 (Co, B ; Z ) = 0. (3.112) 
Conversely, for theorem 4 where we first apply theorem 2 we see that we need to solve 

D 2 (( 0l B ;Z ) = 0; (3.113) 

and 

D l (C Z ; A ) ee Z D l (C ; A ) + D 12 (( ; Z , A ) = 0. (3.114) 
For theorem 5 we pick Zq arbitrarily, and solve the single ODE 

D 1 (CoZ ; A ) + D 2 (Co, B ; Z ) = 0. (3.115) 
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This is a single first-order linear inhomogeneous ODE for A , and hence solvable. (In 
particular this makes it clear that theorem 5 is an inhomogeneous version of theorem 
1 with a carefully arranged "source term" -D 2 (£ , B Q ; Z ). While theorem 5 is not 
"idempotent" it does satisfy the important formal property that: 

Lemma 3. 

T 1 oT 5 ±T 5 (3.116) 

which in particular tells us that the output from theorem 5 is never a seed metric. 

Proof. Applying theorem 5 we need to solve 

D 1 (( Z ; A ) + £> 2 (Co, B ; Z ) = (3.117) 

in order to map 

{( ,B }^{( Z ,B + A } (3.118) 

Now apply theorem 1 to {(oZq, Bq + Aq}, this means we have to solve the homogeneous 
ODE 

ACCo^Ai) =0 (3.119) 
But then, using properties of first-order ODEs 

{£>i(Co^o; A ) + £> 2 (Co, B ; Z ) = 0} © {L>i(Co^o; Ax) = 0} (3.120) 

=> {D^CoZo; A + Ai) + D 2 (( , B ; Z ) = 0} (3.121) 
which is the ODE from theorem 5 back again. □ 

(The net result of this observation, as we shall see in the next section, is that 
theorem 5 can be used to connect one seed metric with the "descendants" of another 
seed metric.) 
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3.4 Classifying perfect fluid spheres 



We will now see the power of these transformation theorems (solution generating 
theorems) by using them in a number of different ways: 

• To generate several new perfect fluid spheres 

• To relate various perfect previously known fluid spheres to each other. 

• And also to classify the geometries we encounter. 

First some minor warnings: Despite comments made in jU], Kuch 2 1 = Tolman V; 
and Kuchb I b is a perfect fluid for general values of its parameters. Furthermore 
RR I = Einstein static; RR-V = Tolman V (n = —5/4). If we had not noted these 
degeneracies then our tables below would have been more extensive, but would have 
conveyed no extra information. 

Starting with the metric for any known perfect fluid sphere and successively ap- 
plying theorem 1 and theorem 2 numerous times will supply us with endless "new" 
perfect fluid sphere solutions. There are so many types of "new" solutions such as: 

• Some of these "new" solutions might already be findable in the literature, some 
of them might be truly novel. 

• Some of these solutions can be written down in a fully explicit form. 

• Some solutions are explicit but not elementary, in the sense that while the metric 
components can be exhibited as specific and explicit integrals, these integrals 
cannot be done in elementary form. 

• Some solutions are so complex that present day symbolic manipulation programs 
quickly bog down. 

(For specific symbolic computations we have used a vanilla installation of Maple. See 
appendix C for some sample worksheets). 

We present several tables and diagrams to summarize the situation. Two tables 
are used to provide the names and explicit metrics for the perfect fluid spheres we 
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consider. Two other tables are used to describe the inter-relationships of these perfect 
fluid spheres under T 1; T2, T 3 , and T 4 . Recall that a seed metric is one for which 
theorem 1 and theorem 2 both yield metrics distinct from the seed: T 1 (g) ^ g ^ T2(g). 
In contrast for non-seed metrics one or the other of these theorems is trivial, either 
Ti(#) = g or T 2 (g) = g. 

In these tables, we can see the notation "[integral]" means that it is definitely a 
novel perfect fluid solution. However, the metric components involve an explicit inte- 
gral that does not appear to be do-able by elementary methods, and is so complicated 
that it does not seem worthwhile to even write it down. 

By considering that theorem 1 and theorem 2 have "idempotent" properties , and 
the fact that theorem 3 and theorem 4 can be expressed in terms of the first two 
theorems, it is possible to structure, and therefore to graphically visualize the rela- 
tionship between all metrics generated from a given seed metric. We demonstrate 
this behaviour starting with Minkowski spacetime as seed metric where we are start- 
ing with this trivial seed to check if it is possible to create endless new perfect fluid 
spheres. The first few steps can be carried out explicitly, and show that Minkowski 
space generates several well known interesting perfect fluid models. 

In figure I3~%1 all the ideas from the previous sections are used. Each box represents a 
specific metric (perfect fluid sphere) while the arrows correspond to the application of 
the different theorems. The horizontal arrows correspond to an application of theorem 
1 and the vertical arrows to an application of theorem 2. In addition, the vectors 
pointing along the diagonals can either indicate an application of theorem 3 or theorem 
4. A dotted arrow corresponds to the application of theorem 3 while the dashed arrow 
represents an application of theorem 4. 

Figure 13.21 shows that after applying theorem 1 to the Minkowski seed metric, we 
get the Einstein static universe. By the idempotence of theorem 1 n applications of 
Ti to the seed metric still results in the Einstein static. Similary, any number of 
applications of theorem 2 after the first (which leads to the K-0 III solution) does 
not give us any further new solutions (see figure 13. 4j) so we can see form this figure 
that theorem 2 is "idempotent" . In addition, we can also see that the first and second 
generating theorems are not commutative. Application of theorem 1 and theorem 2 in 



91 



seed metric:Minkowski 



Th1 



□ Id solution: 
Einstein static 
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□ Id solution: 
Schwarz. Int. 



Th1 



old solution: 
K-0 III 



Th1 



new solution: 
Martin 3 



new solution: 
P1 



Th1 



new solution 
[ [integral] J 



new solution 
[ [integral] J 



new solution 
[integral] j 



new solution 
[integral] j 



Th1 



. new solution: 
, [integral] j 



Figure 3.8: Structure graph for Minkowski space used as seed metric. 



that order to the Minkowski seed metric results in the Schwarzschild Interior geometry, 
whereas application of theorem 1 after theorem 2 gives us the Martin 3 solution, which 
is quite distinct from the three-parameter Schwarzschild interior solution. Indeed, it 
seems as if it is possible to create endless new solutions for a perfect fluid sphere out 
of the Minkowski metric (or any other of the seed metrics). After several iterations 
the calculations become more complex, and can no longer be carried out analytically. 
We then resort to numerical computation to find out whether theorem 3 and theorem 
4 have some sort of numerical limit, a numerical fixed point, or not. Depending on 
the choice for (A, a, e) both theorems converge very quickly. Specifically, we used 
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Matlab to numerically analyze the evolution of £(r) and B(r) after applying theorem 3 
several times to the Minkowski metric. Numerical results are summarized in figure 13*31 
This figure indicates that theorem 3 and theorem 4 both appear to have well defined 
numerical limits, though we have no direct analytic solution for the fixed point metric. 



seed metric:Minkowski 



Th1 



old solution: 
Einstein static 



old solution: 
K-0 III 



Till 



old solution: 
Schwarz. Int. 



Th1 



new solution: 
Martin 3 



new solution: 
PI 



Thl 



new solution 
[ [integral] J 



new solution: 
, [integral] 



new solution 
[integral] j 



new solution: 
[integral] 



Th1 



new solution: 
. [integral] } 



fixed point: 
Theorem 3 



fixed point: 
Theorem 4 



Figure 3.9: Numerical implementations of Theorem 3 and Theorem 4 can be used to 
investigate the behavior for n applications, if the number of applications goes to infinity. 
For specific choices of parameters Aj, (Xj, and e^, multiple reapplications of Theorem 3 and 
Theorem 4 approach a limit, in the sense that both metric components are converging towards 
there fixed-points. 
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3.5 Discussion 



We have developed several transformation theorems that map perfect fluid spheres to 
perfect fluid spheres using Schwarzschild coordinates. We used these transformations 
as a basis for classifying different types of perfect fluid sphere solutions. Applying 
these theorems on a known perfect fluid sphere, different solutions are often obtained. 
Some of the solutions are already known in the literature but most of them are novel. 
Moreover, we can classify which metrics can be produced by our theorems, and which 
ones cannot. The first theorem we presented is in a slightly different formalism, and 
demonstrates an independent way of proving it. Indeed, the second theorem we pre- 
sented is a new transformation theorem. The idempotence property of Theorem 1 and 
Theorem 2 can divide the class of perfect fluid spheres into seed and non-seed met- 
rics. Classifying seed and non-seed metrics are important because we have developed 
a tool to generate new solutions for a perfect fluid sphere, which does not require us 
to directly solve the Einstein equations. Indeed, the whole procedure was set up in 
such a way that we are implicitly and indirectly solving the Einstein equations, but 
the utility of the transformation theorems is that one does not have to go back to first 
principles for each new calculation. 

We established two new theorems by composing first and second generating theo- 
rems. They are in general distinct, which can be traced back to the fact that Theorem 
1 and Theorem 2 do not, in general, commute. Indeed, Theorem 3 and Theorem 4 
are, in general, not idempotent. That means when we apply Theorem 3 or Theorem 
4 more than once, further solutions will in general be obtained. In addition, we can 
see the power of these transformation theorems as we can generate several new perfect 
fluid spheres and also can classify the geometries we encountered. 

Furthermore, we have also established several relationships among the generating 
theorems. Previously, all metrics seemed to have nothing more in common than pre- 
senting a perfect fluid sphere. We have developed the concept of a seed metric, which 
is one that cannot be generated by our theorems. Starting with a seed metric and 
applying our theorems, it is possible to structure, therefore to visualize, the relation- 
ship between all metrics generated from a given seed metric in a graph. Based on this 
example it is possible to create endless new solutions out of the Minkowski metric. We 
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also used a numeric program to investigate whether Theorem 3 and Theorem 4 have 
fixed point limits. Both seem to converge very quickly. The situation is summarized 
graphically in figures (3.1) to (3.9), and in tabular form in tables (3.1) to (3.5). 
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Some selected perfect fluid solutions 



Name 
Minkowski 
Einstein static 
de Sitter 

Schwarzschild Interior 
Schwarzschild Exterior 



Metric 

-dt 2 + dr 2 + r 2 dn 2 

-dt 2 + (1 - r 2 /E?)- l dr 2 + r 2 dQ 2 

-(1 - r 2 /R 2 ) dt 2 + (1 - P/R 2 )- 1 dr 2 



A- B^Jl-r 2 /R 2 ) dt 2 + {l-r 2 /R 2 )- 1 dr 2 + r 2 dtt 2 
(1 - 2m/r) dt 2 + (1 - 2m/r)- 1 dr 2 + r 2 dQ 2 



SI 

K-0 III 
Kuchl lb 
B-L 



-r 4 dt 2 + dr 2 + r 2 dQ 2 

-(A + Br 2 ) 2 dt 2 + dr 2 + r 2 dVL 2 

-{Ar + Br\nr) 2 dt 2 + 2dr 2 + r 2 dfi 2 

-A(r 2 /a 2 )dt 2 + 2(1 + r 2 /a 2 )- x dr 2 + r 2 dfi 2 



Tolman IV 
Tolman IV (b -> oo) 



^ 2 (i + . 2 / a 2 )dt 2 + (1 _^y; 

l+2r 2 /a 2 



-B 2 (1 + r 2 /a 2 ) dt 2 + dr 2 + r 2 df! 



dr 2 + r 2 dn 2 

2 



Tolman V 
Tolman V (A -> 0) 
Tolman VI 

Tolman VII 



Tolman VIII 



_£2 r 2(l+„) dt 2 + (2 _ n 2) [! _ i4r 2(2-n")/(2+n)]-l dr a + r 2 d ft2 
_ B 2 r 2(l+„) d< 2 + (2 _ n 2 )dr 2 + r 2 df7 2 

-(Ar 1 "" + Br 1+n ) 2 dt 2 + (2 - n 2 )dr 2 



2 dfl 2 



B 2 cos < In 



^/l-2r 2 /a 2 +r i /b i + r 2 /a 2 - b 2 /a 2 



1/2 



+ (1 - 2r 2 /a 2 + r 4 /6 4 )-Mr 2 + r 2 dfi 2 



_^2 r 2(n-l)(>-4)/n 



(n 2 -4n+2)(ra 2 -8n+8) 
(n 2 -4n+2) (ri 2 -8n+8) 



_|_ ^ r -(n 2 -8n+8)/n _|_ £y-2(n 2 -4n+2)/n^ 



+ 5 r -(™ 2 -8ri+8)/n _|_ (j r -2{n 2 -4n+2)/n~j dr _2 _|_ , 



2 dr> 2 



- (1 - 2m/r) dt 2 + [(1 - 2m/r) (l + C(2r - 2m) 2 )] ~ X dr 2 + r 2 dQ 2 

- (^A^/l - 2m/r + B r 2 /m 2 + 5r/m - 30 + 15^1 - 2m/rln{l - r/m - y/r(r - 2m) /m} 

+ (1 - 2m/r)- 1 dr 2 +r 2 dVL 2 
-Ar(r -a)+ x ]^ jcii dr 2 + r 2 dtt 2 

-r 2 [A + B arctanh(a/Va 2 + r 2 )] 2 dt 2 + 2(1 + r 2 /a 2 ) _1 dr 2 + r 2 dn 2 
-(1 + ar 2 ) 3 dt 2 + [(1 + ar 2 )/(l - ar 2 /2)]dr 2 + rW 



Kuch 68 II 
Kuch 68 I 

M-W III 
Kuch I b 
Heint Ila (C=0) 

Hcint Ila 



df 



-(l + ar 2 ) 3 dt 2 + 



1 - 



3ar 2 



2(l+ar 2 ) + (l+ar 2 )Vl+4ar 2 



Cr 2 



dr 2 + rW 



Table 3.2: Some well-known perfect fluid spheres and their coresponding metrics. Note 
that we have often reparameterized these metrics to make them easier to deal with, and so 
their appearance (but not the substance) may differ from other sources pi llllOU] . Metric SI 
is a special case of K-0 III, Tolman V, and Tolman VI, notable for its extreme simplicity. 
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Some apparently new perfect fluid solutions 



Name 



Metric 



Martin 1 
Martin 2 
Martin 3 



-(Ar + Brlnr) 2 dt 2 + 2 ^f± 2 f^g^ dr 2 + rW 

-Ar(r - «)^ 2 + I (l - £ - ' dr 2 + r W 

-(1 + ar 2 ) 2 dt 2 + [1 - 6r 2 /(l + Sar 2 ) 2 / 3 }' 1 ^ 2 + r 2 dfl 2 



PI 
P2 

P3 

P4 
P4 

P5 
PG 
P7 
P8 



(1 + ar 2 ) 2 



(Ar + Br Inr) 2 



B J —j= 



A /l_h r 2/( 1 + 3ar 2)2/3 



rdr 



d£ 2 + [1 - 6r 2 /(l + Sar 2 ) 2 / 3 ]-^^ + r 2 dft 2 

-i 2 



a + e / (Ar + Br In r)" 2 rdr 



-1/2 



(If 2 



O 2A+2B In r+B j 2 , „2jn2 
" Z 2A+2Blnr+B-Cr^ Clr 

p (5+2ar 2 )^l-ar 2 /2 

A-I--D ( 1 + ar 2)3/2 



(1 + ar 2 ) 3 



di 2 + [(1 + ar 2 )/(l - ar 2 /2)]dr 2 + r 2 d0 2 
^(a + BJ — \ dt 2 + (1 + Ar 2 / 3 )- x dr 2 + r 2 dfi 2 

V r 3 -y/l+Ar 2 / 3 / 

r 4 (A + S [ifA^anh-^l/Vl + Ar 2 / 3 ) - ^ Vl + Ar 2 / 3 (8r~ 2 - lOAr" 4 / 3 + 15A 2 r- 2 / 3 )]) 
+ (l + Ar 2 / 3 )- 1 dr 2 + r 2 dfi 2 



dt 2 



(A^/7~ a + B VrTa) + ^.V^ dr 2 + 



-r(r-a) \A + B J 



dr 



(r-a), 1-5--B 



-air 7 / 3 



7 A r 2 B (r-aV 3 \ 1 j 2 
+ 4 1/ -° (4r-3o)*/3 J ar + 



(4r- 3a) 4 /3 



B 2 (l + r 2 /a 2 ) A + B/ ^g^ ^rdr 



di 2 



l+2r 2 /a 2 



dr 2 + r 2 dft 2 



(1 + ar 2 ) 3 



A + Bf- 



r dr 



(1+ar 2 ) 2 ./! 



(l-r 2 /6 2 )(l+r 2 /a 2 ) 
2 



dt 2 



- + 



2(l + a r 2) (1 + ar 2 )v / 1 + lar 2 



CV 2 



2(l+ar 2 ) (l+ar 2 )Vl+4ar 2 



dr 2 + r 2 dft 2 



Table 3.3: Some apparently new prefect fluid spheres and their coresponding metrics. 
Sometimes the relevant integrals cannot be done in elementary form. When they can be 
done they are explicitly shown. 
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Some seed geometries and their descendants 



Seed 


Theorem 1 


Theorem2 


Theorem 3 


Theorem 4 


Minkowski 


Einstein static 


K-0 III 


interior Schwarzschild 


Martin 3 


exterior Schwarzschild 


Kuch68 II 


Kuch 86 I 


[integral] 


[integral] 


de Sitter 


Tolman IV 


interior Schwarzschild 


P7 


[integral] 


Tolman V (A = 0) 


Tolman V 


Tolman VI 


Wyman III 


Wyman Ha 


SI 


Tolman V (n = +1) 


K-0 III 


P4 


Martin 3 


M-W III 


Martin 2 


P5 


P6 


[integral] 


Heint Ha (C=0) 


Heint Ha 


P3 


P8 


[integral] 



Table 3.4: Seed solutions and their generalizations derived via theorems 1-4. The notation 
"[integral]" denotes a metric so complicated that explicitly writing out the relevant integral 
is so tedious that it does not seem worthwhile. 



Some non-seed perfect fluid geometries and their descendants 



Base 


Theorem 1 


Theorem2 


Theorem 3 


Theorem 4 


Tolman IV 
B-L 

Heint I la 


Tolman IV 
B-L 

Heint Ha 


P7 

Kuchb I b 
P8 


P7 

Kuchb I b 
P8 


[integral] 
[integral] 
[integral] 


Tolman VI 
Kuchl lb 
K-0 III 


Wyman I la 
Martin 1 
Martin 3 


Tolman VI 
Kuchl lb 
K-0 III 


[integral] 

P2 

PI 


Wyman I la 
Martin 1 
Martin 3 



Table 3.5: Non-seed solutions and their generalizations. 
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Chapter 4 



Solution generating theorems — 
other coordinates 

Consider other coordinate systems: This chapter generalizes the theorems which we 
derived for Schwarzschild coordinates to a number of other coordinate systems. This 
may allow us to develop extra relations between the known solutions. 

In all cases we give the metric and the pressure isotropy condition G?? = Gg§. 
Finch and Skea ^T] estimate that about 55% of all work on fluid spheres is carried out 
in Schwarzschild coordinates, that isotropic coordinates account for 35% of research, 
and that the remaining 10% is spread over multiple special cases (see below). 



4.1 General diagonal coordinates: 

Consider the metric: 

dv 2 

ds 2 = -C(r) 2 dt 2 + ^— + R(r) 2 dQ 2 (4.1) 
B(r) 

and then assume it represents a perfect fluid sphere. That is, Gg§ = G?? = G^. The 
equation = Gi^ is always fulfilled due to spherical symmetry. 
Calculate 

_ CB(R') 2 + 2BRCR' - C (A0 , 
103 



and 

1 2B('R + 2(R"B + (B'R + 2R( "B + R( 'B' 
Gee ~ "2 j% ( 4 - 3 ) 

We set Gff = Gq§, which supplies us with an ODE. This ODE reduces the freedom to 
choose the three functions in equation (j4.1J) to two: 

[R(R()']B' + [2RR"( + 2R 2 (" - 2RR'(' - 2(R') 2 QB + 2( = 0. (4.4) 

This is first order-linear non- homogeneous in B(r), and second-order linear homoge- 
neous in ((r). Note in particular it is linear in both £(r) and B(r) and the comment 
below equation (16.10) in "exact solutions" [20] (while true) completely misses the 
mark. There is nothing particularly nice about the R(r) dependence. 



4.2 Schwarzschild (curvature) coordinates: 

This is the basis of the theorems which were previously derived above. 
Consider the metric 

ds 2 = -C(rf dt 2 + ^-+ rW (4.5) 
B(r) 

and the ODE 

[r(r()']B' + [2r 2 C" - 2«)']£ + 2( = 0. (4.6) 

This is first order-linear non- homogeneous in B(r), and second-order linear homoge- 
neous in C( r )- Solving for B(r) in terms of ((r) is the basis of the [IS], (and is the 
basis for Theorem 2 in chapter 3). 

Sometimes people prefer to write B(r) in terms of m(r) 

ds 2 = -C(r) 2 dt 2 + ^-T + r W ( 4 - 7 ) 

j _ 2m(r) 

or even /i(r) = m/r 3 , but there is no qualitative gain in doing so. 
If we re-group things in terms of ( we find 

2r 2 C" + (r 2 B' - 2rB)(' + (rB' -2B + 2)( = 0, (4.8) 
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which is a linear homogeneous 2nd order ODE. This is the basis of Theorems 1 and 3 
in the above section. 



4.3 Isotropic coordinates: 

Consider the metric 

d. 2 = -C(r) 2 dt 2 + c(r) 4 (r)2 ^ 2 + rW } ( 4 - 9 ) 

Calculate 

G f r = (C) 2 B 2 - (B') 2 C 2 + IB'BCjr (4.10) 

and 

G §§ = -(B'fC 2 + ( 2 BB" - {C'fB 2 + B'BCjr (4.11) 

We set Grf = G^, which supplies us with an ODE. This ODE reduces the freedom to 
choose the two functions in equation ()4.9j) to one: 

() 2B y ' 

There are several ways of improving this. If we write £(r) = exp( J g(r)dr) then we 
have an algebraic equation for g(r): 



^) = ±\I B " J''" - ( 4 - 13 ) 
Conversely, the isotropy condition can be written in terms of B(r) as: 

B" - B'/r - 2g 2 B = 0. (4.14) 

This is about as simple as you can make things. 

There is an improvement in writing B(r) = exp(2 J h(r)dr) so that 

g{rf = 2/i(r) 2 + h'{r) - h{r) /r (4.15) 

which is the basis of |18j . 
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4.3.1 The seventh and eighth theorems 

Theorem 7. In Isotropic coordinates, if {((r), B(r)} describe a perfect fluid then so 
does {C( r ) _1 > B(r)}. This is the Buchdahl transformation in disguise. 
That is, the geometry defined by holding B (r) fixed and setting 

ds 2 = dt 2 + ^rr^dr 2 + r W} (4.16) 

Coin B Q {ry 

is also a perfect fluid sphere. Alternatively, the mapping 

7>:{Co,5o}^{Co~\£o} (4.17) 
takes perfect fluid spheres into perfect fluid spheres. 

Proof for Theorem 1: Assuming that {Co( r ) _ \ Bo(r)} solves equation (j4.12j) . 
We consider L.H.S. of the equation ()4.12|) . 

— (i)WW^HfHf) 2 ™ 

This shows that {((r), B(r)} describe a perfect fluid then so does {((r)" 1 , B(r)}. 

The transformation T 7 defined in Theorem 7 is a "square root of unity" in the 
sense that: 

T 7 o T 7 = /, (4.19) 
To see this, consider the sequence 

{Co, B } » {Ci, B } i * {( 2 , B }... (4.20) 

In this section we consider a two-step composition to check that T 7 is a square root of 
unity. In particular, to be more explicit about this note that at step one 

Ci = Co" 1 (4-21) 

while at the second step 

C 2 = Cr 1 = (Co" 1 )" 1 = Co (4.22) 
For reasons of convenience and without loss of generality, we can write for C( r ) ; after 
applying thorem 1 n-times, as : ((r) = C( r ) - \ for n — odd number and £(r) = C( r ); 
for n = even number. □ 
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Theorem 8. Let {( (r) , B (r)} describe a perfect fluid sphere. Define 

Then for all a and e, the geometry defined by holding Co( r ) fixed and setting 

ds 2 = -Co(r) 2 dt 2 + ^ 2R \ 27M2 {dr 2 + rW} (4.24) 
Co{r) 2 B (r) 2 Z (r) 2 

is also a perfect fluid sphere. That is, the mapping 

T 8 (a, e) : {Co, B } >-> {Co, B Z (B )} (4.25) 
takes perfect fluid spheres into perfect fluid spheres. 

Proof for Theorem 8. The proof of theorem 8 is based on the technique of "reduction 
in order". Assuming that {Co( r ), Bo(r)} solves equation (|4.14j) . write 

B x (r) = B {r) Z (r) . (4.26) 

and demand that {Co( r ), Bi{r)} also solves equation ()4.14|) . We find 

(B Z )" - (B Zo)'/r - 2g 2 (B Z ) = 0. (4.27) 

we can expand the above equation to 

{B» Z + 2B' Z' Q + B ZD - (B'o Z )/r - (B Z' )/r - 2g 2 (B Z ) = 0. (4.28) 

we can also re-group this same equation as 

{B' ' - B'Jr - 2g 2 B ] Z + 2B' Z' + B Q Z' ' - (B Q Z' )/r = 0. (4.29) 

A linear homogeneous 2nd order ODE for B now simplifies to 

(2B' -Bo/r)Z'o + B Zo" = 0, (4.30) 

which is an ordinary homogeneous second-order differential equation, depending only 
on Z' Q and Z' Q '. (So it can be viewed as a first-order homogeneous order differential 
equation in Z' ', which is solvable.) Separating the unknown variable to one side, 

^L = -2^ + -. (4.31) 
Z' Bo r { } 
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Re-write Zq/Z' q = dln(ZQ)/dr, and integrate twice over both sides of equation (|4.31jl . 
to obtain 

Zo(r) = f + e /w}' (4 ' 32) 

depending on the old solution {( (r), B (r)}, and two arbitrary integration constants 
a and e. □ 

"Idempotence" Proof for Theorem 8. To see that the transformation T 8 defined in 
Theorem 8 is "idempotent" we first show 

T 8 o T 8 4 T 8 , (4.33) 

and then iterate. The precise two-step composition rule is 

T 8 (fT 2 , e 2 ) o T 8 (cri, ei) = T 8 ^cr 2 e>"i, e x a 2 + ^-J • (4.34) 

In this section we consider a two-step composition rule to find idempotence. In par- 
ticular, the most outstanding features of these steps are: 

• To see "idempotence", note that for fixed Co( r ) equation (j4.14j) has a solution 
space that is exactly two dimensional. 

• Since the first application of T 8 takes any specific solution and maps it into the 
full two-dimensional solution space, any subsequent application of T 8 can do 
no more than move one around inside this two dimensional solution space - 
physically this corresponds to a relabelling of parameters describing the perfect 
fluid metric you are dealing with, not the generation of new solutions. 

To be more explicit about this note that at step one 

B ^B 1 = B {a 1 + e l I }, (4.35) 



while at the second step 



B 1 ^B 2 = B 1 {a 2 + e 2 \ ). (4.36) 
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That is: 



B 2 = B < <ti + ei 



r dr 



x < 



cr 2 + e 2 



r dr 



5 (r) 2 [a 1 + e 1 / r dr/5 (r) 2 ] 2 



But this can be rewritten as 



(4.37) 



Therefore 



That is 



Bo 



e 2 / dZ 
° 2 7 / ~Z? 



£2 



1 1 

Z (Ti 



-- + 

ei 



e 2 1 
cr 2 H 

6i (Ti 



Zr 



B 1 = -^ + 
ei 

from which the composition law follows: 



a 2 



62 1 

ei cri 



T 8 (a 2 , e 2 ) o T 8 (o"i, ei) = T 8 ( a 2 cr ll e x a 2 + — 

0"! 



(4.38) 

(4.39) 
(4.40) 

(4.41) 
□ 



Note that this almost identical to the computation performed in the Schwarzschild 
coordinates because idempotence of the theorem is a result of "reduction of order" not 
something special to general relativity. 
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4.3.2 Two linking generating theorems: 

Having now found the first and second generating theorems, is now useful to find the 
connection between them. We start by taking a perfect fluid sphere solution {Co, -Bo} 
and applying theorem 7 onto it which gives us a new perfect fluid sphere {Ci, Bq}. The 
new Ci = £q . We now continue with applying theorem 8, again we get a new solution 
{Co ,Bi}, where B\ now depends on the B . All together we can consider this as a 
single process by introducing the following corollary: 



Corollary 3. If {Co, Bq] denotes a perfect fluid sphere, then for all a and e, the 
two-parameters geometry defined by 

ds 2 = --^ dt 2 + {dr 2 + r W} (4.42) 

Co(r) 2 B (r) 2 Z (rY 



Z (r) = U + e I -^-A. (4.43) 



is also a perfect fluid sphere, where Z$ is 
That is 

T 8 oT 7 : {Co, B } » {Co" 1 , B } » {Q 1 , B Z (B )} (4.44) 

Instead of starting with theorem 7 we could first apply theorem 8 on {Co, -Bo}. This 
gives us a new perfect fluid sphere {Co, -E>i}, where B\ = B Z (B ) is given by equa- 
tion (|4.32|) . We now continue with theorem 7 which leads to {Ci, -Bi}, where Ci = C -1 - 
Again, we can consider this as a single process by introducing the following corollary: 



Corollary 4. If {Co, -Bo} denotes a perfect fluid sphere, then for all a and e, the 
two-parameters geometry defined by 



1 ,.2 , Co(r 



,2 



ds 2 = --7-T77 dt 2 + - />"; ' - Idr 2 + r W} (4.45) 



Co(r) 2 B (r) 2 Z (r 

is also a perfect fluid sphere, where Z is 

rdr 
BJr 



Z (r) = {a + e \ - 777 -^ \ . (4.46) 



110 



That is 

T 7 oT 8 : {Co, B } » {Co, B Z (B )} ^ {Q\ B Z (B )} (4.47) 

There are general comments about the relationship between both transformation 
theorems. 

• These two transformation operations are independent of each other. 

• The transformation T7 defined in Theorem 7 is a "square root of unity" . 

• Theorem 8 is in general idempotent. That means when we apply theorem 8 more 
than once no further solution will be obtained. 

• Theorem 7 and theorem 8 are in general commute. 

In addition, these comments can be described in terms of the equations below. 

T 7 #T 8 ; T 7 oT 7 ^J; T 8 oT 8 ^T 8 ; T 7 o T 8 4 T 8 o T 7 . (4.48) 



4.4 Gaussian polar coordinates: 

Consider the metric 

ds 2 = -C(r) 2 dt 2 + dr 2 + R(r) 2 dn 2 (4.49) 

and then assume it represents a perfect fluid sphere. That is, G§g = G?? = Now 
Gqq = Gi^ is always fulfilled due to spherical symmetry. 
We calculate 

G„ (4 . 50) 

and 

Gee = Y( ( } 
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Seed metric 



Th7 



Th7 



new solution 



new solution 



Th7 



Th7 



new solution 



Th8 



ThB 



Figure 4.1: This structure shows that Theorem 7 and Theorem 8 are in general distinct. 
When we apply Theorem 7 to seed metric, we get a new solution, while applying Theorem 8 
to seed metric, we get the other new solution. Furthermore, the diagram shows that Theorem 
7 and Theorem 8 in general commute. 

We set Grf = G§§, which supplies us with an 2nd order ODE for C( r ); which reduces 
the freedom to choose the two functions to one: 

C"R 2 - £RR' + C[l - (R') 2 + R"R] = (4.52) 
We can rewrite this equation as 

e-C« +( { l -W + *« }=0 (4,3) 
4.4.1 The ninth theorem 

Theorem 9. Suppose {(o{r), Ro{ r )} represents a perfect fluid sphere. Define 
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Then for all a and e, the geometry defined by holding Ro(r) fixed and setting 

ds 2 = -Co(r) 2 A (r) 2 dt 2 + dr 2 + i? (r) W (4.55) 
is also a perfect fluid sphere. That is, the mapping 

T 9 (a, e) : {Co, R } ^ {Co A (Co, Ro), Ro} (4.56) 
takes perfect fluid spheres into perfect fluid spheres. 

Proof for Theorem 9. The proof of theorem 9 is again based on the technique of "re- 
duction in order". Assuming that {Co( r )? Ro( r )} solves again equation (|4.52|) . write 

Ci(r) = Co(r) A (r) . (4.57) 
and demand that {Ci( r )> Ro( r )} also solves equation (J4.52j) . We find 

(Co A )" - (Co A )'^ + (Co A ) { — (gZ±^jg | = (4.58) 

we can expand the above equation to 

(Co' A + 2C A' + Co AS) - (Co A + Co A' ) ^ + (Co A ) { 1 Z (g£+ gjg | = (4.59) 

we can also re-group this same equation as 

{Co - Cof + Co { 1 ' {R 'r2 +R " R }} A o + 2C A'o + Co K - Co A' | = 0. (4.60) 
A linear homogeneous 2nd order ODE for C now simplifies to 

(2C - Co f )A' + Co K = , (4.61) 

which is an ordinary homogeneous second-order differential equation, depending only 
on A' and A '. (So it can be viewed as a first-order homogeneous order differential 
equation in A' , which is solvable.) Separating the unknown variable to one side, 

T? = -2? + £ ■ (4-62) 
A o so tu$ 
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Re- write Aq/A' = dln(Ag)/dr, and integrate twice over both sides of equation (|4.62j) . 
to obtain 

Ao(rH f + e /w}' (4 ' 63) 

depending on the old solution {( (r), R (r)}, and two arbitrary integration constants 
a and e. □ 

"Idempotence" Proof for Theorem 9. To see that the transformation Tg defined in 
Theorem 3 is "idempotent" we first show 

T 9 o T 9 = T 9 , (4.64) 

and then iterate. The precise two-step composition rule is 

T 9 ((72, e 2 ) o Tg(eri, ei) = T 9 ^cr 2 e>"i, e x a 2 + • (4.65) 

In this section we consider a two-step composition rule to find idempotence. In par- 
ticular, the most outstanding features of these steps are: 

• To see "idempotence", note that for fixed Rq{t) equation ()4.53j) has a solution 
space that is exactly two dimensional. 

• Since the first application of T 9 takes any specific solution and maps it into the 
full two-dimensional solution space, any subsequent application of T 9 can do 
no more than move one around inside this two dimensional solution space - 
physically this corresponds to a relabelling of parameters describing the perfect 
fluid metric you are dealing with, not the generation of new solutions. 

To be more explicit about this note that at step one 

Co-Ci = Co + * >. (4-66) 



while at the second step 



Ci-C 2 = d {*2 + e 2 l^>- (4-67) 
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That is: 



C2 = Co < 01 + Ei 



Ro(r) dr 



x < 



o" 2 + e 2 



Ro(r) dr 



Co(r) 2 [ax + et / Ro(r) dr/Co(r) 2 P 



But this can be rewritten as 



C2 = Co Z { a 2 + 



£2 / dZp_ 



Co^o^ 2 -|/d(i- 



Co ^0 { 0"2 



£2 



1 1 



Therefore 



That is 



0"2 



£2 l_ 

ei o"i 



6 = -^ 

from which the composition law follows: 



e 2 1 
cr 2 H 



T 9 (a 2 , £2) T 9 (cri, ei) = T 9 ( a 2 (T\, e x a 2 + — 

0"! 



Again, note simularities to what happens in Schwarzschild coordinates. 
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Seed metric: 
Minkowski 



Th9 




■Th9 



Figure 4.2: This structure shows that Theorem 9 is idempotent. When we apply Theorem 
9 to Minkowski as a seed metric, we get the K-0 III as a new solution. 

4.5 Buchdahl coordinates: 

Consider the metric 

ds 2 = -C(r)" 1 dt 2 + C(r)dr 2 + ({r)R{r) 2 dtt 2 (4.73) 

and then assume it represents a perfect fluid sphere. That is, G§§ = G?? = Gii. Note 
Gqq = G^ is always fulfilled due to spherical symmetry. We calculate 

r l f -4C 2 (^) 2 + 4 c 2 + (CO 2 R 2 \ (A7A , 

^-\{ H2R ";p 2R ) <«»> 

We set Gff = G§§, which supplies us with an 2nd order ODE for C( r ) ; which reduces 
the freedom to choose the two functions to one: 

'Cy_ 2[1 - (R'f + RR") ^ (4 76) 

This coordinate system is a sort of cross between isothermal and proper radius. 

Theorem 10. If {((r), R(r)} describes a perfect fluid then so does {((r)" 1 , R(r)}. 
This is the Buchdahl transformation in yet another disguise. 
The geometry defined by holding Ro{r) fixed and setting 

ds 2 = -Co(r) dt 2 + -l^dr 2 + ^fldQ 2 (4.77) 
Co(r) Co(r) 
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is also a perfect fluid sphere. That is, the mapping 

T 10 :{( 0l R }^{Q\R } 
takes perfect fluid spheres into perfect fluid spheres. 



(4.78) 



4.6 Synge isothermal coordinates: 

Consider 

ds 2 = -C(r)~ 2 {dt 2 - dr 2 } + {((ry 2 R(r) 2 dn 2 } (4.79) 

and then assume it represents a perfect fluid sphere. That is, G^ = Gp? = Gii. While 
Gqq = Gii is always fulfilled due to spherical symmetry. We calculate 

Gff = _ 3R *(cr-u«R R CHRre-e (480) 

and 

3R(C) 2 - 2(R'C - 2R ((" + CR" (Asn . 

G§§ = ^ I 4 - 81 ) 

We set Gff = Ggg, which supplies us with an 2nd order ODE for C( r ) ; which reduces 
the freedom to choose the two functions in equation 

,-,£-c{i^£}-a 

This is very similar equation to Gaussian polar. 

Theorem 11. Suppose {( (r) , R (r)} represents a perfect fluid sphere. Define 



R (r)dr 

c^ 7 



Mr)-{c + .r-^\. (4.83) 



Then for all a and e, the geometry defined by holding Ro(r) fixed and setting 

ds 2 = -—J—— {dt 2 - dr 2 } + -J^f—dQ 2 (4.84) 
Co(r) 2 A)(r) 2 1 J Co(r) 2 A (r) 2 v ; 
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is also a perfect fluid sphere. That is, the mapping 

Tu(a, e) : {Co, Ro} » {Co 4)(Co, Ro), Ro} (4.85) 
takes perfect fluid spheres into perfect fluid spheres. 

Proof for Theorem 11. The proof of theorem 11 is based on the technique of "reduction 
in order". Assuming that {( (r) , Ro(r)} solves equation ([4.82)1 . write 

Ci(r) = Co(r) Mr) ■ (4-86) 
and demand that {(i(r) , Ro(r)} also solves equation (|4.82j) . We find 

(Co A )" - (Co A )'| + (Co A ) { 1 ~ {R '^ RR " } = (4.87) 
we can expand the above equation to 

(Co' A) + 2C A' + Co K) ~ (Co A + Co A' ) ^ + (Co A ) { 1 ~ ( ^ 2 + — } = (4.88) 
we can also re-group this same equation as 

{C - Cof + Co { l ~ {R 2 R t RR " }} A « + 2 Co < + Co K ~ Co = 0. (4.89) 

A linear homogeneous 2nd order ODE for C now simplifies to 

(2C - Co f K + Co < = , (4.90) 

which is an ordinary homogeneous second-order differential equation, depending only 
on A' and A' Q '. (So it can be viewed as a first-order homogeneous order differential 
equation in A' Q , which is solvable.) Separating the unknown variable to one side, 

A" P R' 

^ = -2f + f- (4-91) 
^o ^o Mo 

Re- write Aq/A' q = dln(y4' )/dr, and integrate twice over both sides of equation ()4.91|) . 
to obtain 

depending on the old solution {(o{r) , Ro(r)} , and two arbitrary integration constants 
a and e. □ 
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4.7 Exponential coordinates: 



You can always choose to effectively use g u itself as one of the coordinates, by picking 
the functional form of gu to be any arbitrarily specified function f(z). For instance, 
choose 

g tt = -exp(-2z) (4.93) 

then a useful choice is 

f dz 2 1 

ds 2 = - exp(-2z) dt 2 + exp(+2z) <^ — — + R(zf dQ 2 } (4.94) 



,B(z) 

and then assume it represents a perfect fluid sphere. That is, = G^z = G^j. Note 
Gqq = G^ is always fulfilled due to spherical symmetry. Calculate 

(BR 2 -(R>) 2 B + l)eM-2z) 
^zz - ™ (4.yt>J 

^ r 1 (B'R' + 2R"B + 2BR)exp(-2z) 

We set Gzz = G§§, which supplies us with an 1st order ODE for B, which reduces 
the freedom to choose the two functions in the metric to one for which the isotropy 
constraint becomes 

B'[RR'] + B[AR 2 - 2(R') 2 + 2RR"} + 2 = (4.97) 
Theorem 12. Suppose {Bq(z), Rq(z)} represents a perfect fluid sphere. Define 

t ^(1f-{- 4 /1W- (498) 

Then for all X, the geometry defined by holding Rq(z) fixed and setting 

f dz 2 1 
ds 2 = - ex V (-2z)dt 2 + exp(+2^) { — — — — — + i^) W \ (4.99) 

is also a perfect fluid sphere. That is, the mapping 

Ti 2 {(7, e) : {B , Ro} ^ {B T (i2o), Ro} (4.100) 
takes perfect fluid spheres into perfect fluid spheres. 
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Proof for Theorem 12. Assume that {B (z), Rq(z)} is a solution for equation ([4.97)1 . 

We want to find under what conditions does {Bi(z), Rq(z)} also satisfy equation 
(|4.97|) ? Without loss of generality, we write 

B 1 (z) = B (z) + XT (z). (4.101) 

Equation (|4.97|) can now be used to determine T (z). Substitute Bi(z) in Equation 

flOZD 

(B + A To)'^^] + (Bo + A T )[4i2g _ 2 (R' ) 2 + 2RqIQ + 2 = (4.102) 

we can also re-group this same equation as 

[RoR' }B' + [4Rg - 2(R' ) 2 + 2R R' r \B + 2 

+ [R R' ] T' + [4^ - 2(it" ) 2 + 2R R' r \ T = (4.103) 

That ordinary inhomogeneous first-order differential equation in B now simplifies to 

[R R' } T' + [ARl - 2(R' ) 2 + 2R R' r \ T = , (4.104) 

which is an ordinary homogeneous first-order differential equation in T . 
Now we can calculate Tq by following, 



[R R' ] T' = - [ARl - 2(R' ) 2 + 2R R' r ] T . (4.105) 



Rearrage it into 



JJ, = - [AR 2 - 2(R' ) 2 + 2Rp]Q 

= -4§ + 2§-2§ (4106) 

JXq rto rL Q 

A straightforward calculation, including an integration by parts, leads to 

t °m = (S^{- 4 /tH- (4107) 

□ 
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4.8 Discussion 



We have developed several transformation theorems that map perfect fluid spheres 
to perfect fluid sphere using other coordinates. We derived the seventh and eighth 
theorems using Isotropic coordinates. The transformation in Theorem 7 is a "square 
root of unity". Indeed, Theorem 8 is alomost identical to the calculation performed 
in the Schwarzschild coordinates. In addition, the idempotence of this theorem is a 
result of "reduction of order" , which is not special to general relativity. We are able to 
classify the connection between both transformation theorems. The two transforma- 
tion operations are independent of each other. Furthermore, both theorems in general 
commute. Figure l4~Tl shows that after applying Theorem 7 to the seed metric, we get a 
new solution while applying Theorem 8 to seed metric, we get the other new solution. 
By the idempotence of Theorem 8 n applications of Ts to the seed metric still results 
in the same non-seed metric. Furthermore, the diagram shows that Theorem 7 and 
Theorem 8 in general commute. 

In Gaussian polar coordinates, we established another transformation theorem, 
which has many similarities to the results in Schwarzschild coordinates. Figure 14.21 
shows that Theorem 9 is idempotent. When we apply Theorem 9 to Minkowski as a 
seed metric, we get the K-0 III as a new solution. 

Furthermore, we have also established the transformation theorem using Buchdahl 
coordinates. Theorem 10 is the Buchdahl transformation in yet another disguise. 

In Synge isothermal coordinates, we presented Theorem 11 is a transformation the- 
orem. However, this transformation is almost identical to the computation performed 
in the Schwarzschild coordinates. 

The last solution generating theorem we presented using exponential coordinates. 
Indeed, this Theorem 12 also map perfect fluid spheres to perfect fluid spheres. 
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Chapter 5 

TOV equation - related theorems 



5.1 Introduction 

The Tolman- Oppenheimer- Volkoff (TOV) equations describe the hydrostatic equilib- 
rium of perfect fluid stars without rotation in general relativity. We can solve the 
TOV equations by constructing the equation of state for the entire region of nonzero 
density starting from the higher density at center to the surface density. In this chap- 
ter, we can rephrase Theorems 1, 2, 3 and 4 directly in terms of the TOV equation. 
In addition, we also can rephrase all these theorems, which apply in Schwarzschild 
coordinates directly in terms of the pressure profile and density profile. 

5.2 TOV equation 
5.2.1 Interior equation 

In general relativity [21] , continuous matter distributions and fields again are described 
by a stress-energy tensor T a b. The stress tensor of a perfect fluid is given by 

Tab a b , „ / ab , a b\ /r i\ 

= p u u + p [g + u u ) (5.1; 

where p(r) is density of mass-energy distribution in rest-frame of fluid, u is 4-velocity 
of fluid, and p is isotropic pressure in rest-frame of fluids. 
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A* = U = (5-6) 



We know that 

G^SttG/vT^ (5.2) 

We set Gtv = 1 for convenience, and note that: 

V a G ab = (5.3) 

We know that as contracted Bianchi identity. 
The equation of motion for the matter is 

V b T ab = (5.4) 

We define 

dr a 

u« = ^ (5.5) 

This is the 4- velocity of a fluid element. We also define 

d 2 x a du a 

d^~ = ~dV 

This is the 4-accerelation. 

In the rest frame of the fluid element u a = (c, ) and, A a = (0, ~a). 
It follows that u a A a = 0. 

It is an invariant quantity. In other words, the 4-acceleration of a particle is always 
orthogonal to its 4-velocity. Writing 

Ts = {uV) (5 ' 7) 

We see 

= -^(-1) = ^-(g ab u a u b ) =2u a ^-u a = 2u a (uV)u a = 2u a A a = (5.8) 
ds ds ds 

Now consider VftT afe ; we can write this as 

= (V b p)u a u b + (W b p) (g ab + u a u b ) 

Hp + p)(v b u V + ( P + p)u a (y b u b ) 

= [(V b p)u b + (p + p)(V-u)]u a 

+ [Go + p)A a + V b p(g ab + u a u b )] 
= (5.9) 
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Writing out equation (|5.4j) in terms of p, p, and u b , and projecting the resulting 
equation parallel and perpendicular to u a , we can spilt up the equation into two 
equations: 



[(V h p)u b + (p + p)(V-u)]u a = 0, 
[(p + p)A a + V bP (g ab + u a u h )] = 0. 



(5.10) 
(5.11) 



with u a u r , 



-1, and A a u a = 0, because they are orthogonal to each other. 



if 



[V b ( itlO iP + ( itlO + p ) ( V U ) )| Li 5 



Figure 5.1: The diagram shows that when we write out equation (|5.4|) in terms of p,p, and 
u b , and projecting the resulting equation parallel and perpendicular to u a , we can spilt up 
the equation in to two equations. 

Now use (Vbp) g ab u a = (Vbp)u b , and (V&p) u a u b u a = —(Vbp)u b . This means 
(VbP) g ab u a and (VbP) u a u b u a have the same magnitude but opposite direction. We 
conclude 

V • (pu) +pV -u = (5.12) 

and 

Aa = (g ab + u a u b ) V b p 
p + p 
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5.3 A perfect fluid sphere in Schwarzschild coordi- 
nates 



In the Schwarzschild cooridinates |24| . the metric of an arbitrary static, spherically 
symmetric spacetime can be put into the simple form 

ds 2 = -f(r) dt 2 + h(r) dr 2 + r 2 (d6 2 + sin 2 6 d</> 2 ) (5.14) 

when Gn = 1, then G a b = 87rT a b. 
The Einstein equations yield, 

87rp = Goo = Roo + \ W + Ri 1 + R2 2 + Rs 3 ) 

= (rh 2 )- 1 ti + r- 2 (l-h- 1 ), (5.15) 



lirp = Gn = 

-1 t i 



{Ro° + Ri 1 + R2 2 + -R3" 



(5.16) 



fcrp = G» = i(//i)- 1/2 ^:[(/^- 1/2 /1 



(5.17) 



We assume /i is 



2m(r) 



n -1 



(5.18) 



and if we write 



/ = exp(20), 



(5.19) 



the equation ()5.16|) becomes 



d0 [m(r) + 4 7rp(r) r 3 
dr r 2 [1 — 2m(r)/r] 



(5.20) 
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Proof for hydrostatic equilibrium . We substitute equations (|5.18|) and (|5.19|) into equa- 
tion ()5.16|) we get, 

8tt P = (rfhY 1 -r" 2 (l - h' 1 ) 

f i (i-O 



/ rh r 2 

2m(r) 2m(r) 



d0 1 r r 



2 



d0 r — 2m(r) 2m(r) 

,y2 i^«3 



d(^> r — 2m(r) 2m(r) 
2 — ^ = 8vrp 1 



m(r) 



47rp + 



dr r — 2m(r) 
t 

[m(r) + 4 irp r 3 ] 



(5.21; 



r 2 [1 — 2m(r)/r] 

as previously illustrated. □ 

To see what happens to the pressure and density it will be useful to consider the 
TOV equation. 

dp ~[p( r ) + p( r )} \m{r) + 4irp(r) r 3 } 



dr r 2 [1 — 2m(r)/r] 



(5.22) 



Proof. Consider the stress energy of a perfect fluid is in the equation ()5.1|) and it 
satisfies the equation ()5.4|) which yield, 

(p + p)u a W a u b + (g ab + u a u h )V a p = (5.23) 

Now consider, 

(e 1 ) b V a T ab = 0, (5.24) 
where, (ei)^ = /i 1 / 2 (dr) b = h^ 1 / 2 {d/dr) b 

We can reduce the equation ()5.24j) into the equation ()5.25|) because (e\) b is pointing 
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in the r-direction, which is perpendicular to u a . 

(e 1 ) b [(p + p)u a V a u b +(6 b + u aU b )V a p} = (5.25) 
h- 1 / 2 (0,l,0,0) b [(p + p)u a V a u b +(5 b a + u a u b ) V>] = 

^(p+pjKvA + r 1 /^ = o 

dr 

Now we consider, 

(u a V a u b ) = u a [d a u b + T b ca u c ] = u a T b ca u c 

= j = — -g bc [—gu,e + 9te,t + 9et,t] 

= VV=^ (5.26) 
dr 

Now we substitute equation ()5.20|) into equation (j5.25|) we get, 

h^( P + P )^ + h^^ = 
dr dr 

dp = r p + p j<ty 

dr dr 

dp = [ m (r) + 4^(r)r»| 

dr r z [l — 2m{r)/r\ 

Equation (|5.27|) is known as the Tolman-Oppenheimer-Volkoff equation of hydrostatic 
equilibrium. □ 

This equation applies only for a perfect fluid (p r = pt) in Schwarzschild coor- 
dinates. There are generalizations appropriate to other coordinate systems. In the 
non-relativistic Newtonian approximation this simplifies to (p <C p, 2m(r)/r <C 1): 

dp = _|pjrW^ / g 2g N 

dr r 2 

which is something that is very easy to derive in Newton's theory of gravity. 

If we assume p(r) and hence m(r) are given, then the TOV is a nonlinear first-order 
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ODE for p{r) called a Riccatti equation. Given p(0) a unique solution p(r) necessarily 
exists, and can be used to reconstruct the function £(r) via: 

Proof. Firstly, start with some static spherically symmetric geometry in Schwarzschild 
(curvature) coordinates 

Hr 2 

ds 2 = -C(r) 2 dt 2 + —— + r 2 dQ 2 (5.30) 
B(r) 

and then assume it represents a perfect fluid sphere. 
When 

2BC'r -C + BC 



Inp = G r 
iirp 



r2 C 

2B C 1 B 

2B C IB 
— — = 8np + — - — 

C 4npr (1 - B) 



C B 2Br 

when B = 1 — 2m(r)/r, then we get 

£' 47rpr 1 — (1 — 2m(r)/r) 



(5.31; 



£ (1 — 2m{r)/r) 2 (1 — 2m(r)/r)r 
£' m(r) + 47rp r 3 

C r 2 (1 — 2m(r)/r) 
dln(C) m(r) + 47rp r 3 

dr r 2 (1 — 2m(r)/r) 

^ J \r 2 (l-2m(r)/r) / 

exp(ln(C)) = exp [J r2 {1 _ 2m(r)/r] dr 

= f /" Hr)+47rpr 3 ] d J 

\ J r 2 [1 — 2m(r)/r] J 

as required. □ 



129 



In this way, starting from p(r) you can construct first the function m(r) = 4tt J p(r)r 
then (via the TOV) the pressure profile p(r), and finally via the above the gravitational 
redshift C( r )- 

Suppose in contrast that the pressure profile p(r) is given, then the TOV equation 
can be rearranged into a nonlinear first order PDE for m(r): 

1 dm(r) r 2 [1 — 2m(r)/r] dp 

Pv) ~ r„^„\ , a- „^„\ „ai ~JZ. (5.33) 



Airr 2 dr [m(r) + Air p(r) r 3 ] dr 

Proof. We substitute equations (|5.18jl and (|5.19jl into equation (|5.15j) to get 

dm(r) 



dr 

and rearrange in term of p(r) so, 



4irp(ry, (5.34) 



1 dm(r) , . 



rearrange the equation ()5.22|) into, 



r 2 [1 -2m(r)/r] dp 

p(r) = -p( r ) - rV^TX" ^ 5 - 36 ) 

[m[r) + 47r p(r) r 6 \ dr 



now substitute equation ()5.35|) into equation (|5.36jl . we get 



1 dm(r) _ _ r 2 [1 - 2m(r)/r] dp ^ ^ 



47rr 2 dr [m(r) + 47r p(r) r 3 ] dr 



□ 



This is an Abel equation (2nd type, class A). There is no simple general solution. 
However given p(r) and the initial data m(0) = this can in principle be solved to 
determine m(r). Once this is done £(r) is again reconstructed via 



= | f [m(r)+ 47^(^3] | 
sv ' F \J r 2 [l - 2m(r)/r] J V ; 
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5.4 Solution generating theorems 



5.4.1 Four theorems 

The first four theorems we presented in chapter 3 were first explicitly published in 
0. We now rephrase the theorems in slightly different formalism, and demonstrate 
an independent way of proving them. Indeed, we can re-write these theorems directly 
in terms of pressure and density profiles p(r) and p(r). 

Firstly, start with some static spherically symmetric geometry in Schwarzschild 
(curvature) coordinates 

ds 2 = -Co(r) 2 dt 2 + + r 2 dtt 2 (5.39) 



1 



2m (r ) 



and then assume it represents a perfect fluid sphere. Our original four theorems can 
be written in the form 

Theorem 1. Suppose {(o{ r ) , m o{ r )} represents a perfect fluid sphere. Define 

CoW V , f fa(r) Co(r)-r^(r) 



mi(r)=( — "2 ' ^ ex P 2 / 7TT 7TTT ? dr ' ( 5 - 40 ) 
\(o{r)+rQ{r)J { J Co (r) Co (r) + r Q(r) J 

Then for all X, the geometry defined by holding Co( r ) fixed and setting 

ds 2 = -Co(r) 2 dt 2 + -. . ™ + r 2 dQ 2 (5.41) 



^ 2[m (r) + A mi(r) 
is also a perfect fluid sphere. 

Theorem 2. Let {Co( r ) ; m o( r )} describe a perfect fluid sphere. For all a and e, the 
geometry defined by holding m (r) fixed and setting 

ds 2 = -Co(r) 2 I a + e / - r df I d« 2 + + r 2 d« 2 (5.42) 

is also a perfect fluid sphere. 
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Theorem 3. Let {Co( r )> m o( r )} describe a perfect fluid sphere. For for all a, e, and 
X, the three-parameters geometry defined by 



ds' 



<o(r) 



r dr 



a + e 



j(r)a v / 1 _ 2K(r) + Am 1 (r)l 



dt 2 



dr 2 



2[mo(r) + A mi(r) 



+ r 2 dfi 2 



(5.43) 



z's a/so a perfect fluid sphere, where mi(r) z's defined as above. 



Theorem 4. Lei {Co( r ) ; m o(^)} describe a perfect fluid sphere. For all a, e, and A, 
the three-parameters geometry defined by 

2 



ds 2 = -Co(r) 



r dr 



a + e 



Co(r) 



1 - 



2m (r) 



dt 2 + 



dr s 



^ 2[m (r) + A mi(r) 



is also a perfect fluid sphere, where nil (r) is defined as 



mAr) 



X(r) + r C'(r) 
depending on ( = (od, where 

Ci(r) = 



r exp < 2 



C(r) C(0 + rC'(0 



dr 



(5.44) 



(5.45) 



r dr 



Co(r) 



1 - 



2mo(r 



(5.46) 



5.4.2 Four theorems in terms of pressure and density 

We can re- write these theorems in terms of pressure and density profiles p(r) and p(r). 

Consider Theorem 1. By definition ( is unaffected under the transformation. But 
we have a formula for ((r) in terms of p(r) and m(r). Let m2(r) = mo(r) + Ami(r) 
and P2{r) be the mass profile and pressure profile after applying theorem 1. Then 

[mo(r) + 47rpo(r) r 3 ] [m2(r) + 4Trp2(r) r 3 ] 



r 2 [1 — 2m (r)/r] r 2 [1 — 2m2(r)/r] 
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(5.47) 



which we can rearrange to yield 



m 2 (r) + Airp 2 (r) r = [mo(r) + Airpo(r) r ] 



3 [1 - 2m 2 (r)/r] 



[1 — 2m (r)/r] 



(5.45 



Furthermore, we can rearrange it into 



A7rp 2 {r) r 3 = -m 2 {r) + [m (r) + 4vrp (r) r 3 ] E j^^Tl ( 5 - 49 ) 

[1 — 2mo{r)/r\ 

Now we substitute m2(r) = mo(r) + Ami(r) into equation ()5.49j) . 

/ \ s / n / \ r / n / \ 3i [1 _ 2(m (r) + Xm 1 (r))/r] 

Airp 2 (r) r 6 = -m (r) - \m\{r) + [m (r) + 4:irp (r) r 3 ] — 



[1 — 2m (r)/r] 



(5.50) 



After a bit of algebra: 

3 -Ami(r) + 47rp (^) — 87rp (?") f 2 ^o( r ) ~ 8A7rp (?") f 2 mi(r) 



4irp 2 (r) r £ 



[1 — 2m (r)/r] 
-Ami(r)(l — 87rpo(^) ^ 2 ) 47rpo(f) ^ 3 [1 — 2mo(r)/r] 



= 4 7rp (^)^ 
rearrange it in terms of p 2 (r), 



[1 — 2mo(r)/r] [1 — 2mo(r)/r] 

3 Ami(r)(l + 87rp (^) ^ 2 ) 



[1 — 2mo(r)/r] 



(5.51) 



= PoW - Ami(r) 



47rr 3 



^7rj9 (r) r 2 



1 — 2m (r)/r 



(5.52) 



This allows us to rewrite Theorem 1 as follows: 



Theorem 13 (Theorem IB:). Suppose we are given a pair of functions Po(r) and 
Po(r) that satisfy the TOV. 

We use these functions to construct mo(r) , Co( r ); and m\{f) as follows: 



m (r) =4:% Po(r) r 2 dr 
Jo 



(5.53) 



Co(r) = exp 



[wio(r) + 4np (r) r 3 ] 
r 2 [1 — 2mo(r)/r] 
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(5.54) 



mi(r) 



Co(r) 



r exp < 2 



g(r) Co(r) - r Q(r) 
Co(r) Co(r)+r^(r) 

1 + 87rp (^)^ 2 



.CoO) + r Co( r ) 
T/ien for all X, the pair of functions 

P2(r) = Po(r) - Ami(r) 

P2 Po ) + 4.^2 ^ r 
a/so satisfies the TOV and so defines a perfect fluid sphere. 



dr 



47rr 3 



2mo(r)/r 



X dmi(r) 



(5.55) 



(5.56) 



(5.57) 



Similar logic can be applied to theorem 2. Note that in that case you are considering 
the transformation 

m o( r ) — * ^( r ) = r ^o( r ) (5.58) 



Co(r) - Ca(r) = Co(r) < a + 6 / == )> 



2mo(r| 
— 



So that mo(f) remains fixed. Now from the Gf? component of the Einstein tensor: 
consider the equation below 



£' m(r) + 47rp(r) r 3 
C r 2 (1 — 2m{r)/r) 

and rearrange this equation in term of p(r) 

— r 2 (1 — 2m{r)/r) = m(r) + Anp{r) r 3 



when 



(5.60) 



47rp(r) r 3 = — m(r) + r 2 (1 — 2m(r)/r) — 

p(r) = + — (1 - 2m(r)/r) ^ (5.61) 



m(r) _ m(r)/3 _ 1 
4vrr 3 " (4vr/3) r 3 ~ 3 P{T) 
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Now we have p(r) equal to 

p(r) = 

Then in particular 



so that 



1 1 C'(r) 

3 P{r) + ^Tr {1 - 2m(r)/r) W) 



(5.63) 



l l C'( r ) 

P2O) = -- po{r) + - — (1 - 2mo(r)/r) 

3 47rr C2(, r J 



(5.64) 



J*(r) = -1 po(r) + JL (l-^gM 



Co(r) < (7 + e 



r dr 



Co(r) 2 \A 



2 / 1 2m (r) 



Co(r) {<y + e 



r dr 



a 1 \i 2m (r) 



(5.65) 



That is 



1 / n 1 A 2m (r)\ CoW 



1 2mn(rl 



47rr 



1 - 



/• i \2 A 2m (r) 
Co(r) 2 \/l fM 



(5.66) 



cr -h e 



r dr 



1 - 1 C'(r) 

we know p(r) = — 3 p(r) + (1 — 2m(r)/r) ^Jyj > so that 



f J_^ 1 _2n^rl Co(r) _ 2 < ) a + e 



r dr 



Co(r) Vl - ^ 



(5.67) 



This leads us to: 
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Theorem 14 (Theorem 2B:). Suppose we are given a pair of functions po(r) and 
Po{r) that satisfy the TOV. 

Use these functions to construct m (r) and Co( r ) as follows: 

m o( r ) = 47T / po(f) f 2 df (5.68) 
Jo 

T/ien /or all a and e, the pair of functions 

p 2 (r) = p (r) (5.70) 



-i 



e / 1 2mo(r) A , N _ 2 / r dr 



r 



a/so satisfies the TOV and so defines a perfect fluid sphere. 

Having now found these theorems IB and 2B, is now useful to define a new theorem 
by composing them. All together we can consider this as a single process by introducing 
the following theorem: 

Theorem 15 (Theorem 3B:). Suppose we are given a pair of functions po(r) and 
Po(r) that satisfy the TOV. 

Use these functions to construct m (r) , Co( r )> and m\{r) as follows: 

m (r) = Air Po(^) ^ 2 df (5-72) 
Jo 

, i x ff [m (r) +4vrp (r) r 3 ] \ 

C ° (r) = ^ \J [l-2mo(r)/r] ^ J ^ 
Mr) = (,,?° {r l,X r 3 exp / 2 / ^ ^- r ^ ( 1 dr) (5.74) 



Co(r) + r C (r)J " [ 7 Co(r) Co(r) + r Co(r) 
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Then for all a, e, and X, the pair of functions 



p 2 (r) = p (r) - Xmi(r) 



Airr 3 



]np (r) r 2 



1 — 2nio(r)/r 



+h V 1 - ^ ^ 



r dr 



a + e 



Co(0 



2 L _ 2m 2 (r) 



(5.75) 



P2 = Po(r) + 



A dmi(r) 



4iir 2 dr ^ 
where rri2{r) = mo(r) + Ami(r), a/so satisfy the TOV and so define a perfect fluid 

sphere. 

Theorem 16 (Theorem 4B:). Suppose we are given a pair of functions po(r) and 
Po(r) that satisfy the TOV. 

Use these functions to construct mo(r), Co( r ); an d mi(r) as follows: 



(5.77) 
(5.78) 
(5.79) 



m o( r ) = 47r / Po( r ) r2 dr 
Jo 

[roo(r) + Airp (r) r 3 ] 



mAr) 



Coir) = exp jy 

C(0 



dr 



X(r) +r C'(r) 
depending on ( = CoCi? where 



r exp < 2 



Ci(0 = < cr + e 



r 2 [1 — 2mo(r)/r] 

C'(r) {(r)-rC(r) 
C(r) C(r)+r(>{r) 

r dr 



dr 



Then the pair of functions 

P2{r) = Po{r) - Xmi(r) 



1 + 87rpo( r ) r2 



(5.80) 



47rr 3 



1 — 2m (r)/r 



P2 = Po(r) + 



X dmi(r) 



4 7rr 2 ^ r 

also satisfy the TOV and so define a perfect fluid sphere. 



(5.81) 
(5.82) 
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Some seed geometries and their descendants in term of pressure and density via Theorem IB 



Seed 

Minkowski 

exterior Schwarzschild 
de Sitter 

Tolmaii V (A = 0) 
SI 

M-W III 

Heint Ha (C = 0) 



Mr) 




3 

WrrW 



(»+l) 2 

8?' 2 7r(n 2 -2) 
1 

2-rrr' 2 

(q-6r) (a+2r) 
56 7T r 2 a 2 

9a(-l+ar 2 ) 
167r(l+ar' 2 ) 2 



PThi(r-) 

3500 A (r- 



Po(r) " 

Po(r) - 
Po{r) + 

Po(r) ~ 



■of 

A(,B 2 -3r 2 ) 
47r(-i? 2 +2r 2 ) 

/ (2r» + „ 2 ) ' 

ArV 



(« + 2) 



5A 



4 7r(n+2) 2 



(2w+3) 



36 7rr< 4 /3) 

^r^Cl) (r-a)(2a+3r) 



9(a+r)(-4r+3a) ( 3> 
A(7ar 2 + 1) 



12\/3 7r(l+aT- 2 ) 2 V-4ar 2 -l 



Po(0 





3 

KivK 2 



+ 

8r 2 7r(-2+n 2 ) 



3(a 2 +4r 2 ) 
56 7rr 2 a 2 

3a (3+ar 2 ) 
167r(l+ar' 2 ) 2 



PThi (r) 



3A 
in 



3500 A (i--a)(3r-5a) 



A (7fl 2 r 2 -3-R 4 -6r 4 ) 



ArV 



PoW 



(™ + 2 ) 



5A 



4 7r(n+2) ;i 



(3+2n)(2-n) 



108 7rr( 4 / 3 ) 



(-4r+3a) ( 3 ) 



PoW - 



A(l+3ar 2 ) 



4V37r(l+ar 2 ) 2 (-4ar 2 -l) ( 2> 



Table 5.1: Seed solutions and their generalizations derived via theorems IB in terms of 
pressure and density profile. 



Some seed geometries and their descendants in term of pressure and density via Theorem 2B 



Seed 



Prh2(r) 



PTh2(r) 



Minkowski 

exterior Schwarzschild 
de Sitter 

Tolman V (A = 0) 
SI 

M-W III 
Heint Ila (C=0) 



47T [a+l 



2 7r(r 2 +5ra-30 a 2 + 2 " ^ r< f -^+15 A /r(r-2a) a 2 ln(r-a+y/r(r-2a))) 

Po(r) + 



Po{r) + 

Po{r) + 
Po(r) + 

Po(r) + 



2Trr 2 [2B 2 nr( 1 + 



n) /_L 
V 2-n 



4 7rr 4 a- 



2\/7 nAraa(r~a) . ~ 

— ~T ' 7 

Ya 2 — r 2 



g7r (l + ar 2 )( 7 / 2 ) / 2e (5+2ar 2 ) 

\JT-2ar 2 \ 27a ( 1 + ar 2 ) ( 3 / 2 ) \Z4-2ar 2 






PoO") 

Po(r) 


PoM 
Po(r) 



Table 5.2: Seed solutions and their generalizations derived via theorems 2B in terms of 
pressure and density profile. 
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Some seed geometries and their descendants in term of pressure and density via Theorem 3B 



Seed 

Minkowski 



PTh3(V) 

Po(r) 



\/l-2Ar 



exterior Schwarzschild 



[integral 



de Sitter 



M _ A(i? 2 -3r 2 ) 
POV ) 47r(-_R 2 +2r 2 ) 



, r 2 2Ar 2 (-H 2 +r 2 ) 
(-K 2 +2r 2 ) 



47r(fl 2 -r 2 ) 2 



' y fl 2 (-H2 +2 i-2) 



dr 



<l + A)(-2„) 

Ar ( 2 +"> (3+2n) 
4 7i-(2+ra) 2 



Tolman V (A = 0) 



Po(r) 



1 2Ar ( 2 + ") 

T + ,„ , ..,2 



(2-n 2 ) (2 + n) 2 



47rB 2 r (l + ra )^ 



r (-l-2r.) 



15 2 



> dr 



2(2-2ti2) 

1 , 2ir ( 2 + "> 
(2 ti 2 ) (2 + n) 2 / 



SI 



Po(r) 



12r( 4 / 3 > • 



27rr 4 



V9-2Ar( 2 / 3 > 



L.-L I 1 | A / 1 



P°V'J 9 V5/ (-4r+3a)( 4 / 3 > 



(473) T 9 I 2r (2/3) 



3V9-2 Ar( 2 / 3 ) 



arctanh 



V9-2A r ( 2 / 3 > 



M-W III 



(q 2 -r 2 ) 2Ar(V 3 ) (r-a)5( 2 / 3 ) 
H 45(-4r + 3a)< 4 / 3 > 



27rAr(r — a) 



cr+e 



., , / c 2 , 2 , , ;.V- 1 - -0)5(2/3) 
7a 2 + 45(-4r+3a)( 4 / 3 ) 



dr 



A(7ar 2 + 1)~ 



Heint Ila (C=0) 



Po(r) 



^/^(W 2 )- 



3\/-12ar 2 -3 



47r(l+ar 2 ) 3 



: dr 



4Ar 2 
3\/-12ar2_3 



Table 5.3: Seed solutions and their generalizations derived via theorems 3B in terms of 
pressure and density profile. The notation "[integral]" denotes a metric so complicated that 
explicitly writing out the relevant integral is so tedious that it does not seem worthwhile. 
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Some seed geometries and their descendants in term of pressure and density via Theorem 4B 



Seed 




PTh4(r) 


Minkowski 


PThlO") 


A(5er^+6cr) 


exterior Schwarzschild 


PThl( r ) 


[integral] 


de Sitter 


PThl( r ) 


[integral] 


Tolman V (A = 0) 


PThl( r ) 


[integral] 


SI 


Ptm(0 


\(10 err* -3 e) 
47r(6<rr 2 -e)( 5 / 3 ) 


M W III 


PThlO) 


[integral] 


Heint Ha (C=0) 


PThl( r ) 


[integral] 



Table 5.4: Seed solutions and their generalizations derived via theorems 4B in terms of 
pressure and density profile. The notation "[integral]" denotes a metric so complicated that 
explicitly writing out the relevant integral is so tedious that it does not seem worthwhile. 

5.4.3 Discussion 

In this chapter we have rephrased the theorems 1, 2, 3, and 4 of chapter 3 directly 
in terms of the TOV equation, and consequently in terms of the pressure profile and 
density profile. We have done this by rephrasing theorem 1 in a slightly different 
formalism, and developing an independent proof that is easy to rewrite in terms of 
the pressure and density profiles — pir) and p(r) — resulting in the theorem we call 
theorem IB. Similar logic, involving the transformation T2, can be applied to theorem 
2 — now leading to theorem 2B. With these new theorems IB and 2B in hand, we 
can compose them in various ways to obtain theorems 3B and 4B. 

We present several tables that exhibit seed solutions and the geometries generated 
from them by combining all the various theorems. Table (|5.1jl shows some seed ge- 
ometries and their descendants as obtained via theorem IB, and similarly table (|5.2jl 
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shows some seed geometries and their descendants obtained via theorem 2B. Finally 
table (|5.3|) and (|5.4|) deal with theorems 3B and 4B. We emphasise that theorem 2B 
does not affect the density profile, and only the pressure profile changes. 
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Chapter 6 
Conclusions 



As the range of problems tackled is so diverse, it is hard to think of a concise way of 
summarising all of them. Indeed, to a large extent much of the work has already has 
summarised them at the end of each section. Instead of explaining the problems again 
and how far through completion they are, I would like to leave with list of things that 
I believe are useful to understand the overall concept of the thesis. 

This thesis has been written with the goal of making it accessible to people with 
a basic background in general relativity, especially in coordinate conditions, perfect 
fluid spheres in general relativity. 

6.1 The main analysis: Structure of the thesis 

• The thesis summarized the basic concept of special relativity to make general 
relativity easier to understand. General relativity is the relativistic theory of 
gravitation. It is required from any new physical theory to satisfy the so called 
principle of correspondence. 

• Coordinate conditions are important for the long time scale simulations of rel- 
ativistic systems. We demonstrated the importance of coordinate conditions in 
general relativity and proposed a method of looking for nice solutions of the 
Einstein equations. In addition, a number of examples are given to illustrate the 
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method. 

• We introduced "algorithmic" techniques that permit one to generate large classess 
of perfect fluid spheres in a purely mechanical way. Furthermore, we extended 
these algorithmic ideas, by proving several solution-generating theorems of vary- 
ing levels of complexity. We have developed several transformation theorems 
that map perfect fluid spheres to perfect fluid spheres using Schwarzschild co- 
ordinates. Furthermore, we used these transformations as a basis for classifying 
different types of perfect fluid sphere solutions. 

• In chapter 4, we have developed several transformation theorems that map per- 
fect fluid spheres to perfect fluid sphere using other coordinates. 

• We rephased Theorems 1, 2, 3 and 4 in chapter 3 directly in terms of the TOV 
equation. 

6.2 Further interesting issues 

I would like to leave with the list of things that I believe are interesting to continue to 
analyze. In addition, all the lists below are further issues that might be interesting to 
looking at. 

• For a Riemannian geometry in (2+1) dimensions, we can calculate the Ricci 
Scalar, the Ricci tensor and the Riemann tensor. However, all results are too 
difficult to analyze in detail. This is one of the physically interesting cases, and 
this is a case that I would like to examine in more detail. 

• Similar for a Riemannian geometry in (3+1) dimensions. This is also one of 
the physical interesting case, and this is the case that will be most difficult to 
analyze in detail. This is a challenging case that I would like to work through in 
more detail. 

• We analysed a generic spacetime using both ADM and reverse- ADM (Kaluza- 
Klein) decompositions. However, we found that for the general ADM case co- 
ordinate techniques are too complicated to contemplate (page 514 of fZ|). The 



144 



inverse ADM decomposition can calculated the Christoffel symbols from equa- 
tion (|2.178|) . This is so interesting to find the techniques for contemplation. 

• The other things that would be interesting is consider other coordinate systems: 
Generalize the theorems which we derived for Schwarzschild coordinates or co- 
ordinates in Chapter 4 to a number of other coordinate systems. This may allow 
us to develop extra relations between the known solutions. 

• For Theorem 8 in Isotropic coordinate, it would be interesting to see what would 
happen when we feed this theorem a Minkowski space seed. In addition, this is 
interesting to find simple ways of generating interesting solutions. 

• For Theorem 9 in Gaussian polar coordinate, Theorem 10 in Buchdahl coordi- 
nates, Theorem 11 in Synge isothermal coordinates, and Theorem 12 in Expo- 
nential coordinates, it would be interesting to see what would happen when we 
feed this theorem a Minkowski space seed. In addition, this is intersting to find 
the way that we can do to R(r) to make its ODE look nicer. 

All the above suggestions would be interesting and feasible, although some of them 
would be more tedious to extend to standard results in general relativity. However, 
this thesis would provide a platform for better understanding coordinate conditions 
and the perfect fluid sphere in general relativity. This thesis developed a way of looking 
for nice new solutions of the Einstein equations. Furthermore, the goal of this thesis 
was to classify the different types of perfect fluid sphere solutions. This thesis has 
hopefully served as a brief introduction the "Algorithms" ideas by showing how to 
prove several solution-generating theorems of varying levels of complexity. Finally, we 
illustrated the formalism to develop extra relations between the known solutions for 
Schwarzschild coordinates and other coordinates. The generating theorems in perfect 
fluid spheres lead to additional useful "exact solutions" . 



145 



146 



Bibliography 



[1] R. Adler, M. Bazin, and M. Schiffer. Introduction to general relativity. McGraw- 
Hill Book company, 1965. 

[2] Various authors. Special relativity. 2005. [http://en.wikipedia.org/wiki/Special 
relativity] . 

[3] S. Berger, R. Hojman, and J. Santamarina. General exact solutions of Ein- 
stein equations for static perfect fluids with spherical symmetry. J. Math. Phys., 
28:2649, 1987. 

[4] H. Bondi. Spherically symmetrical models in general relativity. Mon. Not. Roy. 
Astron. Soc, 107:410, 1947. 

[5] H. Bondi. Massive spheres in general relativity. Mon. Not. Roy. Astron. Soc., 
282:303-317, 1964. 

[6] P. Boonserm, M. Visser, and S. Weinfurtner. Generating perfect fluid spheres in 
general relativity. Phys. Rev., D71 (124037), 2005. [arXiv:gr-qc/0503007]. 

[7] H. A. Buchdahl. General relativistic fluid spheres. Phys. Rev., 116:1027-1034, 
1959. 

[8] H. A. Buchdahl. General relativistic fluid spheres II: General inequalities for 
regular spheres. Ap. J., 146:275-281, 1966. 



147 



[9] M. S. R. Delgaty and K. Lake. Physical acceptability of isolated, static, spher- 
ically symmetric, perfect fluid solutions of einstein's equations. Comput. Phys. 
Commun., 115:395, 1998. [arXiv:gr-qc/9809013]. 

[10] P. A. M. Dirac. General theory of relativity. Princeton University Press, 2 edition, 
1975. 

[11] M. R. Finch and J. E. F. Skea. A review of the relativistic static fluid sphere. 1998, 
unpublished, [http:/ /www. dft.if.uerj.br/usuarios/JimSkea/papers/pfrev.ps]. 

[12] V. Fock. The Theory of Space, Time, and Gravitation. Pergamon, New York, 
1964. 

[13] K. Lake. All static spherically symmetric perfect fluid solutions of Einstein's 
equations. Phys. Rev., D67: 104015, 2003. [arXiv:gr-qc/0209104]. 

[14] L. D. Landau and E. M. Lifshitz. The classical theory of fields, volume 2. 
Butterworth-Heinemann, Oxford, 4 edition, 1975. 

[15] D. Martin and M. Visser. Bounds on the interior geometry and pressure profile of 
static fluid spheres. Class. Quant. Grav., 20:3699, 2003. [arXiv:gr-qc/0306038]. 

[16] D. Martin and M. Visser. Algorithmic construction of static perfect fluid spheres. 
Phys. Rev., D69: 104028, 2004. [arXiv:gr-qc/0306109]. 

[17] C. W. Misner, K. S. Thorne, and J. A. Wheeler. Gravitation. (W. H. Freeman 
and Company, San Francisco), 1973. 

[18] S. Rahman and M. Visser. Spacetime geometry of static fluid spheres. Class. 
Quant. Grav., 19:935, 2002. [arXiv:gr-qc/0 103065]. 

[19] R.D'Inverno. Introducing Einstein's Relativity. (Clarendon Press, Oxford), 2002. 

[20] H. Stephani, D. Kramer, M. MacCallum, C. Hoenselaers, and E. Herlt. Exact 
solutions of Einstein's field equations. (Cambridge University Press), 2003. 

[21] J. L. Synge. Relativity: the General Theory. North-Holland, Amsterdam, 1964. 



148 



[22] M. Visser. Math 465: Notes on general relativity. Victoria University, 2005. 
[http:/ / www.mcs.vuw.ac.nz/courses/MATH465/2005T2/Lecture-Notes/]. 

[23] M. Visser, C. Barcelo, and S. Liberati. Analogue models of and for gravity. Gen. 
Rel. Grav., 34:20, 2002. [arXiv:gr-qc/0111111]. 

[24] R. M. Wald. General Relativity. The University of Chicago Press, 1984. 

[25] M. Wyman. Radially symmetric distributions of matter. Phys. Rev., D75:1930, 
1949. 



